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The  objective  of  this  study  was  to  carry  out  a tradeoff  analysis  among  systems  parameters 
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OBJECTIVE 


Investigate  the  potential  and  suitability  of  optical  communication  links  with  anti- 
jam and  anti-intercept  capabilities  for  task  force/group  operations  at  greater  than  line-of-sight 
ranges.  Specifically,  perform  a tradeoff  analysis  among  systems  parameters  for  two  types  of 
extended  line-of-sight  optical  communications  links;  optical  forward  scatter  from  the  atmos- 
phere and  remotely  piloted  vehicle  relay  links. 


CONCLUSIONS 

Based  upon  the  analytical  model  developed  for  over-the-horizon  optical  propagation 
and  a systems  tradeoff  analysis,  the  following  conclusions  were  drawn:  the  operating  data 
links  can  use  the  same  shipboard  system;  significant  performance  advantages  can  be  achieved 
by  exploiting  decreasing  aerosol  concentrations  with  elevation;  the  primary  factor  which 
determines  communication-link  availability  at  a given  range  anu  data  rate  is  meteoro- 
logical visibility  and,  based  upon  statistical  studies  of  visibility,  communication  ranges  for 
bit  rates  of  2400  bits  per  second  and  for  rates  of  75  bits  per  second  were  determined  for 
availabilities  of  both  85  and  70  percent.  The  communication  ranges  are  considerably  greater 
at  night  or  when  low  clouds  occur  to  provide  a scattering  layer. 


RECOMMENDATIONS 

The  analytical  model  developed  in  this  study  should  be  experimentally  tested  to 
assess  its  validity.  Simultaneous  pulsed  laser  propagation  measurements  at  two  wavelengths 
(1.06  and  0.53  pm)  should  be  conducted.  The  performance  gain  achieved  by  the  use  of 
elevated  beams  and  exploitation  of  clouds  should  be  evaluated.  Fading  rate  and  other  data 
should  be  obtained  to  allow  evaluation  of  time,  spatial,  and  receive  angle  diversity  receiver 
techniques. 

It  is  strongly  recommended  that  experiments  be  conducted  to  verify  the  meteoro- 
logical effects  upon  propagation  predicted  by  the  models.  An  integrated  and  systematic 
capability  for  measuring  and  predicting  the  meteorological  conditions  predominating  over 
the  path  should  be  developed. 
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SECTION  1.  INTRODUCTION 


The  wide  dispersal  of  elements  of  a modem  task  force  requires  coordinated  intra- 
element,  real-time  communications.  At  present,  the  Fleet  cannot  operate  with  positive 
Command  Control,  and  Communications  (C^)  without  making  itself  vulnerable  to  enemy 
detection  and  countermeasures.  The  Naval  Telecommunications  Architect  has  stipulated 
that  the  post- 1985  Naval  Telecommunications  System  must  have  survivability  equal  to  the 
forces  served  in  the  face  of  both  physical  and  electromagnetic  attack.  It  is  planned  that 
this  goal  will  be  achieved  through  media  diversification  and  anti-jam,  low-probability-of- 
intercept  ( AJ/LPI ) capabilities.  For  intra  task  force  (group)  netting,  hf  systems  with  AJ/ 

LPI  capabilities  have  been  specified  as  the  prime  service,  but  hf  will  require  considerable 
augmentation  to  accomplish  the  full  capability  needed  for  the  Naval  Tactical  System  (NTS). 

Optical  systems  are  strong  candidates  for  this  augmentation  role  for  a number  of 
reasons.  These  include  reduced  susceptibility  to  jamming,  intercept,  spoofing,  and  direction 
finding;  a reduction  in  spectrum  crowding;  and  potentially  high  data  rates.  Optical  links 
also  have  potential  limitations  including  meteorological  effects,  multiplatform  transmissions, 
and  range  limitations  imposed  by  the  horizon. 

The  true  potential  of  electro-optic  communications  systems  needs  to  be  better 
defined  in  order  to  determine  its  proper  future  (1985  era)  use  in  the  NTS  intra-task  group 
application.  This  use  is  to  be  scoped  in  relation  to  other  candidate  future  systems  such  as 
hf  and  JTIDS  (Joint  Tactical  Information  Data  System).  Systems  concept  engineering  is 
needed,  as  is  more  extensive  data  on  performance  capabilities  and  limitations.  The  candidate 
concepts  are  being  explored  in  a variety  of  projects  within  the  electro-optical  communica- 
tions program  at  NELC.  The  particular  aspect  addressed  by  “Extended  Line-of-Sight  Opti- 
cal Communications  System  Exploratory  Development”  is  to  examine  and  demonstrate 
the  feasibility  for  maintaining  optical  links  to  distances  well  beyond  the  horizon.  This  re- 
port represents  a theoretical  analysis  of  the  concept  of  optical  communications  to  beyond 
line-of-sight  ranges. 


OBJECTIVES 

The  primary  objective  of  this  study  was  to  carry  out  a tradeoff  analysis  among 
system  parameters  for  two  types  of  beyond  line-of-sight  optical  communication  links: 
Forward  Scatter  Channel  and  Remotely  Piloted  Vehicle  (RPV)  Relay.  The  general  appli- 
cation of  such  communications  links  is  for  the  control  of  dispersed  task  force  units. 

The  tradeoff  study  for  the  forward-scatter  channel  originally  intended  to  use 
available  scatter  models.  However,  discrepancies  were  noted  between  the  models  and  data. 
This  type  of  discrepancy  was  particularly  apparent  in  the  observed  small  spot  size  on  the 
horizon  compared  to  the  large  diffuse  illumination  predicted  by  many  models.1  Thus,  an 
added  objective  was  to  generate  a forward-scatter  model  which  would  allow  for  realistic 
engineering  estimates  of  scattered  power  at  ranges  beyond  the  horizon. 


1 Naval  Research  Laboratory  Report  6152,  Experimental  Observations  of  Forward  Scattering  of  Light  in 
the  htwer  Atmosphere,  by  JA  Curcio  and  LF  Drumrncter,  Jr,  30  September  1964 
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SCOPE 


The  tradeoff  study  was  broadly  divided  into  two  parts:  analysis  of  the  optical  chan- 
nel, or  link  analysis,  and  equipment-dependent  factors.  The  link  analysis,  embodying  con- 
sideration of  background  noise  and  transmission  loss,  is  more  complex  analytically  and  is 
presented  in  a form  such  that  it  usually  need  not  be  repeated  for  a change  in  system  para- 
meters. The  results  of  the  link  analysis  are  presented  as  signal-to-noise  ratios  (S/N)  for  a 
group  of  nominal  system  parameters.  The  resulting  S/N  for  a change  in  a parameter,  eg 
radiated  power,  can  usually  be  rapidly  determined. 

A forward  scatter  propagation  model  which  appears  consistent  with  previous  field 
data  was  generated.2  The  model  also  revealed  novel  strategies  for  reducing  path  loss,  eg, 
the  use  of  vertical  fan  beams  and  the  use  of  an  optimum  elevation  angle  for  pencil  beams. 
This  model  was  used  for  the  scatter  channel  path  loss  computations.  The  scatter  link  analy- 
sis predicted  that  the  near  infrared  (IR)  wavelength  would  provide  the  best  performance. 
Consideration  of  state  of  the  art  sources  and  detectors  indicates  that  1 .06  microns  is  the  pre- 
ferred wavelength.  Therefore,  the  most  complete  data  for  both  communication  links  are 
presented  for  that  wavelength.  Complete  link  and  equipment  analysis  at  all  wavelengths 
would  be  prohibitive.  However,  since  the  final  definition  of  an  over-the-horizon  system  or 
systems  is  likely  to  be  an  iterative  analytical  and  experimental  process,  computer  programs 
were  generated  to  provide  additional  data  beyond  those  presented  in  this  report.  The  pro- 
grams were  designed  tor  ease  of  altering  channel  parameters. 

Characteristics  of  various  optical  system  components  and  possible  tradeoffs  are  dis- 
cussed. Emphasis  is  on  the  near  IR  but  some  discussion  is  included  for  ultraviolet  through 
IR  wavelengths.  A discussion  of  modem  techniques  and  required  S/N  is  included. 

An  attempt  has  been  made  to  include  information  which  might  be  of  use  as  design 
concepts  change,  since  this  is  apt  to  happen  in  any  new  technology.  Thus,  besides  discuss- 
ing the  system  concepts  presented,  additional  information  of  the  design  handbook  type 
has  been  included. 


PLAN  OF  REPORT 

Tiie  general  system  performance  requirements  are  given  in  section  2.  Section  3 dis- 
cusses the  concepts  for  the  two  over-the-horizon  links  and  summarizes  the  system  parameters 
and  expected  performance.  Section  4 discusses  signal  modems  and  S/N  requirements  for 
system  acquisition.  Communication  data  rates  are  also  considered. 

Transmission  loss  data  tor  both  link  types  are  given  in  section  5.  This  is  combined 
with  background  radiation  to  provide  predicted  S/N  for  nominal  system  parameters.  The 
data  presented  in  section  5 are  based  on  models  and  computer  programs  described  in  detail 
in  the  appendixes. 

The  characteristics  of  various  system  components  such  as  lasers,  optical  filters,  and 
detectors  are  given  in  section  6.  Tradeoffs  between  the  component  parameters  are  discussed 
briefly  in  section  7.  Covert  and  antijamming  factors  are  briefly  considered  in  section  8. 
Section  l>  presents  conclusions  and  recommendations;  the  more  detailed  summary  is  pre- 
sented in  section  3. 


2 Megatek  Report  R2005-039-IF-1,  Propagation  Model  for  a Laser-Type  Bevond-thc-Horizon  Communica- 
tions Link , by  PH  Levine,  15  December  1975 
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SECTION  2.  GENERAL  REQUIREMENTS 

The  main  use  of  an  over-the-horizon  optical  link  would  be  to  communicate  between 
naval  task  groups.  Communication  within  a group  would  generally  be  on  a line-of-sight 
basis,  although  slightly  extended  line-of-sight  might  be  considered  within  a group  for  in- 
creased dispersal.  The  range  of  interest  for  the  links  is  25  to  1 50  nautical  miles. 

The  present  Navy  Tactical  Data  System  (NTDS)  typically  uses  a rate  of  2400  bits 
per  second.  This,  then,  becomes  the  desired  information  rate  for  the  system.  There  are 
other  tactical  requirements  which  can  be  satisfied  with  a rate  of  75  bits  per  second.  A high 
degree  of  jamming  and  intercept  protection  is  expected,  and  a general  requirement  is  for 
the  system  to  be  much  less  susceptible  to  interference  and  intercept  than  high-frequency 
radio  systems.  In  summary,  the  requirements  are: 

Range:  25  to  1 50  nautical  miles 

Data  Rate:  75  bits  per  second  (minimum) 

2400  bits  per  second  (desired) 

Covertness:  High  anti-jamming  margin 

Low  probability  of  intercept 


SECTION  3.  SYSTEM  CONCEPTS  AND  PERFORMANCE  SUMMARY 

This  section  discusses  the  concepts  for  two  over-the-horizon  (OTH)  optical  commu- 
nications links.  Modulation  and  detection  considerations  are  summarized.  Link  acquisition 
and  communication  data  rates  are  discussed.  Finally,  systems  composed  of  state-of-the-art 
components  and  subsystems  are  defined,  and  performance  predictions  are  given  or  each 
type  of  link. 


SYSTEM  TYPES 
LINKS 

Two  different  types  of  OTH  optical  communications  systems  are  consul*  red.  They 
are:  Forward-Scatter  Channel  and  Active  RPV  Relay.  The  forward-scatter  systc  m makes 
use  of  scattering  from  atmospheric  “haze”  and  aerosols  and  does  not  depend  u > hi  the 
presence  of  clouds.  If  clouds  were  present,  scatter  from  them  could  often  be  i;  d,  depending 
upon  their  altitude. 

An  active  relay  would  receive  an  optical  line-of-sight  signal  from  one  sr  i.  would 
demodulate  the  signal,  and  would  use  it  to  modulate  a laser  on  the  RPV.  Thi  in  optical 
receiver  and  transmitter  and  a pointing  system  would  be  located  on  the  RPV  \ power 
source  would  be  required  for  the  transmitter,  receiver,  and  pointing  system 

Both  links  would  use  high-power  shipboard  transmitters  with  pointin  nd  tracking 
systems.  The  receiver  would  be  boresiglited  to  the  transmitter  for  use  with  ' relay  link 
but  might  be  capable  of  pointing  and  tracking  independently  for  the  scattc  annel. 


A particular  Fleet  application  might  use  any  one  or  a combination  of  links.  For 
example,  for  a certain  ship  spacing,  a scatter  channel  might  be  used  at  night  and  during  high- 
visibihty  daytime  conditions.  If  the  visibility  began  to  decrease,  as  indicated  by  a decrease  in 
signal  level,  an  RPV  could  be  launched. 


HETERODYNE  VERSUS  DIRECT  DETECTION 

Heterodyne,  or  coherent,  optical  detection  systems  have  been  developed  at  10.6/um 
which  approach  theoretical  sensitivity  limits  very  closely.  At  10.6  jum,  a heterodyne  system 
ls  about  40  dB  more  sensitive  than  a direct-detection  receiver  limited  by  300°K  background 
ources  which  peak  out  at  10  gin.  However,  at  wavelengths  shorter  than  about  1 . 1 gim,  de- 
tec  tors  are  available  with  low  noise  internal  gain  which  give  dramatic  improvement  in  direct 
detection,  as  compared  to  available  direct  detectors  at  10  /am.  In  addition,  both  practical 
problems  and  inherent  limitations  in  coherent  detection,  especially  relative  to  a scatter  chan- 
nel. are  compounded  at  shorter  wavelengths.  For  the  reasons  discussed  here,  a direct  detec- 
tion system  is  assumed  (wavelength  less  than  1.1  /am). 

The  receiving  aperture  and  the  field  of  view  in  a heterodyne  receiver  are  not  inde- 
pendent. but  arc  related  by  the  approximate  relation.3 

A f2r  = 

or  A-j0^=x2 

4 r 

where 

A = receiving  aperture, 

= solid  angle  field  of  view. 

Of  = plane  angle  field  of  view,  and 

X = wavelength. 

Note  that  for  a fixed  field  of  view,  the  ap.  area  varies  as  the  square  of  the  wavelength, 
for  a field  ot  view  (FOV)  of  1 degree,  such  as  might  be  used  in  a scatter  channel,  the  max- 
imum aperture  at  I gm  is  about  4 X 109  square  metres.  Even  for  a FOV  of  0 1 
degree,  which  might  be  used  in  a line-of-sight  (LOS)  link,  the  aperture  area  is  about  4 X 10~7 
square  metres.  A direct-detection  receiver  does  not  have  to  satisfy  the  coherence  require- 
ment and  may  have  a large  aperture. 

Spatia1  coherence  is  degraded  due  to  atmospheric  turbulence ; 4 the  coherence  diameter 
vanes  asX  IS.  I bus,  turbulence  is  more  of  a problem  at  short  wavelengths.  The  coherence 
properties  of  the  scatter-channel  field  have  not  been  analyzed  but  it  is  expected  that  coher- 
ence would  be  lower  at  shorter  wavelengths. 

Heterodyne  systems  are  more  complex,  in  that  a very  stable  local  oscillator  is  re- 
quired (this  may  be  the  unmodulated  transmitter  source  for  a homodyne  system).  There 
are  also  stringent  alignment  tolerances  required  to  keep  the  signal  and  local  oscillator  in 
phase  over  the  detector  surface. 


Sicgman,  AE,  The  Antenna  Properties  of  Optical  Heterodyne  Receivers.  Proceedings  of  the  II 1 1 v S4 
p 1350  to  1356,  1966 

Pratt.  WK,  Laser  Communications  Systems.  John  Wiley  & Sons,  1969 
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Based  upon  the  above  factors  and  the  fact  that  the  path  loss  is  excessive  at  10.6  Mm, 
a direct-detection  system  is  assumed  throughout  this  report.  If  the  atmospheric  window 
at  3.8  Mm  becomes  available  through  the  development  of  new  lasers  in  the  future,  hetero- 
dyne receivers  should  be  considered. 


MODULATION 

For  a direct-detection  system,  short,  high-power  pulses  enhance  the  ability  of  a re- 
ceiver to  discriminate  against  a fixed  background  by  allowing  time  gating.  The  signal-to- 
noise  ratio  (S/N)  expression  given  in  section  4 shows  the  advantage  quantitatively.  As  will 
be  discussed  in  section  5,  shortening  the  pulse  much  below  the  time  dispersion  in  a scatter 
channel  does  not  continue  to  improve  the  S/N.  However,  on  a LOS  link,  it  appears  that 
the  pulse  duration  (constant  energy)  should  be  as  short  as  practical. 

The  use  of  pulse  modulation  has  a number  of  potentially  practical  features. 

The  space  between  pulses  can  be  monitored  for  interference  on  any  of  the  link  types. 

The  pulses  can  be  used  in  a light-detection-and-ranging  (L1DAR)  system  for  possible  de- 
tection of  approaching  missiles,  or  for  monitoring  the  location  of  an  RPV. 

Some  types  of  pulse  modulation  and  expected  data  rates  are  discussed  in  section  4. 
Pulse-modulation  systems  have  generally  been  assumed  throughout  the  report,  although 
many  results  (eg.  path  loss)  apply  independent  of  equipment  types.  Pulse  rates  of  about 
200  pulses  per  second  (pps)  are  required  for  a data  rate  of  2400  bits  per  second. 

Polarization  and  intensity  modulation  appear  possible  for  the  relay  link.  Phase  mod- 
ulation should  also  be  possible  but  could  not  be  detected  with  a non-heterodyne  system. 
Pulse-rate  requirements  for  a pulsed  interrogation  beam  are  discussed  in  section  4.  A mini- 
mum 'ate  of  480  pps  would  be  required  for  a data  rate  of  2400  bits  per  second  assuming  a 
1 >-dB  S/N  ( 5 bits  per  pulse). 


ACQUISITION 

FORWARD-SCATTER  LINK 

An  acquisition  subsystem  similar  to  that  used  in  the  OCCULT  system  is  assumed  for 
the  forward-scatter  link.  The  pointing  system  on  all  ships  would  rotate  continuously  and 
synchronously  through  360  degrees  of  azimuth.  A ship  initiating  acquisition  would  flip  the 
transmitting  beam  180  degrees  in  azimuth  and  begin  transmitting  with  a 4 degree  vertical 
fan  beam.  Because  of  the  synchronous  rotation,  the  receiver  FOV  of  another  ship  would  be 
pointing  at  the  transmitting  ship  at  the  time  the  beam  swept  the  ship. 

The  single  scattering  out  of  the  beam  will  appear  to  the  receiver  as  an  approximately 
vertical  line  source  (see  appendix  E).  The  forward  scatter  is  peaked  in  the  forward  direction 
but  has  some  angular  width,  so  that  the  line  (starting  at  the  horizon)  can  be  “seen”  over 
several  degrees  of  azimuth.  Assume  a horizontal  FOV,  0|„  for  the  receiver.  The  FOV  of  a 
single  receiver  will  be  limited  to  about  4 degrees  by  the  interference  filter  constraints  (assume 
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1A  width;  see  section  6).  The  receiver  would  see  the  line  source  a fraction,  0.50^/360°,  of  a 
rotation  of  the  acquisition  system.  The  number  of  pulses,  Np,  entering  the  receiver  would 
be: 


N = — — 

P 720A’ 

where 

0h  = horizontal  field  of  view  in  degrees, 

A = Azimuthal  scan  rate  - revolutions  per  second,  and 
R = Pulse  repetition  rate  - pulses  per  second. 

As  an  example  to  indicate  approximate  scan  rates,  consider  a 4 degree  horizontal 
FOV  (0^  = 4°)  and  require  4 pulses  to  occur  within  the  receiver  FOV.  A range  of  acquisition- 
response  criteria  could  then  be  used.  For  normal  operations,  a single  detected  pulse  might 
initiate  a response,  while  a “quiet”  ship  might  require  four  pulses,  or  even  more  than  one 
group  of  four  pulses,  at  the  proper  time  interval,  before  triggering  a response.  Reference 
curves  useful  for  determining  detection  probability  and  false  alarm  probability  tradeoffs 
are  given  in  section  4.  When  more  than  one  pulse  is  detected,  the  “center  of  gravity” 
direction  would  be  used  for  the  transmitted  response.  If  the  laser  was  operated  at  the  200 
pps  nominal  rate  required  for  communication  at  2400  bits  per  second,  the  azimuthal  scan 
period  would  be  3.6  seconds  per  revolution.  Normally,  up  to  2 revolutions  (7.2  seconds) 
would  be  required  for  azimuthal  acquisition. 

After  azimuthal  acquisition,  the  first  ship  could  switch  from  a fan  beam  to  a pencil 
beam  and  scan  the  transmit  beam  in  elevation.  The  received  signal  level  at  the  second  ship 
would  be  fed  back  to  the  first  ship  to  enable  it  to  select  the  optimum  beam  elevation.  The 
second  ship  would  then  go  through  the  same  procedure  to  complete  acquisition. 

The  S/N  required  for  acquisition  depends  on  the  identification  friend  or  foe  (IFF) 
code  and  probability  requirements.  As  discussed  in  section  4,  1 5-dB  S/N  should  cover  most 
IFF  requirements.  Many  of  the  S/N  curves  in  this  report  are  for  a FOV  of  1 degree,  while 
the  example  assumes  a FOV  of  4 degrees.  Increasing  the  vertical  FOV  increases  both  the 
received  energy  and  the  received  background  noise.  As  is  discussed  in  section  5,  the  S/N 
increases  with  increased  vertical  FOV  until  the  time  dispersion  is  about  equal  to  the  pulse 
duration  and  there  remains  constant  if  the  receiver  electrical  bandwidth  is  adjusted  correctly. 
For  typical  pulse  lengths  (0.02  Msec)  at  1.06  micrometres,  the  S/N  does  not  generally  im- 
prove for  an  elevation  FOV  greater  than  1 degree.  Increasing  the  horizontal  FOV  from 
1 to  4 degrees  increases  the  background  noise  by  6 dB.  The  energy  for  a single  pulse  is  not 
increased  as  it  was  for  a vertical  FOV.  Four  pulses  are  received  and  filter  techniques,  such 
as  summing  the  pulses  by  use  of  a tapped  delay  line,  can  produce  essentially  the  same  S/N 
as  for  a FOV  of  I degree.  However,  it  would  probably  be  easier  to  use  a FOV  of  1 degree 
when  a single-pulse  detection  criterion  is  used.  When  a 4-pulse  IFF  code  with  a 1 5-dB  S/N 
requirement  is  used,  the  2 1 -dB  S/N  level  on  a FOV  curve  of  1 degree  would  be  used  for 
coverage  prediction. 

The  S/N  and  resulting  IFF  detection  can  be  enhanced  at  the  expense  of  acquisition 
time.  Rather  than  the  example  given  above,  assume  a FOV  of  1 degree  with  a resulting 
6-dB  lower  noise  level.  If  the  azimuthal  scan  rate  is  reduced  by  a factor  of  4,  4 pulses  will 
occur  within  the  receiver  FOV.  Increasing  the  pulse-repetition  rate  does  not  improve  the 
detection  since  the  laser  has  an  average  power  limitation  so  that  the  peak  power  is  reduced. 


If  approximate  ship  locations  are  known  from  a previous  acquisition  or  by  other 
means,  the  acquisition  transmission  could  be  applied  to  a restricted  azimuthual  region.  This 
would  reduce  the  probability  of  intercept.  If  the  pointing  systems  on  both  ships  were  re- 
stricted in  azimuth,  the  acquisition  time  could  be  reduced. 


RPV  RELAY  LINK 

The  acquisition  technique  for  an  RPV  relay  link  would  be  very  similar  to  that 
described  for  the  scatter  link.  The  launch  ship  would  activate  an  azimuthal  scan  with  a 
vertical  fan  beam  on  the  RPV  when  it  was  at  a sufficiently  high  altitude.  Figure  25  in  sec- 
tion 5 shows  that  a beam  divergence  of  1.8  degrees  would  provide  a range  capability  of  a 
full  ( 1 50  nautical  miles  (RPV  at  midpoint)  for  a visibility  of  20  kilometres.  The  fan  beam 
would  be  approximately  7 degrees  by  0.25  degree.  Requirements  of  from  1 to  4 pulses  for 
acquisition  result  in  4 by  360  to  16  by  360  pulses  per  azimuthal  scan.  Use  of  a rate  of 
200  pps  gives  a scan  time  of  between  7.2  and  29  seconds.  If  the  launch  ship  were  to  know 
the  approximate  location  of  the  other  ship,  the  RPV  transmissions  could  be  restricted  in 
azimuth  to  minimize  the  probability  of  intercept  and  to  conserve  energy. 


COMMUNICATION  RATES 

The  communication  rate  for  a pulse-position  modulation  system  is  discussed  in 
section  4.  It  is  shown  that  a rate  of  1 2 bits  per  pulse  is  not  far  from  the  theoretical 
information  rate  when  pulse-time  dispersion  (for  the  scatter  link)  is  taken  into  account. 

A data  rate  of  2400  bits  per  second  requires  a 200-pps  rate  using  1 2-bit  words. 
For  a data  rate  of  75  bits  per  second.  1 2-bit  words  would  still  be  used  and  the  pulse  rate 
would  be  reduced.  The  reduction  in  pulse  rate  would  be  converted  into  a higher  peak- 
pulse  power  with  a resultant  extension  in  operating  range. 


STATE-OF-THE-ART  SYSTEM  PERFORMANCE 

SYSTEM  PARAMETERS 

The  link  analysis  given  in  section  5 clearly  shows  the  advantage  of  the  near  infrared 
for  both  a torward-scatter  and  a LOS  optical  system  Consideration  of  the  state  of  the  art 
in  laser  sources  and  detectors  shows  that  1 .06  micrometres  is  the  preferred  wavelength. 

In  this  section,  we  consider  systems  operating  at  1.06  micrometres  and  maximize 
all  system  parameters,  consistent  with  available  components  as  discussed  in  section  6.  The 
parameters  tor  the  systems  defined  in  this  section  differ  from  the  nominal  system  parameters 
used  in  section  5.  The  nominal  system  parameters  (section  5)  were  selected  to  be  somewhat 
typical  ol  a large  span  of  wavelengths  and  were  chosen  to  be  powers  of  10  where  reasonable. 

The  system  concept  synthesized  in  this  study  uses  the  same  shipboard  equipment 
for  both  data  links,  but  with  different  operating  parameters,  as  indicated  in  table  1.  Most 
of  the  parameters  are  based  upon  components  listed  in  section  6.  Significant  advances  cur- 
rently being  made  in  the  development  of  plastic  optics,  are  not  discussed  in  this  report. 
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TABLE  1.  1.06  /jm  SHIPBOARD  SYSTEM  PARAMETERS. 


Scatter  Link 

Active  RPV 

Parameter 

Constant 

Acqui- 

sition 

75 

b/s 

2400 

b/s 

i 

Acqui- 

sition 

2400 

b/s 

Transmitter 

Radiated  peak  power  (MW) 

40 

1280 

*** 

40 

40 

40 

Pulse  rate  (pps) 

200 

25 

••• 

200 

200 

200 

Beam  divergence  (mrad) 
Pulse  width  (psec) 

0.015 

1°  beam 
or  4°  fan 

1.5 

1.5 

1.5 

1.5 

Receiver 

Field  of  view  (deg) 

Aperture  area  (m3) 
Optical  bandwidth  (pm) 
Detector  responsivity 

0.2 

10-4 

0.018 

1° 

1° 

1° 

7°  X 
V4°  fan 

2° 

(amp/watt) 

Gain 

Dark  current  (amp) 

* 

10*13 

Electrical  bandwidth  (Hz) 

** 

*• 

** 

7 X 107 

i 

7 X 107 

Notes: 


’Photomultiplier  gain  sufficient  to  give  dark  current  (or  background)  limited  operation. 

••Bandwidth  adjusted  to  be  equal  to  the  inverse  of  the  received  pulse  duration  at  planned  maximum  range. 
•••See  text. 


However,  the  use  of  plastics  is  expected  to  reduce  large-aperture  costs  from  a few  thousand 
to  between  ten  and  a few  hundred  dollars.  A 50.8  cm  diameter  aperture  (0.2  m^)  is  assumed 
for  the  system. 

The  Nd : YAG  laser  is  assumed  to  have  1 20  watts  average  power,  and  a tradeoff  be- 
tween pulse  rate  and  peak  power  is  made  for  different  operating  conditions.  The  tradeoff 
to  1 280  megawatts  (6  pps),  listed  for  a data  rate  of  75  bits  per  second,  may  be  impractical. 
However,  if  it  cannot  be  achieved,  it  is  expected  that  receiver  techniques  which  average  over 
many  pulse  intervals  may  produce  similar  effects  on  S/N. 

A beam  divergence  of  1.5  milliradians  (mrad)  is  listed  for  data  transmission  on  the 
links.  The  divergence  is  not  critical  on  the  scatter  link  as  long  as  it  is  less  than  about  17  mrad 
( 1 degree).  In  the  scatter  link  performance  prediction,  it  is  assumed  that  the  beam  is  pointing 
at  the  optimum  elevation  angle. 

The  shipboard  transmitter  parameters  for  the  active  RPV  link  are  listed  as  being  the 
same  as  for  the  scatter  link.  This  represents  much  more  power  than  required  for  the  uplink. 
The  power  would  normally  be  reduced  to  a value  just  sufficient  for  satisfactory  communica- 
tions, in  order  to  reduce  the  probability  of  intercept.  The  shipboard  receiver  also  uses  a 
vertical  fan-shaped  FOV  for  link  acquisition. 
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The  general  characteristics  of  an  active  RPV  terminal  used  for  this  performance  pre- 
diction are  listed  in  table  2.  The  only  vehicle  parameter  used  explicitly  in  the  performance 
prediction  is  the  altitude.  The  other  parameters,  listed  for  convenience,  give  the  expected 
capability  over  the  next  few  years  (see  section  6 for  a discussion  of  RPVs). 


TABLE  2.  RPV  TERMINAL. 


Parameter 

Value 

Remotely  Piloted  Vehicle 

Altitude 

2.5  km 

Loiter  time 

6 hours 

Payload 

50  pounds 

Transmitter 


Radiated  power 

1 MW 

Pulse  rate 

200  pps 

Beam  divergence: 

narrow  beam 

1 mrad 

wide  beam 

1.8° 

The  ship-to-RPV  link  is  not  explicitly  considered  here  since  the  range  with  fairly 
simple  receivers  exceeds  the  down-link  range.  Two  ditferent  examples  of  beam  divergence 
are  considered  for  the  active  RPV  transmitter.  The  wide-beam  case  (1.8°)  will  provide  less 
coverage  but  will  greatly  simplify  the  pointing  system. 


PERFORMANCE 

Table  3 gives  the  performance  predictions  for  daytime  background  and  2 visibility 
values.  Tlte  contour  maps  in  section  5 show  that  10-km  visibility  is  typically  available 
about  85  percent  of  the  time.  Data  indicate  that  a 20-km  visibility  would  be  available 
about  70  percent  ot  the  time.  The  link  availability  neglects  the  presence  of  clouds  which 
can  either  degrade  or  enhance  a channel  as  discussed  in  section  5. 

A 1 5-dB  S/N  is  taken  as  the  requirement  for  data  communications.  This  S/N  should 
Provide  low  error  rates  but  includes  little  margin  for  signal  fading. 

The  channel  acquisition  is  based  on  acceptance  of  a single  pulse  with  a 1 5-dB  S/N. 

If  IFF  codes  are  used  to  minimize  false  alarms,  the  acquisition  times  will  be  increased  over 
those  listed  in  the  table.  For  the  forward-scatter  link,  the  use  of  a pencil  beam  of  1 degree 
and  a tan  beam  of  4 degrees  are  considered.  The  estimate  for  the  fan  beam  is  conservative 
relative  to  that  of  the  pencil  beam  since  it  assumed  that  only  the  energy  at  the  best  1 -degree 
elevator  segment  is  used.  In  reality,  other  portions  of  the  beam  will  contribute  energy,  but 
then  pulse  time  dispersion  effects  must  be  included.  Other  tradeoffs  between  acquisition 
time  and  range  are  possible,  and  this  area  should  receive  additional  attention. 
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TABLE  3.  PERFORMANCE  PREDICTION.  DAYTIME. 


•Takes  into  account  loss  due  to  dispersion,  effect  negligible  at  other  ranges. 


It  is  noted  that  an  active  relay  channel  provides  the  full  150-nmi  communication 
coverage  with  an  85-percent  availability.  The  scatter  link  gives  a 100-km  communication 
range  with  70-percent  availability.  This  range  can  be  extended  to  165  km  by  reducing  the 
data  rate  to  75  bits  per  second. 

Table  4 shows  that  considerably  better  performance  is  provided  at  night.  The  scat- 
ter link  can  provide  a data  rate  of  75  bits  per  second  over  the  1 50-nmi  range. 
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SECTION  4.  SIGNAL  MODEM  CONSIDERATIONS 

The  purpose  of  this  section  is  to  describe  genera)  considerations  for  IFF  and 
communication  signal  design,  and  to  establish  signal-to-noise  ratio  requirements.  Some 
examples  are  given,  but  the  final  design  and  specification  of  parameters  will  depend  upon 
the  specific  Fleet  application.  As  discussed  earlier,  the  emphasis  is  on  pulse-modulation 
systems.  Information  is  presented  on  IFF  design  to  allow  for  tradeoffs  between  false- 
alarm  and  detection  probabilities.  The  effect  on  S/N  of  several  system  parameters  are 
considered. 


NOISE  STATISTICS  AND  SIGNAL-TO-NOISE  RATIO 

Photodetection  is  most  correctly  analyzed  as  a Laquerre  counting  process.5  However, 
under  most  conditions  of  interest,  Poisson  statistics  apply  and  the  Poisson  detection  process 
has  been  analyzed  in  some  detail.6 

Following  Karp,  the  average  number  of  electrons  produced  by  a received  optical 
pulse  of  peak  power  Pc  and  duration  T is  given  by 

KS  = J?PcT/h», 

where  77  is  the  photodetector  quantum  efficiency,  v is  the  optical  Irequency,  and  h is  Planck’s 
constant.  Similarly,  the  number  of  electrons  produced  by  background  radiation  in  a time  T 
is  given  by 


Kn  = r?PbT/hv, 

where  P^  is  the  received  background  average  power.  Note  that  Kn  is  proportional  to  T for  a 
constant  background,  so  that  the  use  of  short  pulses  and  corresponding  short  signal  intervals 
decreases  the  amount  of  interfering  radiation. 

The  number  of  electrons  “counted”  in  a time  slot  that  contains  the  signal,  is  a 
Poisson-distributed  random  variable  with  mean  Kn  + Kg,  and  the  number  of  electrons 
counted  in  a slot  containing  only  noise  is  a Poisson  random  variable  with  mean  K\.  The 
ratio  of  the  square  of  the  electron  count  (ie  square  of  detector  current),  when  no  noise  is 
present,  to  the  variance  of  the  count,  when  noise  is  present,  has  many  of  the  characteristics 
of  a S/N  and  is  usually  referred  to  in  this  way.  For  Poisson  statistics  this  becomes. 

S/N  = K“  / (Ks  + KN). 
s 

The  error  probability  for  Poisson  statistics  depends  upon  both  signal  and  noise  energy, 
and  not  upon  merely  their  ratio.  An  example  of  an  error  probability  curve  is  shown  in  figure 
I.  It  is  noted  that  the  error  probability  approaches  the  familiar  Gaussian  case  for  a large 
background  count.  It  appears  that  the  familiar  Gaussian  statistics  can  be  used  for  daytime 

5 G»8li»nh,  RM,  “The  Effect  of  Timing  Errors  in  Optical  Digital  Systems.”  IEEE  Transactions  on  Communi- 
cation Technology,  COM-20,  p 87  to  93,  1972 

‘Gagliardi,  RM  and  S,  Karp.  "Mary  Poisson  Detection  and  Optical  Communications,"  IEEE  Transactions 
on  Communication  Technology , COM- 1 7,  p 208  to  216,  1969 
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background  conditions  while  the  more  precise  Poisson  statistics  would  apply  at  night.  Since 
a system  would  generally  be  designed  for  worst-case  conditions,  the  emphasis  is  on  Gaussian 
statistics. 


The  S/N  expression  can  be  cast  in  more  useful  engineering  terms  by  noting  that  the 
signal  photocurrent.  Ic,  is  the  product  of  the  electron  charge,  q,  and  the  number  of  photons 
per  unit  time.  The  current  may  also  be  expressed  as  the  product  of  the  detector  respon- 
sivity,  Sd,  (ampere  per  watt)  and  the  incident  power  or. 


Ic  = (nq/h^)Pc  = SdPc. 


A similar  expression  holds  for  the  background  current.  Consider  a photodetector  with  inter- 
nal current  gain  G,  feeding  an  amplifier  presenting  a load,  Rl-  It  can  be  shown5  that 


S/N 


sa  p;g: 


2qB  (Sd(Pb  + Pc)  + Id)  G2  + 4kTB/RL 


where 


Sd  is  the  detector  responsivity  (ampere/watt), 
q is  the  electronic  charge  (coloumbs), 

Pc  is  the  receiver  carrier  power  (watts), 

B is  the  electrical  bandwidth  of  the  system  (hertz), 

Pb  is  the  background  power  incident  on  the  detector  (watts), 
Id  is  the  detector  dark  current  (amperes), 
k is  Boltzmann’s  constant  (joules/K), 

T is  the  thermal  noise  temperature  (K),  and 
Rl  is  the  load  resistance  (ohms). 
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For  practical  amplifiers,  T will  be  greater  than  the  ambient  temperature.  This  expression, 
neglecting  the  thermal  noise  term  (second  term  in  the  denominator),  is  used  in  the  link 
analysis  given  in  section  5.  For  sufficiently  high  gain  G,  the  S/N  is  independent  of  the 
thermal  noise. 

Practical  gain  mechanisms  include  current  gain  within  photomultipliers  and  ava- 
lanche photodiodes  and  heterodyne  detection.  As  discussed  earlier,  a heterodyne  system 
is  not  considered  a candidate  for  a scatter  link  because  of  the  operating  wavelength  and 
possible  loss  of  coherence  in  a scattered  field.  Practical  photodetector  parameters  are 
given  in  section  6. 

For  sufficiently  high  gain,  the  thermal  noise  term  can  be  neglected.  When  the 
background-generated  current,  SjPj,,  is  large  compared  to  the  dark  current,  the  expression 
for  background-limited  detection  is  obtained  as 


S/N  (background  limited)  = — 

2qB  (Pj,  + Pc) 

For  systems  operating  in  large  backgrounds,  Pj^  may  be  considerably  larger  than  Pc.  In  this 
case,  a decrease  in  optical  filter  bandwidth.  B0  (note  that  P^  °c  B0),  gives  a proportional 
improvement  in  S/N:  reduction  of  bandwidth  by  a factor  of  2 improves  S/N  by  3 dB 
(A  reduction  in  the  electrical  bandwidth,  B.  has  the  same  effect,  but  B must  be  wide  enough 
to  respond  to  the  pulse  envelope.)  An  increase  of  2 in  Sj  (or  quantum  efficiency)  improves 
S/N  by  3 dB.  The  S/N  is  increased  6 dB  by  a doubling  of  transmitter  power  (or  Pc),  rather 
than  3 dB  as  might  be  expected  intuitively.  A doubling  of  the  receiving  antenna  aperture 
doubles  both  Pc  and  Pj,  and  thus  increases  Pf/Pb  by  3 dB. 

The  S/N  expression  indicates  that  it  is  advantageous  to  concentrate  the  laser  pulse 
energy  in  a short  pulse  in  order  to  obtain  high  peak  power.  Cutting  the  duration  in  half 
gives  a 3-dB  S/N  advantage,  since  P^is  increased  6 dB  while  the  noise  is  increased  3 dB.  due 
to  a doubling  ot  the  required  electrical  bandwith.  B.  Pulse  lengthening  due  to  propagation 
time  dispersion  will  limit  this  improvement. 


MODULATION 

Pulse  modulation  has  been  assumed  for  both  the  scatter  and  active  RPV  links.  Common 
types  of  pulse  systems  are  pulse  amplitude,  pulse  position  (PPM),  and  pulse  internal  (P1M ) 
modulation.  Because  of  possible  amplitude  fluctuations  in  the  scatter  link,  and  primarily 
because  of  the  low  pulse  rate  (for  high  peak  power),  pulse-amplitude  modulation  is  not  con- 
sidered for  the  active  links.  Pulse-position  modulation  requires  accurate  synchronization 
(for  example,  by  the  use  of  precision  frequency  standards)  but  has  some  advantage  over  PIM 
in  S/N  and  covertness,  in  both  communications  and  IFF. 

Pulse-time  dispersion  will  limit  the  width  of  independent  or  orthogonal  intervals  for 
both  PPM  and  PIM.  For  a scatter  link,  the  time  dispersion  can  be  as  large  as  2 microseconds 
(table  5)  for  a fan  beam  and  a field  of  view  of  4 degrees.  However,  for  likely  system  para- 
meters the  dispersion  should  be  considerably  less  than  0.5  microsecond,  and  this  value  will 
generally  be  considered  for  the  independent  interval. 
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IFF 


IFF  coding  tradeoffs  between  false-alarm  rates  and  probability  of  detection  will  vary 
with  Fleet  vehicles  and  tactics.  For  example,  a submarine  wants  to  have  a very  low  proba- 
bility of  being  triggered  into  an  unnecessary  transmission.  False-alarm  probability  changes 
very  rapidly  with  S/N  so  that  a reasonable  example  establishes  the  approximate  S/N  require- 
ment. 

An  example  of  an  IFF  code  consisting  of  three  pulses  has  been  analyzed7 . The 
required  S/N,  for  a false  alarm  rate  of  0.5  per  hour  with  a detection  probability  of  0.99,  was 
calculated  for  a range  of  pulse  rates  from  10  to  200  pps,  and  from  1 to  4096  time-code 
intervals  per  pulse.  A S/N  of  12  to  13.5  dB  covered  the  range.  A 15-dB  S/N  appears  to  be 
a conservative  requirement. 

Curves  of  detection  probability  as  a function  of  false  alarm  and  S/N  are  included  to 
aid  in  the  design  of  specific  IFF  codes.  Figure  2 is  based  on  a Rayleigh  noise  distribution, 
or  the  envelope  of  Gaussian  noise,  and  would  apply  for  daytime  background  conditions. 

Pn  in  the  figure  is  the  probability  that  the  noise  exceeds  a threshold,  while  Pj  is  the  prob- 
ability that  the  signal  (plus  noise)  exceeds  the  same  threshold.  This  curve  is  for  single 
pulse  probabilities.  (See  Anderson7  for  computations  involving  more  than  one  pulse.)  If 
an  IFF  code  is  to  be  created  specifically  for  night  use.  Poisson  statistics  would  probably 
apply.  Probability  curves  for  this  case  are  given  in  figures  2 and  3. 


c:«ute  2.  Probability  of  detection  versus  S/N  ratio. 


Naval  I iectronics  Laboratory  Center  Technical  Document  447,  Peris<<>t>e  Aircraft  Submarine  Optical 
Communication s System  (Pi  RASOCSj:  Preliminary  Definition  and  Analysis,  by  R Anderson,  D Adrian, 
and  WR  Stone,  1 September  1 ‘>74 
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COMMUNICATION  RATE 


When  the  system  switches  from  acquisition  to  communications,  the  tradeoff  between 
false  alarm  and  detection  probability  changes.  That  is,  additional  false  indications  are 
acceptable  if  the  pulse  detection  can  be  enhanced  enough  to  reduce  the  overall  error  rate. 

In  a fixed-threshold  system,  the  threshold  level  would  have  to  be  altered.  However,  another 
strategy,  such  as  maximum-likelihood  detection,  should  probably  be  used. 

The  information  rate  for  a pulse  position  modulation  system  has  been  shown  to  be 
approximately:9’10 


H = (log2M)/MAT  - PE  j (log2M)MAT| , 

where 

M = number  of  slots, 

AT  = duration  of  slot,  and 
PE  = probability  of  error. 

The  PE  term  is  small  for  the  S/N  considered  here  for  communications  and  is  ignored  in  ob- 
taining approximate  information  rate.  The  information  rate  for  PIM  should  be  close  to 
that  of  PPM  with  some  capacity  used  for  synchronization.  To  obtain  the  maximum  infor- 
mation rate  for  a given  pulse  rate,  assume  that  the  period  between  pulses  is  divided  into 
slots  of  0.5-microsecond  duration.  This  duration  is  dictated  by  pulse-time  dispersion  in  a 
scatter  link.  For  a LOS  link,  the  slot  duration  could  probably  be  shortfv.  Table  5 gives  the 
information  capacity  for  several  pulse  rates.  It  is  noted  that  the  capacity  increases  with  a 
decrease  in  M (increase  in  pulse  rate).  The  information  rate  is  also  given  assuming  one  1 2- 
bif  word  p>  - nulse.  Optical  communication  systems  using  12-bit  words  have  been  imple- 
mented1 1 and  it  is  noticed  that  very  little  is  generally  gained  by  going  to  more  than  4096 
slots.  ^ 


TABLE  5.  INFORMATION  RATfc. 


Pulse  Rate 

Period 

(msec) 

M 

(no) 

H 

(bits/sec) 

Data  Rate  Using 
12  bits/pulse 
(bits/sec) 

10 

too 

200.000 

176 

120 

20 

50 

i 00 .000 

332 

240 

50 

20 

40.000 

265 

600 

100 

10 

20.000 

1430 

1200 

200 

5 

10.000 

2660 

2400 

" RCA,  Electro-Optics  Handbook . i 074 

' Katp,  S ami  RM  Gagliardi.  "The  design  ot  a Pulse-Position  Modulated  Optical  Communication  System," 
I LET.  Transactions  on  Communication  Technology , COM-17,  p 670  to  676,  1 060 
‘“Karp.  S,  EL  O’Neill  and  RM  Gagliardi,  “Communication  Theory  for  the  Free-Space  Optical  Channel, 
Proceedings  of  the  IEEE,  v 10,  p 1611  to  1625,  1970 
1 1 Naval  Electronics  Laboratory  Center  Technical  Note  2672,  Voice  Communication  with  a Elashktmp- 
/:  riled  Ulue-dreen  Dye  Laser,  by  JH  Celto.  EJ  Schimitschck  and  JA  Trias,  24  April  1974* 

•NfcLC  technical  notes  aie  informal  documents  intended  primarily  for  use  within  the  Center. 
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SECTION  5.  LINK  ANALYSIS 


The  link  analysis  is  separated  into  four  parts.  First,  the  transmission  loss  as  a function 
of  range  and  wavelength  is  given  for  both  the  scatter  and  direct  (relay)  propagation  links. 

The  optical  background  noise  is  also  described.  The  S/N  is  detennined  as  a function  of  range 
and  wavelength  for  nominal  system  parameters.  The  value  of  this  computation  is  that  it  indi- 
cates the  preferred  operating  wavelengths  and  allows  a rapid  assessment  of  the  S/N  resulting 
from  a change  in  system  parameters.  Heavy  rain,  dense  fog,  and  clouds  can  produce  losses 
ot  from  20  to  30  dB  per  kilometer  and  obliterate  an  optical  channel.  Estimates  of  link 
availability  are  also  given. 


TRANSMISSION  LOSS 

In  the  visible  portion  of  the  spectrum,  the  major  loss  mechanism  is  aerosol  scattering. 
The  aerosol  scattering  also  permits  OTH  communications  as  discussed  quantitatively  in 
appendix  E.  Near  the  low  wavelength  end  of  the  visible  region,  molecular  (Rayleigh)  scat- 
tering begins  to  contribute  to  the  loss.  In  the  ultraviolet  (UV)  region,  Rayleigh  scattering 
becomes  more  important,  and  absorption  by  ozone  causes  very  large  losses  for  wavelengths 
below  0.3  pm.  Because  of  these  losses,  the  UV  region  has  not  been  considered  in  detail. 

In  the  1R  region,  the  effects  of  absorption  by  water  vapor  and  carbon  dioxide  need 
to  be  considered.  These  losses  are  minor  for  windows  in  the  near  IR  but  become  pronounced 
at  larger  wavelengths.  The  absorption  losses  are  superimposed  upon  the  scatter  losses.  IR 
absorption  loss  is  first  considered,  and  then  OTH  propagation  in  the  visible  and  IR  regions 
is  discussed. 


IR  ABSORPTION  LOSS 

In  :he  absence  of  haze  or  clouds,  the  major  loss  mechanism  in  the  IR  region  is  ab- 
sorption by  atmospheric  components,  particularly  water  vapor  and  carbon  dioxide.  The 
computation  of  radiative  transfer  is  very  complicated  and  numerous  approximate  methods 
have  been  applied.1 2 For  a single  frequency,  the  transmittance  over  a homogeneous  path 
is  of  the  same  general  form  as  for  scatter,  ie, 

T(X)  = exp  (-k(X)w), 


where 


T is  the  transmittance, 
k is  a coefficient, 

X is  wavelength,  and 
w is  the  amount  of  absorber  in  the  path. 


' 2 LaRocca,  AJ,  “Methods  of  Calculating  Atmospheric  Transmittance  and  Radiance  in  the  Infrared,” 
Proceedings  of  the  IEEE,  v 6.3,  p 75  to  94,  1975 
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However,  considering  absorption  over  an  entire  absorption  “line”  and  the  effects  of  u xcent 
absorption  lines  leads  to  various  expressions  for  transmittance.  For  weak  absorption  lines  or 
short  distances,  the  above  form  is  often  applicable.  For  strong  absorption  lines  or  long  dis- 
tances the  following  form  is  often  used:1 4 

T(X)  = exp  (-k(X)w)l/,J . 

This  expression  is  often  used  for  water  vapor  in  part  of  the  1R  spectrum,  and  the  following 
expression  is  used  for  other  1R  wavelengths.1 4 

T(X)  = exp  [-k(X)w/(  1 + 2k2(X)w)l/:il  . 

For  large  w,  this  approximates  the  previous  form  for  transmittance. 

Transmittance  tables  for  H2O  and  CO2  are  given  by  Valley  (tables  10-13  and  1 0 1 4 >. 1 3 
The  table  for  CO2  covers  the  range  of  absorption  concentration  of  interest  in  the  present 
study.  The  CO2  concentration  is  about  0.294  pascal  (30  atm-cm)  per  kilometer  of  sea  level 
path  and  decreases  with  altitude.  A value  of  0.245  pascal  (25  atm-cm),  corresponding  to  an 
average  altitude  of  1 km,  was  used  to  obtain  the  absorption  for  transmission  windows  be- 
tween 1 and  4.6  pm.  For  CO-i  laser  transmission  at  10.6  pm.  the  laser  transition  occurs 
in  reverse  in  the  atmospheric  CO->.  This  absorption  is  combined  with  water  vapor  absorp- 
tion. A typical  loss  value  is  about  1 dB/1.6km. 

Absorption  loss  in  water  depends  upon  the  total  water  vapor  concentration  over  the 
path  (scatter  depends  on  drop  size).  The  relationship  between  precipitable  water  and  rela- 
tive humidity  and  temperature  is  shown  in  figure  4.  In  determining  absorption,  we  assume  a 
temperature  of  20°C  and  50-percent  relative  humidity.  The  0.9  cm/km  is  reduced  to  0.84 
to  account  for  a drop  in  water  content  with  altitude.  The  transmittance  table  in  Valley  is 
only  applicable  to  50  cm  of  water  (60-km  path).  The  transmittance  values  for  longer 
ranges  were  generally  extrapolated  assuming  an  exponential  (kw>1/2  variation.  When 
the  variation  in  the  table  was  still  close  to  exponential  (w)  at  maximum  w.  this  variation  was 
used  to  2 w and  then  the  w 1/2  variation  was  assumed.  With  the  extrapolation  used,  loss 
values  for  ranges  beyond  about  100  km  are  somewhat  suspect. 

The  absorption  loss  as  a function  of  range  for  various  transmission  windows  is  plotted 
in  figure  5.  The  loss  is  similar  for  wavelengths  close  to  those  indicated  but  generally  becomes 
extremely  large  at  intermediate  wavelengths.  The  1R  absorption  loss  is  combined  with  the 
scatter  loss  in  the  next  section  for  the  different  types  of  propagation  channels.  Absorption 
loss  for  wavelengths  below  I pm  is  considered  negligible  and  is  ignored. 


,3Unitcd  States  Air  Force,  Handbook  of  Geophysics  and  Space  Environment,  SL  Valley,  Editor,  1965 
1 4 Naval  Electronics  Laboratory  Center  Technical  Note  2714.  Radiated  Background  Noise  in  the  10-45 
GHz  Band,  by  DB  Sailors,  DJ  Adrian  and  Pll  Levine  24  June  1974* 

*NLLC  technical  notes  are  informal  documents  intended  primarily  for  use  within  the  Center. 
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figure  4.  Water-vapor  concentration  as  a function  of  temperature  and 
relative  humidity  (ref  13). 


f igure  5.  Infrared  absorption  loss. 


FORWARD-SCATTER  LINK 


A forward-scatter  optical  communication  channel,  where  the  scattering  arises  from 
atmospheric  “haze"  aerosols,  has  been  analyzed.2  The  referenced  report  is  included  as 
appendix  E for  convenience.  As  illustrated  in  figure  6,  the  two  primary  mechanisms  for 
OTH  propagation  (with  absence  of  clouds)  are  single  and  multiple  scattering  from  atmospheric 
aerosols.  Ignoring  absorption  and  noting  that  scattering  tends  to  be  peaked  in  the  forward 
direction,  it  is  meaningful  to  think  in  terms  of  an  “aura”  of  multiply  (forward)  scattered 
photons  surrounding  the  laser  beam.  This  aura  will,  in  general,  attenuate  less  rapidly  with 
distance  than  the  direct  beam  and,  were  it  not  for  the  effect  of  the  nearby  Earth’s  surface, 
would  be  the  dominant  source  of  radiance  at  long  ranges.  However,  the  effect  of  the  Earth's 
surface  is  to  remove  energy  from  the  aura  by  a combination  of  absorption  and  diffuse  scat- 
tering so  that,  for  near-surface  transmitter/receiver  links,  the  aura  will  attenuate  more  rapidly 
with  range  than  that  of  the  singly  scattered  component.  For  this  reason,  the  computed  trans- 
mission loss  is  for  the  singly  scattered  component.  It  should  be  borne  in  mind,  however,  that 
at  ranges  less  than  about  100  kilometers,  the  aura  can,  under  favorable  conditions,2  be  the 
dominant  component  (by  about  10  dB). 

A BASIC  language  program  (see  appendix  B)  was  written  to  compute  the  scatter 
transmission  loss  as  a function  of  visibility,  wavelength,  distance,  and  transmitter  vertical  beam 
width.  The  scatter  model  was  developed  in  terms  of  the  extinction  coefficient.  As  discussed 
in  the  appendix,  the  commonly  used  simple  relationship  between  extinction  coefficients  and 
visibility  or  meteorological  range,  underestimates  the  extinction  coefficients  by  as  much  as  an 
order  of  magnitude  for  some  wavelengths  and  visibilities.  However,  it  has  been  possible  to 
combine  the  semi-empirical  model  of  Elterman  with  other  theoretical  and  experimental 
results  to  obtain  a model  of  sufficient  accuracy  for  the  range  of  visibilities  and  wavelengths 
required  for  the  system  tradeoff  studies.  This  is  discussed  further  in  appendix  A. 


Figure  6 Beyond-thehori/on  propagaiion  mechanisms  (cloud-free  case). 
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Ihe  computed  transmission  loss  does  not  include  the  aerosol  scattering  “gain"  in  the 
forward  direction  but  it  is  included  in  the  S/N  computations  presented  later.  The  absorption 
loss  in  the  IK  band  has  been  combined  with  the  scatter  transmission  as  a multiplicative  factor 
in  the  data  presented  here.  Dotted  lines  connect  the  IR  windows  for  ease  of  viewing  the  data 
Absorption  between  windows  is  generally  such  as  to  preclude  communications  at  any  signi- 
ficant distance. 

Figure  7 shows  the  path  loss  versus  wavelength  and  distance  for  a pencil  beam.  The 
vertical  lapse  rate  of  the  aerosol  concentration  is  1.3  km  for  the  figure  This  figure 
indicates  the  advantage  of  a system  operating  in  the  near  IR.  The  background  noise  varies 
with  wavelength  as  shown  in  the  next  section  and  the  resulting  S/N  indicates  an  even  greater 
advantage  lor  near  IR  (eg.  1 .06  ^m)  compared  to  visible  (eg.  0.53  pm).  Experimental  com 
panson  of  scatter  at  1 .06  and  0.53  jam  would  appear  to  be  an  important  step  in  evaluating 
the  forward-scatter  model 

Figure  8 shows  the  loss  assuming  a 4 degree  vertical  fan-beam  transmitter.  The  other 
parameters  are  the  same  as  for  figure  7.  Comparison  of  the  two  figures  shows  the  signal  enhance- 
ment produced  at  long  ranges  by  the  fan  beam.  As  discussed,  (in  LINK  S/N  RATIOS)  this 
advantage  will  be  effected  by  time  dispersion  for  short-pulse  systems.  It  is  also  noted  that 
the  fan  beams  tend  to  shift  the  wavelength  for  minimum  path  loss  towards  the  lower  IR 
range.  Considering  available  laser  sources.  1 .06  „m  is  clearly  the  most  probable  operating 
Tht  path  ,OSS  at  l 06/im  as  a function  of  visibility  and  range  is  shown  in 
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Figure  9.  Path  loss  at  1 .06  /im,  scatter  link. 


LINE-OF-SIGHT 

A line-of-sight  (LOS)  propagation  path  is  involved  in  the  active  relay  link.  A program 
to  compute  the  transmittance  as  a function  of  visibility,  wavelength  and  distance  is  described 
in  appendix  C.  The  program  also  computes  the  S/N  when  the  applicable  system  parameters 
are  specified.  Transmittance  as  a function  of  wavelength  and  distance  is  shown  in  figure  10 
for  a 20-km  visibility.  The  distance  indicated  is  the  surface  distance  while  the  transmit- 
tance is  computed  for  a path  between  a surface  terminal  and  an  airborne  terminal  at  the 
minimum  altitude  necessary  to  maintain  a line-of-sight  path.  The  IR  absorption  loss  is  super- 
imposed upon  the  scatter  loss.  Geometric  loss  due  to  beam  divergence  is  not  included  in  the 
path  loss  but  is,  of  course,  included  later  in  the  S/N  calculation.  It  is  seen  that  the  near  IR  is 
preferable  to  the  visible  for  LOS  links  as  well  as  for  the  scatter  link  described  in  the  previous 
section. 

Figure  1 1 shows  the  path  loss  to  an  RPV  at  2.5  km  altitude  as  a function  of 
range  and  visibility.  It  is  noted  that  the  higher  RPV  altitude  gives  considerable  smaller  path 
loss  than  shown  in  figure  10.  As  range  is  increased,  the  loss  incre;  ses  more  rapidly  since  a 
lower  elevation  angle  is  required  (constant  RPV  altitude),  thus  increasing  the  fraction  of  the 
path  in  the  tow-altitude,  high-loss  region. 
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Figure  1 1 . Path  loss  as  a function  of  range  at  1 .06  /urn. 


BACKGROUND  NOISE 

Noise  and  noise  sources  may  be  characterized  by  source  distribution,  spectral  distri- 
bution, and  statistical  characteristics.  Statistical  characteristics  of  noise  in  optical  systems 
were  discussed  in  the  section  on  signal  modems.  Sources  and  spectral  distributions  are  dis- 
cussed here. 


SOURCE  DISTRIBUTIONS 

Some  contrasting  characteristics  of  source  location  are:  internal  or  external  (back- 
ground); uniformly  distributed  or  localized;  and  stationary  or  nonstationary. 

Internal  noise  may  be  of  two  types.  Thermal  noise  originates  within  the  amplifiers 
and  load  resistors  of  the  receiver.  Internally  generated  quantum  noise  is  caused  by  dark 
current  in  a photoemissive  or  photovoltaic  detector.  The  quantrm  or  photon-fluctuation 
noise  produced  by  external  radiation  is  proportional  to  the  powt  incident  on  the  detector. 
Thus,  both  the  signal  and  background  radiation  contribute  to  the  noise.  The  photon  noise 
associated  with  the  signal  (self  noise)  is  included  in  the  program  listed  in  the  appendix 
although,  at  long  operating  ranges,  it  is  negligible  compared  to  the  background. 
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The  importance  of  the  distributed  versus  discrete  source  has  to  do  with  the  measure 
of  background.  If  the  source  fills  the  receiver  field  of  view,  the  spectral  radiance,  N(X), 

(watt  m-“  steradian-'  /im-' ) is  an  appropriate  measure.  However,  if  the  source  does  not 
fill  the  field  of  view  (discrete  source)  the  irradiance,  H(A),  (watt  m~-  /im-* ) is  more 
appropriate. 

When  a source  is  contained  within  the  receiver  field  of  view,  the  received  background 
power,  Pb,  in  terms  of  H is  given  by 

Pb  = H(\)Art0taB0  , 

where 

Ar  = receiving  aperture, 

B0  = optical  bandwidth, 
tQ  = transmittance  of  the  receive  optics, 
ta  = transmittance  of  the  atmosphere,  and 
X = wavelength. 

When  a source  completely  fills  the  receiver  field  of  view,  the  received  power  may  be  expressed 
in  terms  of  the  radiance,  N,  as 

Pb  - N(X)ftrArt0taB0  , 


where 


£2r  = solid  angle  receiver  field  of  view. 

For  a conical  field  of  view,  the  solid  angle  is  related  to  the  plane  angle  of  view,  0T,  by 
S2r  = 2tt(  1-cos  0tI 2) 


wer 


, for  small  angles. 


This  latter  relation  holds  for  systems  considered  in  this  study.  For  a nonconical  field  of 
view.  eg.  cylindrical  optics,  the  effective  solid  angle  needs  to  be  either  determined,  or  roughly 
approximated  by  replacing  Of  in  the  above  expression  by  the  product  of  the  linear  FOVs 
in  the  largest  and  the  smallest  direction. 

The  major  background  radiation  sources  are  the  sun  and  sky.  The  sun  may  be  con- 
sidered a discrete  or  a distributed  source  depending  upon  whether  the  receiver  FOV  is  larger  or 
smaller  than  the  0.5-degree  angle  subtended  by  the  sun.  The  spectral  distribution  for  the  sun 
is  given  in  the  next  section,  but  it  is  generally  assumed  in  this  study  that  the  receiver  does 
not  look  directly  at  the  sun.  The  sky  background  may  be  considered  as  a distributed  uniform 
source  for  any  FOV  considered  in  this  study.  The  relative  locations  of  natural  sources  may 
be  considered  approximately  stationary  when  viewed  from  a ship.  This  may  no  longer  be 
true  for  an  aircraft. 
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Intentional  or  unintentional  man-made  interference  sources  will  generally  be  point 
sources.  However,  scattering  of  this  energy,  eg  by  clouds,  may  generate  a distributed  inter- 
ference source.  Interference  source  locutions  may  well  be  nonstationary. 

SPECTRAL  CHARACTERISTICS 

The  spectral  distribution  for  the  sun  is  shown  in  figure  1 2.' s It  is  noted  that  a black- 
body  at  5999  K fits  the  curve  quite  well,  in  fact,  the  blackbody  approximation  is  reasonably 
good  through  the  millimetre  portion  of  the  radio  band.14  Figure  12  also  shows  the  spectral 
distribution  at  sea  level.  The  effect  of  absorption  by  water  vapor,  carbon  dioxide,  and  ozone 
is  apparent.  The  spectral  irradiamcc  within  the  atmosphere  depends  upon  the  density  and 
composition  of  the  air  along  the  path  of  the  beam.  Figure  13  illustrates  the  depletion  of  the 
solar  beam  in  passing  through  clear  air  to  sea  level  as  a function  of  the  optical  air  mass  and 
zenith  angle.  Detailed  tables  of  solar  irradiance  are  given  in  Valley.1* 

The  daytime  sky  background  radiance  is  due  to  scattering  of  incident  solar  radiation 
and  to  emission  by  atmospheric  particles  heated  by  incident  radiation.  The  visible  radiance 
is  predominately  scattered  radiation,  and  the  infrared  is  mainly  atmospheric  emission  Both 
the  level  and  the  spectral  distribution  depend  upon  the  characteristics  of  the  scatterer. 
Radiance  from  sunlit  clouds  may  be  an  order  of  magnitude  higher  than  from  a clear  sky. 

Sky  background  data  for  several  conditions  are  shown  in  figure  14.  Atmospheric  absorption 
bands,  illustrated  in  figures  I 2 and  13,  are  not  included.  The  dashed  line  is  discussed  in  the 
next  section  of  this  report. 

Night  sky  background  radiance  data  are  shown  in  figure  15.  The  main  source  in 
the  visible  is  scattered  moonlight.  At  night,  aurora  and  airglow  can  be  important  for  some 
wavelengths  near  \ t>  /am.  Vbermul  emission  predominates  for  wavelengths  greater  than  2.5 
/am.  in  contrast  to  a value  of  4 to  5 /am  during  the  day. 


ANALYTIC  AL  APPROXIMATION  TO  BACKGROUND 

An  analytical  approximation  to  the  background  was  needed  for  the  link  performance 
computer  program.  From  figures  14  and  15,  it  is  seen  that  a 283  K blackbody  radiation  curve 
serves  as  a fairly  good  envelope  to  the  measured  values.  Examination  of  curves  in  the  IR 
band'*’*  indicate  that  288  K ma>  be  slightly  better.  In  fact,  the  283  K indicated  for  figures 
14  and  15  may  well  be  a misprint.  The  ruihance  may  be  derived  from  Planck’s  law  for  spectral 
radiant  emittance,  W(X), 


W(X)  = 


2trc-h  I 

\5  expjhc/XkT-  1 


1 Kopcika,  NS  and  J Boldnpna.  “Bin.  kgvoond  Noise  in  Optical  Communication  Systems,”  Proceedings  nt 
ih,  IMF.  v 10.  p 1571  to  1577.  |d70 
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Figure  12.  Spectral  distribution  curves  related  to  the  sun;  shaded  areas 
indicate  absorption,  at  sea  level,  due  to  the  atmospheric 
constituents  shown  (ref  13)  . 
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Figure  13.  Solar  spectral  irradiance  curves  at  sea  level  for  various  optical 
air  masses;  the  value  of  thc^olar  constant  used  in  this 
calculation  was  1322  W m'*(ref  13). 
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Figure  1 4.  Day  sky  background  data  ( ref  15). 


Figure  15.  Night  sky  background  data  (ref  15)  . 
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where 


c = velocity  of  light, 

X = wavelength, 

T = temperature  of  background  radiation  source, 
h = Planck’s  constant,  and 
k = Boltzmann’s  constant 

If  it  is  assumed  that  the  source  is  a diffuse  radiator  (Lambert  surface),  it  can  be  shown  that 1 p 
N(X)=  W (X)/tt  . 

The  expression  for  N with  T = 288  K is  used  for  the  emission  portion  of  the  sky  noise. 

The  scattering  of  sunlight  is  a function  of  wavelength,  so  that  a blackbody  term  alone 
cannot  be  used  for  this  component.  Comparison  of  the  blackbody  radiation  for  5900  K and 
the  measured  data  yielded  an  empirical  multiplying  factor  of 

4 X 10'5/X3/2, 

where  X is  in  /a in.  The  empirical  expression  for  the  daytime  background  is  then 


4.75  X |p4  1.2  X IQ8 

xn/2(e2.44A_,)  * x5(t50M_!) 


(daytime). 


where  X is  in  /am.  The  dashed  curve  in  figure  14  shows  a plot  of  this  expression.  A similar 
approach  for  the  night  sky  background  gives 


N(X)  = 


8.3  X 1 0-3 

X)V2  (,244/X.i) 


1.2  X 10» 
X5 


(night). 


This  expression  is  plotted  as  the  dashed  curve  in  figure  1 5. 


LINK  S/N  RATIOS 


PULSE  DISPERSION  EFFECTS 

The  transmission  loss  given  earlier  is  actually  the  value  for  loss  of  energy.  The  curves 
also  give  power  loss  for  continuous-wave  transmission  or  pulse  transmission  where  the  pulse 
duration  is  long  compared  to  the  time  dispersion  introduced  by  the  media.  However,  time- 
dispersion  effects  must  be  considered  when  obtaining  the  S/N  for  a short-pulse  system.  The 
received  energy  will  be  spread  over  a time  equal  to  the  sum  of  the  pulse  duration  and  dis- 
persion. 

We  first  consider  the  time  dispersion  for  a pencil  transmitter  beam  and  a finite  field 
of  view.  Referring  to  figure  E3  of  appendix  E,  let  the  elevation  angles  be  equal  (ie,  (3|=0  2=0) 
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and  denote  the  dashed  line  by  Z3.  The  difference  hi  propagation  time  within  the  receiver 
field  of  view,  F,  (F  = dpi  in  the  figure)  is: 

AT  fiZ,+Z3-Z2]  . 

sin  F + sin  Ip 
sin(20+F) 


where  c is  the  velocity  of  light.  The  law  of  sines  was  used  in  the  derivation.  Next  we  let: 

I)  = 2Z-i  and 
P = D/2R. 

where  D is  the  distance  between  receiver  and  transmitter  and  R is  the  “4/3  earth”  radius 

(see  appendix  E)  of  8493  kilometres.  After  some  trigonometric  (manipulations  and  replacing 

their  arguments,  we  obtain  the  time  dispersion  for  a pencil  beam: 

/ 

/ 
r 

(pencil  beam). 

dispersion  in  microseconds  for  several  distances. 

TABLE  6.  TIME  DISPERSION. 


the  sines  of  angles  by 


AT  = 


D“F 


4 Re- 

Table  6 lists  the  time 


Time  in  Microseconds 

Path  D 

Pencil  Beam 

Fan,  Beam 

in  km 

FOV 

Fan  Angle  = FOV 

1° 

4° 

1° 

4° 

50 

0.004 

0.017 

0.025 

0.41 

100 

0.017 

0.069 

0.050 

0.81 

200 

0.060 

0.28 

0.10 

1.6 

300 

0.156 

0.62 

0.15 

2.4 

For  a fan  beam,  we  assume  an  equal  fan  and  FOV  angle,  F.  The  difference  in  propagation 
time  between  the  upper  and  lower  edge  of  the  beam  is  then  found  to  be  approximately 
(assuming  lower  edge  as  a fiat  earth): 

AT  = — j-*  -r  (1/cosF-l)  (fan  beam). 

Time-dispersion  values  are  shown  in  table  6.  Some  of  the  dispersion  values  are  significantly 
longer  than  typical  pulse  lengths  (eg  about  0.02  /usee  for  a Nd:Y  AG  laser). 

If  the  dispersion  equals  the  pulse  duration  (eg,  at  100  km,  for  a pencil  beam  and  I- 
degree  field  of  view),  the  peak  power  is  reduced  to  half  that  expected  from  the  path-loss 
computations.  Further  increase  in  the  FOV  would  increase  the  signal  energy  but  not  the 
peak  power.  The  increased  energy  comes  from  increased  pulse  stretching.  Consider  the 
S/N  expression. 
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SNR  (background  limited)  = 

The  required  bandwidth,  B,  is  inversely  proportional  to  the  received  pulse  duration,  T+AT, 
where  T is  the  transmitted  pulse  duration  and  AT  the  dispersion.  Since  Pc  is  also  inversely 
proportional  to  the  received  pulse  duration,  the  S/N  is  approximately, 


Sd  V 

2qB  (Pb+Pc) 


SNR  = 


K 

T+AT 


Once  T is  smaller  than  AT,  there  is  little  advantage  in  reducing  it  further.  It  is  noted  that 
the  bandwidth  must  be  varied  as  AT  varies.  This  could  be  done  stepwise,  for  example,  with 
2-to-l  changes  in  bandwidth  (3-dB  changes  in  noise). 

Now  consider  the  effect  of  varying  the  FOV.  Once  the  FOV  is  large  enough  that  the 
time  dispersion  is  comparable  to  the  pulse  duration,  a further  increase  will  give  no  increase  in 
peak  Pc.  The  background  P^  will  increase  but  B may  be  reduced,  thus  keeping  the  S/N  con- 
stant. It  is  advantageous  to  increase  the  FOV  until  the  background  current  term  exceeds 
the  dark  current.  In  addition,  a wider  FOV  would  relax  pointing  accuracy  requirements. 

The  time  dispersion  or  extended  source  effects  might  conceivably  be  put  to  some 
use.  such  as  angle  diversity.  For  example,  2 detectors  could  be  used,  each  with  a different 
elevation  look  angle.  The  pulse  received  at  the  higher  angle  should  lag  the  other  by  an 
amount  depending  upon  the  transmitter  distance.  This  would  allow  discrimination  against 
interference  (intentional  or  unintentional)  with  a different  time  relationship.  Some  advan- 
tages of  diversity  reception  could  also  be  achieved  in  that  2 scattering  volumes  with  possibly 
different  short-term  characteristics  would  be  viewed. 


FORWARD-SCATTER  LINK 

The  program  described  in  appendix  B was  used  to  compute  the  S/N  data  given  in 
this  section.  The  IR  absorption  loss  discussed  earlier  in  this  section  is  included  in  the  data. 
Dashed  lines  connect  the  IR  windows  for  ease  of  following  the  curves.  In  general,  there  will 
be  excessive  absorption  loss  in  the  wavelength  regions  indicated  by  the  dashed  lines.  The 
computed  values  are  applicable  when  pulse  dispersion  can  be  ignored.  For  some  of  the  data, 
corrections  are  applied  for  time-dispersion  effects  as  described  in  the  last  section. 

Figures  16  through  18  show  the  S/N  variation  as  a function  of  wavelength  and  dis- 
tance for  three  different  visibilities  using  a pencil  transmitter  beam  pointing  at  the  horizon. 
The  noise  level  for  most  of  the  curves  is  the  maximum  noise  (analytic  expression).  The 
results  for  the  noise  curve,  D of  figure  14  is  indicated  in  figure  16.  The  system  parameters 
given  in  table  7 are  used  independently  of  wavelength,  in  order  to  indicate  the  generally 
preferred  wavelength  region.  As  can  be  seen  in  section  6,  there  is  not  generally  a strong  regu- 
lar variation  of  system  parameters  with  wavelength.  The  figures  indicate  the  advantage  of 
the  near  IK.  Consideration  of  available  laser  sources  anil  detectors  indicates  that  1 ,0o  /am 
is  the  best  operating  frequency  at  present.  The  3.N-/um  region  should  be  considered  in  the 
tuture  if  component  advances,  especially  chemical  (deuterium  fluoride)  lasers,  continue. 
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TABLE  7.  NOMINAL  SYSTHM  PARAMETERS. 


Parameter 

Value 

Transmitter 

Radiated  Power 

I06  W 

Beam  Divergence 

4°  fan,  or  pencil  beam 

Height 

25  m 

Receiver 

Field  of  View 

1° 

Aperture  Area 

10~7m7 

Optical  Bandwidth 

I0'4  gim  ( 1 A) 

Detector  Responsivity 

0.065  A/W 

Dark  Current 

I0-14  A 

Electrical  Bandwidth 

7 X I07  Hz 

Height 

25  m 

(Background  or  dark  current  limited  operation  is  used,  based  on  the  use  of 
a detector  with  adequate  gain). 


Figure  1 9 compares  S/N  for  different  aerosol  distributions  using  a fan  beam.  The 
best  results  are  given  for  a marine  aerosol,  which  would  be  expected  over  most  of  the  ocean. 
Continental  hazes,  which  may  exist  near  coast  lines,  are  expected  to  generally  produce  signals 
a few  dB  smaller  than  marine  aerosols.  The  difference  between  the  S/N  for  the  two  aerosols 
is  greater  in  the  visible  than  in  the  near  IR.  The  Deirmerdjian  forward-scattering  function  for 
a maritime  M haze  at  wavelengths  of  0.45,  1.19,  and  10.0  /am  is  given  in  figure  20.  Similar 
curves  for  a Deirmerdjian  cloud  are  given  in  figure  21 . 

Figure  22  shows  the  S/N  versus  range  for  a 4-degree  fan  beam  with  visibility  as  a 
parameter,  for  a 1 .06-gtm  system.  The  system  parameters  used  were  the  same  as  given  in 
table  7.  The  presently  available  detector  responsivity  at  1.06  gin  is.  however,  only  0.018 
ampere/watt  (A/W)  rather  than  0.065  A/W.  The  curves  would  also  apply  to  far  ranges 
(where  Pj,>  Pc)  for  a responsivity  of  0.018  and  a laser  power  of  1 .9  megawatts  (which  is 
presently  available).  The  effect  of  time  dispersion  on  S/N  (assuming  a 0.015-psec  pulse)  is 
indicated  in  the  figure  for  the  20-km  visibility  curve.  The  same  reduction  would  apply 
to  all  curves.  It  is  assumed  that  the  electrical  filter  bandwidth  is  adjusted  to  the  inverse  of 
the  dispersion  time  for  each  distance.  The  pencil  beam  results  at  V = 20  km  are  shown  in  the 
figure  for  comparison.  It  is  seen  that  the  time  dispersion  reduces  the  advantage  of  the  fan 
beam  for  short  pulses  as  discussed  in  the  previous  section. 

Figure  23  shows  the  effect  of  launch  angle  on  the  S/N  for  a pencil  beam.  A detector 
responsivity  of  0.018  A/W  and  a dark  current  of  10"^  ampere  was  used  in  the  computation. 
Other  parameters  are  as  shown  in  table  7.  It  is  noted  that  a beam  elevated  at  about  4 degrees 
gives  a large  enhancement  in  S/N  at  long  ranges.  The  optimum  elevation  depends  on  aerosol 
distribution  and  range.  Thus,  elevation  scanning  should  be  part  of  the  link-acquisition  proce- 
dure. 
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Figure  23.  Effect  of  launch  angle  on  S/N. 
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Examination  of  figures  18  and  23  shows  that  a 15-dB  S/N  (for  Sj  = 0.018  A/W.  P=1.9 
MW.  X=l.06  gm)  could  he  expected  for  ranges  less  than  60  km  and  100  km  for  visibilities  of 
20  and  50  km  respectively.  The  use  of  a 40-megawatt  transmitter  would  extend  these  ranges 
to  80  and  1 60  km.  The  effects  of  parameter  changes  on  S/N  are  discussed  in  section  4 and 
coverage  predictions  for  an  optimum  set  of  parameters  arc  given  in  section  3. 


ACTIVE  RPV  LINK 

The  received  power  for  an  active  link  is  given  by, 

ArT(R)P, 

Pc  = 7r(0D  R/2)2  ’ 


where 

Ar  = receiving  aperture  area. 

T(R)  = path  loss  at  range,  R, 

Pt  = transmitter  power, 

0 1)  = transmitter  divergence,  and 

R = range. 

Loss  in  the  receiver  optics  is  ignored  but  could  be  considered  to  be  incorporated  in  an  effec- 
tive aperture  area.  The  received  power  was  used  to  determine  S/N.  It  is  assumed  that  suffi- 
cient gain  is  used  in  the  photomultiplier  so  that  the  thermal  noise  of  following  amplifiers  can 
be  ignored.  Dark  current  is  included  in  the  computation  but  has  negligible  effect  upon  day- 
time S/N. 

Using  the  transmitter  and  receiver  parameters  given  in  table  7 the  S/N  was  computed 
for  a range  of  visibilities.  It  is  assumed  that  the  path  is  between  a ship  (system  at  7-m  altitude) 
and  an  RPV  (at  2.5-km  altitude).  Figure  24  shows  that  a 150-nmi  communication  range 
(S/N  = 15  dB)  with  an  RPV  located  at  the  midpoint  (140-km  ship-to-RPV  path)  could  be 
achieved  with  visibilities  greater  than  10  km.  S/N  curves  for  a 10-degree  EOV  should  allow 
a simpler  receiving  pointing  control  system  to  be  used  (ship-to-RPV  link). 

I he  RPV  pointing  accuracy  for  transmission  presents  a more  stringent  requirement. 

The  curves  in  figure  24  assume  that  the  receiver  is  within  the  l-mrad  FOV.  Figure  25  shows 
the  S/N  for  a 20-km  visibility.  Curves  are  also  shown  for  a transmitter  divergence  of  1.8 
degrees  (or  0 ^ = 0.06°  and  P^  = 10-^  W).  It  is  noted  that  the  desired  range  can  be  achieved 
with  a pointing  accuracy  of  about  ± degree  on  both  the  transmitter  and  receiver  with  a visibility 
of  20  km. 

LINK  AVAILABILITY 

Weather  affects  the  link  availability  in  2 ways.  First,  the  path  loss,  in  any  of  the  link 
types,  is  a strong  function  of  visibility  or  meteorological  range.  Second,  clouds  can  affect  the 
links  in  different  ways.  On  an  RPV  link,  a cloud  on  the  path  will  attenuate  the  signal.  In  a 


'F 


' HP 


46 


S/N  (dBI 


40  60  80  100  120  140  160  180  200 


RANGE  (km) 

Figure  24.  S/N  for  active  link,  RPV  altitude  =2.5  km. 
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Figure  25.  S/N  for  active  link,  RPV  altitude  = 2.5  km,  visibility  = 20  km. 
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scatter  link,  a cloud  can  either  attenuate  the  signal  or  increase  the  scatter,  depending  upon 
location.  Operational  strategies  could  probably  be-  developed  to  maximize  scatter-link 
availability  under  cloudy  conditions. 

Contour  maps  showing  the  percent  of  time  that  the  horizontal  visibility  is  less  than  5 
nmi  have  been  published.16  Examples  of  these  maps,  giving  the  percent  of  time  that  the 
visibility  is  less  than  5 nmi,  are  shown  in  figures  26  through  29.  We  are  primarily 
interested  in  visibilities  significantly  greater  than  10  km  (almost  5 nmi).  Examination  of 
the  data  samples  taken  from  nine  different  locations  show  that  a visibility  greater  than  20  km 
occurs  between  10  and  20  percent  (usually  about  15  percent)  less  often  than  a visibility 
greater  than  10  km.17  Thus,  in  the  absence  of  a more  definitive  information,  it  can  be 
estimated  that  the  availability  of  a 20-knt  visibility  would  he  I 5 percent  less  than  shown 
on  the  maps. 

Contour  maps  showing  the  percent  frequency  of  low  clouds  are  shown  in  figures  30 
through  33. 16  Low  clouds  are  defined  as  clouds  with  a base  at  2438  m or  lower  To 
determine  the  availability  of  an  RPV  link,  a comprehensive  analysis,  taking  into  account  the 
base  height  distribution  and  distribution  of  cloud  size  and  spacing,  would  be  required.  Some 
information  for  such  a study  can  be  found  in  Cato18  and  Green  and  Greenwell  19  published 
on  1 May  1975.  The  contour  maps  do,  however,  provide  information  on  the  lower  limit  of 
availability.  That  is,  an  RPV  system  could  operate  some,  as  yet  undetermined,  fraction  of  the 
time  even  with  low  cloud  cover.  As  indicated  earlier,  the  effect  of  clouds  on  a scatter  system 
is  complicated  and  can  sometimes  improve  as  well  as  degrade  the  link.  If  the  cloud  base  is 
greater  than  9|  5 or  1 220  metres,  the  scatter  signal  can  be  greatly  enhanced. 


16  McDonnell  Douglas  Report  1 06 3,  Marine  Weather  of  the  World,  June  1468 

1 7 Naval  Weapons  Center,  China  Lake,  Technical  Note  4056- 16.  Weather  Effects  on  Infrared  Systems  for 
Point  Defence,  by  FE  Nicodemus,  May  1972 

'"Electro-Optical  Systems  Report  4440,  Final  III.  Laser  Systems  Study,  ‘‘art  III,  Effects  of  Clouds,  by  GA 
Cato,  L.W  Carrusand  KJ  von  Essen,  14  December  1975 
19  Naval  Electronics  Laboratory  Center  Technical  Note  2923,  Optical  SATCOM:  Effects  of  Clouds  on 
Global  Availability  of  Optical  Satellite  Communications  Links,  by  SW  Green  and  RA  Greenwell,  I May  1975* 
*NELC  technical  notes  arc  informal  documents  intended  primarily  for  use  within  the  Center. 
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Figure  26.  Percent  time  visibility  5 nautical  miles.  North  Atlantic.  July 
(courtesy  McDonnell  Douglas)  (ref  16). 


Figure  27.  Percent  time  visibility  5 nautical  miles,  North  Atlantic.  January 
(courtesy  McDonnell  Douglas)  (ref  16). 
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Figure  28.  Percent  time  visibility  5 nautical  miles.  North  Pacific.  January 
(courtesy  McDonnell  Douglas)  (ref  16). 


Figure  2'h  Percent  time  visibility  5 nautical  miles,  North  Pacific,  June 
(courtesy  McDonnell  Douglas)  (ref  16). 
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Figure  32.  Percent  frequency  low  cloud  amount  6/10  or  more.  North  Pacific,  January 
(courtesy  McDonnell  Douglas)  (ref  16), 


Figure  33.  Percent  frequency  low  cloud  amount  6/10  or  more.  North  Pacific.  July 
(courtesv  McDonnell  Douglas)  (ref  16). 
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SECTION  6.  COMPONENTS  AND  SUBSYSTEMS 


I lie  major  components  involved  in  system  tradeoff  considerations  are  briefly 
described  in  this  section.  The  range  ot  practical  parameters  for  each  component  is  generally 
listed  as  a function  of  wavelength.  Approximate  costs  are  given  where  available  to  aid  in 
tradeoff  considerations.  Components  to  be  discussed  include  the  following:  lasers,  optical 
filters,  detectors,  pointing  control  systems,  and  remotely  piloted  vehicles  (RPVs). 


LASERS 

Laser  action  has  been  demonstrated  in  various  media  including  gases,  liquids,  crys- 
talline solids,  and  glasses.  A huge  number  of  continuous  and  pulsed  sources  have  been 
developed  with  useful  power  outputs  at  wavelengths  which  span  the  optical  spectrum.  How- 
ever. there  are  a relatively  small  number  of  high  peak-power  laser  types.  These  are  discussed 
briefly  in  this  section.  Summaries  of  laser  types  and  characteristics  are  found  in  Geusic20  and 
Rt  A , while  the  performance  of  commercially  available  lasers  is  given  in  Laser  Focus.2 1 

I he  principal  high-power  pulse  lasers  for  different  wavelength  regions  arc  listed  in 
table  X I ulse  parameters,  beam  divergence,  and  cost  are  listed  for  commercial  units.  It 
should  be  noted  that  peak  power  and  pulse-repetition  rate  can  generally  be  altered,  keeping 
the  average  power  constant.  Dye  and  chemical  lasers  are  under  active  development  and 
much  higher  powers  are  likely  in  the  future. 

Principal  high-power  continuous  wave  laser  types  are  listed  in  table  9.  Costs  for 
various  output  powers  are  listed  for  use  in  tradeoff  considerations. 


•"(.euMt.  Jf  WB  Budges  and  Jl  Panklove,  “Coherent  Optical  Sources  for  Communications,"  Proceedings 
of  the  Ih  h l . v 10.  p 1 4 Id  to  14.19,  June  1970 

/v  5 laser  hocus  Buyers' (luide.  Advanced  Technology  Publications,  Inc,  I97S 
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TABU  X.  Pl'LM  LASERS 


OPTICAL  FILTERS 


Background  radiation  noise  may  be  reduced  in  the  frequency  domain  (filtering),  the 
time  domain  (gated  receiver  in  a pulse  system),  o» spatially  (limited  FOV).  Frequency-selective 
filtering  is  considered  in  this  section  hut  the  maximum  FOV  is  generally  related  to  band- 
width tot  narrow-band  filters.  The  major  emphasis  in  this  section  is  upon  interference  fil- 
ters. which  provide  narrow  passbands  and  can  be  designed  over  a wide  range  of  optical  wave- 
lengths. Acoustically  tuned  optical  filters  are  also  discussed.  Some  other  filter  types  are 
briefly  considered. 


INTERFERENCE  FILTERS 

Interference  filters  are  commercially  available  with  typical  standard  bandwidths  of 
10,  30  and  I00A  (about  5200,  $150  and  $ 100  for  6 1 -cm  diameter).  Multilayer  dielectric- 
stack  filters  are  available  with  bandwidths  down  to  about  5A  and  solid-spacer  Fabray-Perot 
filters  have  been  developed  with  an  optical  bandwidth  of  0.5A.22 

Recently  developed  filters  in  the  infrared  provide  bandwidths  of  from  OJA  to  8A  in 
the  I to  5 pm  band.23  The  0.3A-filter  at  I pm  has  a 1-degree  allowable  FOV. 

In  its  simplest  form,  an  interference  filter  consists  of  two  partially  reflecting  surfaces 
separated  by  a spacer  with  refractive  index  n as  indicated  in  figure  34. 


1 lie  signals  will  be  enhanced  when  the  optical  path  length  is. 

2nd  = mA, 

where  in  is  an  integer  (order  of  interference).  Peak  transmissions  occur  at  wavelengths. 

Am  = 2nd/m, 


McIntyre,  C,  el  al,  Op.ical  Components  and  Technology  in  L aser  Space  Communications  Systems,” 
Proceedings  of  the  IEEE,  v 10,  p 1491  to  1503,  1970 

’Roche,  Al  and  AM  Title,  "Tilt  Tunable  Ultra  Narrow  Band  Fillers  for  High  Resolution  Infrared  Photometry,” 
Applied  Optics,  v 14,  p 765  to  769,  1975 
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so  that  multiple  transmission  bands  are  produced.  The  separation  between  transmission  bands, 
or  the  free  spectral  range,  is, 

^s  = 

One  way  of  retaining  very  narrow  bandwidths  but  increasing  the  effective  Xs  is  to  operate  2 
filters  in  tandem.  The  free  spectral  ranges  of  the  2 filters  are  chosen  so  that  undesired  coin- 
cidences occur  in  regions  where  they  can  be  reduced  by  blocking  filters  (eg,  absorption  filters). 

When  light  is  incident  on  the  filter  at  an  angle,  6 , relative  to  the  normal,  the  wave- 
length of  peak  transmission,  A0,  in  each  band  is  lowered  by  an  amount  of  6\  given  by, 

6X  = Ao02/2n2, 

for  small  0.  This  allows  the  filter  to  be  fine  tuned  by  tilting  but  also  restricts  the  FOV  for 
effective  filtering.  Figure  35  shows  a plot  of  5A/X0  versus  0 for  typical  values  of  refractive 
index.  The  figure  shows,  for  example,  that  a I A bandwidth  at  1 .06  pm  is  limited  to 
a ±2  degree  FOV.  The  values  of  n vary  with  the  general  optical  band  in  a direction  to 
reduce  the  direct  dependence  of  bandwidth  on  wavelength. 


Figure  35.  Shift  of  center  wavelength  of  typical 

interference  filters  with  angle  of  incidence. 

The  temperature  stability  of  solid-spaced  interference  filters  depends  both  upon 
the  bulk  linear  coefficient  of  thermal  expansion  of  the  substrate  and  the  temperature 
dependence  of  the  index  of  refraction. 
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I lie  magnitude  of  the  effects  depends  upon  material,  operating  temperature,  and  wavelength. 
While  temperature  control  may  be  necessary,  it  is  not  generally  a difficult  requirement.  For 
example,  in  the  near-to-mid  IK  range,  a temperature  stability  on  the  order  of  ±2  degrees  ('  is 
sufficient  for  less  than  a 1A  shift.23 


AC  OUSTICALLY  TUNFD  OPTICAL  FILTERS 

Electronically  tunable  acousto-optic  filters  (ATOF)  have  been  developed  in  the  past 
few  years.*4  : A tunable  filter  combined  with  a tunable  laser  (eg.  dye  laser),  offers  the 

possibility  ol  spread  optical  spectrum  techniques  with  resulting  countermeasures  advantages. 
Bandwidths  roughly  comparable  to  those  obtained  with  interference  filters  can  be  achieved. 
However,  the  filter  is  polarization  sensitive  and  requires  a microwave  power  source.  Inter- 
ference filters  would  lx*  preferable  unless  the  tuning  capability  is  needed. 

Figure  36  shows  the  bandwidth  as  a function  of  wavelength.  It  is  noted  that  a FOV 
of  ±3  degrees  causes  little  degradation  in  the  bandwidth.  Other  designs  (crystal  at  different 
orientation)  can  produce  very  large  off-angle  sensitivity  (eg.  a 20  times  increase  in  bandwidth 
for  a I -degree  aperture  angle). 


Harris,  SL,  STK  Nieli  and  DK  Winslow,  “Electronically  Tunable  Acousto-Optic  Filter,"  Applied  Physics 
Letter  15,  p 325,  1969 

^'Kustcrs,  JA.  DA  Wilson  and  DL  Hammond,  “Optimum  Crystal  Orientation  for  Acoustically  Tuned 
Optical  Filters,"  Journal  of  the  Optical  Society  of  America,  v 64,  p 434  to  440,  1 974. 

Naval  Electronics  laboratory  Center  Technical  Report  1 990,  Possible  l sc  of  Organic  Dyes  and  Inorganic 
Compounds  in  l iquid  as  Wide  Aperture.  Narrowband.  Optical  i nters,  bs  J A Trias  (in  preparation) 

2 ’Senitzky,  B,  “Narrowband  Ultraviolet  Vapor  Filter.  ” Applied  Optics,  v .4,  p 238  to  244,  1975 
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ABSORPTION  FILTERS 


Absorption  filters  can  be  used  with  a wide  FOV,  but  usually  have  such  a wide  band- 
width that  they  are  generally  used  only  for  gross  filtering  applications.  Recently,  however, 
a moderate  bandwidth  ( 100A ) absorption  filter  has  been  achieved  in  the  blue-green  with 
the  use  of  organic  dyes  having  a field  of  view  of  ±45  degrees.2  6 Absorption  filters  may 
also  be  used  to  suppress  unwanted  transmission  bands  in  interference  filters. 


VAPOR  FILTERS 

Vapor  filters  use  the  phenomena  of  selective  reflection  from  a vapor.  The  vapor  be- 
haves as  a reflector  at  a resonant  frequency  and  transmits  radiation  at  other  frequencies.  A 
mercury  vapor  UV  filter  at  2537  A has  been  constructed.27  The  bandwidth  of  the  filter  varies 
trom'O.  I A to  I A as  a function  of  mercury  vapor  pressure.  An  acceptance  angle  of  7 degrees 
was  achieved.  It  appears  that  vapor  filters  will  only  be  available  at  certain  specific  frequen- 
cies in  the  near  term  and  are  not  likely  candidates  for  the  present  system. 


PHOTODETECTORS 

The  three  parameters  which  define  a photodetector  are  its  responsivity  or  quantum 
efficiency,  its  noise,  and  envelope  bandwidth.  For  a pulse  system,  a sufficient  bandwidth  is 
required  for  the  detector  to  respond  to  a short  pulse.  Pulse  durations  as  short  as  about  0.01 5 
microsecond  are  considered  in  this  study.  The  advantage  of  internal  detector  gain  was  dis- 
cussed in  section  4.  Two  types  of  detectors  with  internal  gain  are  discussed  briefly  in  this 
section.  A detailed  discussion  of  photodetectors  is  given  in  Melchior.2® 


PHOTOMULTIPLIERS 

Commercially  available  photomultipliers  using  conventional  electrostatic  focusing 
have  bandwidths  up  to  about  100  MHz.  which  is  sufficient  for  the  present  application.  The 
use  of  additional  electrodes  and  crossed-fleld  techniques  can  provide  much  larger  bandwidths. 
In  tact,  a dynamic  cross-field  photomultiplier  can  be  used  as  a high-speed  gated  detector 
to  minimize  the  effect  of  background  radiation.29 

As  indicat  *d  in  figures  37,  38  and  39,  high-efficiency  photomultipliers  cover  the 
range  from  the  near  UV  through  the  visible  to  about  0.9  gm.  Typical  responsivity 
in  the  near  UV  is  0.04  to  0.07  ampere/watt  while  0.0b  to  0.08  ampere/watt  is  available  in 
the  visible  spectrum.  Quantum  efficiency  decreases  very  rapidly  above  0.9  micrometre  as 
indicated  in  figure  39.  However,  units  are  now  available  at  1.06  micrometre  with  a 2-percent 
quantum  efficiency  or  0.01 8 ampere/watt  responsivity.2 1 The  gain  of  photomultipliers 
typically  runs  up  to  about  I0ft. 


21  Melchior,  H,  MB  Fisher  and  Fit  Arams,  "Photodelectors  for  Optical  Communication  Systems," 
Proceedings  of  the  IEEE,  v 10.  p 1466  to  I486,  1970 
2 9 Leverenz,  I)J  and  OL  (faddy,  “Subnanosecond  (fating  Properties  of  the  Dynamic  Cross-Field  Photo- 
multiplier,” Proceedings  of  the  IF.  EE,  v I0,p  1487  to  1490,  1970. 
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Figure  38  Visible  photoeiniticr  characteristics  Figure  3l<.  Near  inliarcd-visible  photoemitter  characteristics 

(courtesy  RCAMrefS).  (courtesy  RCA)  (ref  8). 


The  dark  current  in  a photomultiplier  depends  upon  the  type  of  photoemissive 
surface,  cathode  area  and  temperature.  At  room  temperature,  the  range  of  dark  current 
for  different  surface  materials  varies  from  about  10“'  7 to  10"'-*  amperes  per  square 
cm.  Values  for  specific  materials  are  given  in  figure  37.  The  dark  current  decreases 
with  temperature  and  orders  of  magnitude  decreases  are  possible.  Dark-current  values 
for  photomultiplier  tubes  are  usually  given  in  amperes  rather  than  amperes  per  cm- 
and  often  are  given  as  anode  or  amplified  dark  current.  A dark  current  of  10" *4 
ampere  was  used  for  most  of  the  link  analyses.  This  value  is  appropriate  at  room  tem- 
perature for  most  of  the  spectrum,  but  10"'  ampere  (room  temperature)  is  more  typical 
at  1 .06  jam.  For  most  practical  system  parameters,  the  current  generated  by  the  day- 
time background  will  exceed  the  dark  current.  At  night  the  reverse  condition  usually 
exists. 

Commercial  grade  photomultiplier  tube  costs  in  the  visible  vary  from  less  than  one 
hundred  to  a few  hundred  dollars.  Solar-blind  UV  tubes  cost  close  to  a thousand  dollars, 
while  near  UV  costs  are  comparable  to  the  visible.  Prices  in  the  near  1R  ( 1 .06  /am ) range 
up  to  a few  thousand  dollars  but  would  be  expected  to  drop  once  volume  production  was 
obtained. 


AVALANCHE  PHOTODIODES 

The  avalanche  photodiode  is  a p-n  junction  detector  in  which  avalanche  multiplication 
is  obtained  by  operation  with  a high  reverse  bias  voltage.  These  devices  have  high  internal 
gain  and  large  bandwidth  but  have  a small  light-sensitive  area  and  generate  an  internal  noise 
which  is  a function  ot  the  gain.  The  small  junction  aTea  is  required  to  keep  the  diode  capaci- 
tance low  (for  fast  response)  and  leakage  currents  small.  The  sensitive  area  of  avalanche 
photodiodes  varies  from  about  2 X I0"-S  to  2 X IO"-?  cm-  which  makes  focusing  of 
light  onto  the  area  much  more  difficult  than  with  photomultipliers,  which  have  cathode 
areas  on  the  order  of  I cm-. 

Quantum  efficiency  as  a function  of  wavelength  for  several  types  of  photodiodes 
is  shown  in  figure  40.  Other  characteristics  of  avalanche  photodiodes  are  shown  in  table  10. 

In  addition  to  dark  current,  avalanche  photodiodes  have  a multiplicative  noise  term  which  is 
a function  of  gain.28  The  S/N  expression  in  section  4 may  be  written  to  include  this  term  as: 


S3  l*c  G ^ 

S/N  - 

2qB(Sd(Ph  + pc)+  ld)C,2F(G)+  4kTB/Rl 


where  F(G)  represents  the  multiplicative  noise.  Other  quantities  were  defined  in  section  4. 

The  factor,  F,  is  essentially  equal  to  (1  for  germanium  diodes  and  usually  equal  to  about  the 
0.4  to  0.5  power  of  G for  silicon.  The  optimum  gain  is  achieved  when  the  noise  out  of  the 
diode  (first  term  in  denominator)  equals  the  amplifier  thermal  noise.  Even  though  the  dark 
current  is  generally  high  (10"'"  to  10""  A),  the  daytime  background  current  for  X.  near  I /am 
is  expected  to  exceed  it  because  of  the  high  Sd  (about  0.5  A AV). 
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TABU  10.  CHARACTERISTIC  SOP  AVAL  ANCHE  PHOTODIODES  (REF  28). 


0.6  0.8 

WAVELENGTH  |am| 


Figure  40.  Wavelength  dependence  of  quantum  efficiency  and  responsivity 
for  several  high-speed  photodiodes  (courtesy  IEEE)  (ref  28). 

The  S/N  for  an  avalanche  diode  depends  upon  the  amplifier  noise  and  the  form  of 
F(G).  Assuming  a daytime  background,  F = G^  and  a low-noise  amplifier,  an  optimum  gain 
on  the  order  of  100  can  be  obtained.  Under  these  assumptions,  the  resulting  S/N  would  be 
just  slightly  better  than  the  2-percent  quantum  efficiency  multiplier  at  1 .00  pin.  The  choice 
between  photomultipliers  and  avalanche  diodes  requires  a detailed  consideration  of  practical 
design  problems.  Silicon  avalanche  photodiodes  are  in  the  price  range  between  $200  and 
$2000. 30 

POINTING  AND  TRACKING  SYSTEMS. 

Ship-to-ship  and  ship-to-RPV-to-ship  optical  communication  links  require  pointing 
and  tracking  systems.  Furthermore,  as  has  been  demonstrated  in  the  Naval  Electronics 
Laboratory  Center  OCCULT  experiments,  high-rate  reciprocal  pointing  and  tracking  (in 
which  both  terminals  in  a link  track  each  other)  can  produce  orders  of  magnitude  improve- 
ment in  S/N  and  acquisition  speed  over  less  sophisticated  system  ;.  The  purpose  of  this  sub- 
section is  to  briefly  summarize  what  is  currently  available  in  such  systems. 

For  shipboard  use,  the  OCCULT  system  represents  the  current  state  of  the  aTt.  Its 
specifications  call  for  pointing  stability  to  a given  point  with  peak  errors  less  than  + 0.2 
degree  (±  3.5  mrad)  on  a platform  experiencing  20-degree  rolls  with  10-second  periods. 


10  Medved.  D,  “Photodiodes  for  Fast  Receivers,”  /.aser  Focus,  v 10,  p 45  to  47,  1974 
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It  will  be  capable  of  2-axis  (train  and  elevation)  tracking,  at  rates  of  20  degrees  per  second, 
either  axis.  Upon  receipt  of  an  acquisition  pulse,  it  will  stop  searching  and  return  to  the 
correct  position  ± 0.2  degree  in  <0.2  seconds,  using  a 90  degree  per  second  search  rate.  The 
system  is  currently  in  its  prototype  test  and  development  stage,  with  some  shipboard  tests 
having  been  carried  out.  It  appears  possible  that  the  pointing  accuracy  specifications  may 
be  exceeded  by  a factor  of  two.  In  its  prototype  form,  it  requires  2 metres  of  standard  rack 
space  tor  the  control  electronics,  plus  <0.1  cubic  metre  for  the  gimbal  and  mount.  System 
hardware  cost  is  about  S35000. 

Currently  available  systems  for  KPV  use  are  typified  by  an  Army  system  which  is 
par!  of  a video  tracker/laser  designator  payload  weighing  a total  of  14  kg.  It  provides  a 
gyrostabilized  platform  with  tracking  in  2 axes,  azimuth  and  elevation.  It  tracks  at  a 
rate  of  10-degrees  per  second,  with  360-degree  azimuth  and  15-degree  elevation  move- 
ment. Pointing  accuracy  is  not  known;  however,  for  target-designation  purposes  (the  Army 
mission)  accuracy  of  <17  mrad  ( 1°)  would  probably  be  required.  This  system  has  under- 
gone operational  tests. 


REMOTELY  PILOTED  VEHICLES 

The  purpose  of  this  section  is  to  characterize  the  current  and  near-term  future 
capabilities  of  remotely  piloted  vehicles  (RPV)  for  use  as  central  platforms  in  OTH  optical 
communication  links.  Current  RPVs  range  from  1 1 kg  and  2.4  m in  wingspan  to  a remotely 
piloted  f'- 1 02 A aircraft.  Operational  requirements  for  the  Navy  Fleet  and  probable  system 
physical  specifications  for  active  repeater  platforms  for  use  with  OTH  optical  links  imply 
that  the  smaller  RPVs  (often  termed  “mini-RPVs”)  are  adequate.  Figure  41  shows  the  mini- 
mum platform  altitude  (km).  111.  required  to  achieve  LOS  to  the  platform  over  a distance, 
I),  km,  given  a transmitter  at  altitude  11]  (km).  Note  that  the  maximum  total  communica- 
tion distance  is  twice  D. 


figure  41.  Minimum  platform  altitude  required  to  achieve  line-of-sight. 


65 


The  ideal  RPV  would  have  a radar  and  optical  cross  section  of  zero,  a payload  of 
about  ‘>0  kg.  a range  of  150  km.  be  launched  and  recovered  on-board  ship,  have  an  alti- 
tude range  from  0 to  5 km.  and  be  able  to  loiter  on  station  as  long  as  desired.  Although 
extreme,  most  of  these  specifications  are  at  least  approachable  within  the  next  5 years, 
based  on  existing  RPVs  and  associated  technology,  fable  I I summarizes  existing  small 
RPVs  and  those  under  development.  Of  special  interest  for  OTH  optical  communication 
applications  is  the  Teledyne  Ryan  mini-RPV,  called  "STAR"  (Ship  Tactical  Airborne 
RPV)  I his  RPV  was  scheduled  for  flight  tests  in  February  1 070.  It  is  to  be  launched 
from  a ship  via  pneumatic  rail  and  recovered  on-board  by  landing  it  in  a horizontal  net.  It 
can  operate  either  in  a preprogrammed  mode,  thus  eliminating  the  necessity  of  a command  link 
to  the  RPV  (and  increasing  the  covertness  of  the  RPV),  or  in  a fully  controllable  command 
mode.  The  payload  weight  can  vary  from  12  to  about  22  kg.  With  a 1 2-kg  payload, 
the  RPV  can  fly  at  a speed  of  255  km  per  hour  to  a range  of  225  km.  loiter  for  3 hours 
at  a speed  of  1 13  km  per  hour,  and  return.  Shorter  ranges  would  provide  longer  loiter 
times.  As  presently  configured,  STAR  has  a flight  ceiling  of  about  2.5  km  (which,  from 
figure  41.  is  consistent  with  its  flight  range  of  200  to  225  km  for  a reasonable  loiter  time). 

The  current  STAR  payload  contains  a video  tracker  and  laser  designator  in  a 25-cm 
diameter,  hemispherical,  optical  dome.  The  package  is  gyrostabilized.  The  estimated  life- 
time is  10  nights,  with  the  limiting  factors  being  recovery  and  pilot  error  (although  the 
RPV  has  notation)  The  engine  requires  maintenance  every  100  flight  hours.  The  cost 
of  such  an  RPV  is  heavily  dependent  upon  the  payload.  With  the  video  tracking/laser 
designator  payload,  a cost  of  SI  50  000  per  vehicle  in  quantities  of  50  is  estimated.  This 
could  easily  drop  to  S200  000  in  quantities  of  1000.  without  payload. 

The  optical  and  radar  cross  sections  of  mini-RPVs  are  small.  Physically,  the 
human  eye  cannot  resolve  an  object  the  size  of  an  RPV  beyond  a range  of  about  25  km 
(and  it  is  probable  that  the  actual  detection  range  is  less  than  one  fifth  of  this).  Non- 
human  optical  detectors  can  be  made  to  have  even  greater  range  limitations  by  the  proper 
choice  of  RPV  surface  finish.  The  radar  cross  section  can  be  of  the  order  of  two  seagulls 
for  an  RPV  using  plastic  structural  materials  and  skin  (as  is  done  in  STAR  ) 

In  summary,  current  RPV  technology  appears  capable  of  providing  a shipboard- 
launched-and-recovered  platform  which  would  not  limit  an  OTH  optical  communications 
link,  except  for  the  time  it  can  remain  on  station,  assuming  a payload  capacity  of  about 
I I -kg  is  adequate.  It  appears  within  current  and  near-term  future  technological  capability 
to  produce  a shipboard  launched  and  recovered  RPV  with  an  on-station  endurance  of 
()  to  |()  hours  and  a payload  capacity  of  23  kg.  The  other  characteristics  of  such  an  RPV 
would  be  comparable  with  STAR. 
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SECTION  7.  TRADEOFF  CONSIDERATIONS 


The  level  of  effort  in  this  study  has  precluded  a rigorous  system  optimization  pro 
eedure  of  the  type  described  hy  Pratt.32  Such  an  analysis  might  be  premature  anyway, 
especially  for  the  scatter  link,  since  the  system  performance  depends  so  strongly  upon  the 
link  characteristics.  Taking  advantage  of  propagation  characteristics  can  be  expected  to  give 
improvements  in  performance  which  can  be  much  mo  a-  significant  than  the  difference  be- 
tween optimum  and  a>asonable  component  combinations,  l or  example.  Hying  an  active 
RPV  at  a 2.5  km  altitude  rather  than  at  the  minimum  LOS  altitude  gives  a predicted 
30-dB  S/N  improvement  at  a range  of  ‘>0  km  for  a visibility  of  20  km. 

Since  the  link  analysis  indicated  that  the  desired  150  nmi  range  could  not 
usually  be  achieved  under  typical  daytime  conditions,  the  “optimization"  procedure  (sec- 
tion 3)  was  to  select  approximately  the  best  state-of-the-art  component  parameters  and 
determine  system  performance  for  both  links.  An  example  of  the  effect  of  system  para- 
meters on  the  S/N  (see  section  4)  is  given  in  table  1 2.  It  is  assumed  that  the  detector  has 
sufficient  internal  gain  that  the  noise  contribution  of  the  following  amplifier  has  negligible 
effect  on  S/N.  Note  that  peak  power  and  pulse  rate  can  usually  be  traded  off  under  a fixed 
average  power  constraint.  The  relations  given  in  table  I 2 can  be  used  for  rapid  evaluation 
of  tradeoffs  and  then  the  appropriate  program  in  the  appendix  used  for  an  accurate  S/N 
determination. 

Costs  for  the  shipboard  system  are  dominated  by  the  laser  and  pointing  and  track- 
ing subsystem.  This  assumes  that  expected  economies  are  achieved  in  plastic  optic  s.  For 
the  RPV  relay,  the  laser  and  pointing-control  costs  are  expected  to  be  significant. 

A tradeoff  or  selection  of  one  of  the  link  types  over  the  others  does  not  seem 
appropriate  at  this  time.  As  one  progresses  from  a scatter  link  to  an  RPV  system,  both  the 
cost  and  the  communication  range  increase.  Both  could  use  the  same  shipboard  system 
and  be  interchangeable  as  environmental  conditions  change. 


12  Pratt,  WK,  LS  Stokes  and  R Hinckley,  “Optimization  of  Optical  Communication  Systems,"  Proceedings 
of  the  IEEE :.v  10.  p 1737  to  1740.  1970 
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TABLE  12.  EFFECT  OF  COMPONENT  PARAMETERS  ON  S/N  OF  A SPECIFIC  SYSTEM 


Fixed  Parameters 

Value 

Wavelength 

1 .06  nm 

Detector  dark  current 

io- 

A/W 

Internal  detector  gain 

(see  text) 

Variable  Parameter 

Condition 

S/N  effect  proportional  to 

Peak  power,  P 

Signal  power»background 

10  log  P 

Signal  power«background 

20  log  P 

Responsivity,  Sj 

Daytime 

10  log  Sj,  10  log  A 

Aperture  A 

Night,  large  received  signal 

10  log  Sj,  10  log  A 

Night,  small  received  signal 

20  log  Sj,  20  log  A 

Optical  Bandwidth.  B() 

Daytime 

-10  log  ^ 

Night 

no  effect,  reasonable  B,, 

Transmitter  divergence. 

Scatter  link 

Complicated  (see  section  5), 

®D 

small  effect  for  small  0[y 

Line-of-sight  link 

-20  log  0q 

Receiver  FOV 

Scatter  link,  background 

no  effect  in  vertical  plane 

limited 

-10  log  FOV  horizontal  plane 

Line-of-sight  link,  back- 

-20  log  FOV 

ground  limited 

SECTION  S.  COVERT  CONSIDERATIONS 

Covert  considerations  have  been  alluded  to  throughout  this  report.  This  section 
briefly  summarizes  general  factors  which  affect  covertness.  A more  quantitative  analysis  is 
recommended  before  system  implementation. 

Because  of  their  extremely  narrow  divergence  or  beamwidth.  laser  transmissions  are 
difficult  to  detect.  Jamming  is  also  difficult  because  of  the  narrow  FOV  usually  used  in  re- 
ceivers. However,  the  aerosol  scatter  characteristics  present  a limit  on  both  covertness  and 
antijamming:  signal  energy  is  scattered  out  of  the  beam  and  jamming  energy  is  scattered 
into  the  receiver.  For  this  reason,  if  no  other,  the  scattering  characteristics  of  marine  hazes 
-•hould  be  well  defined.  In  the  long  term,  spread-spectrum  techniques,  using  tuneable  lasers 
and  filters,  should  be  incorporated  to  enhance  the  covertness  of  optical  systems. 

From  an  intercept  standpoint,  the  transmitted  power  should  be  adjusted  to  be  just 
sufficient  to  support  the  required  data  rate.  However,  a reserve  of  transmitter  power  is  de- 
sirable to  combat  jamming.  With  this  power  limitation,  an  intercept  receiver  could  not  be- 
at a much  greater  range  than  the  communication  link  For  a LOS  link,  the  intercept  re- 
ceiver would  have  to  be  at  a much  shorter  range  unless  it  was  on-axis  (or  nearly  so)  with  the 
ire  am.  Even  for  a scatter  channel,  being  off-axis  restricts  the  range  because  of  the  scattering 
function.  Scattered  power  is  typically  down  by  10  dB  at  30  degrees  off  axis,  decreasing  to 
-20  dB  at  65  degrees  anil  -30  dB  at  1 10  degrees,  for  marine  aerosols.  In  addition,  the  angular 
source  size  of  the  beam  becomes  much  larger  off-axis  which  requires  a very  large  receiver 
FOV,  significantly  increasing  background  and  pulse  dispersion.  C learly,  this  must  be  further 
characterized  both  experimentally  and  theoretically. 
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SECTION  9 CONCLUSIONS  AND  RECOMMENDATIONS 


CONCLUSIONS 

An  analytical  model  for  OTH  laser  propagation  was  developed,  based  upon  the  single 
scattering  approximation.  The  model  appears  to  be  in  reasonable  agreement  with  limited  field 
data  on  OTH  optical  propagation  at  relatively  close)  ranges  up  to  40  km.  A more  conclusive 
field  test  of  the  model,  particularly  at  longer  ranges,  would  be  necessary  to  accurately  assess 
its  utility  for  OTH  optical-scatter  communications  system  design  Based  upon  the  use  of  sys- 
tems composed  of  state-of-the-art  (1976)  components,  the  following  conclusions  were  draw  n 
from  this  analysis. 

1.  The  operating  wavelength  should  be  in  the  l-to-3  pm  range  for  all  links; 

2.  Both  relay  and  scatter  OTH  data  links  can  use  the  same  shipboard  system; 

3.  Significant  performance  advantages  can  be  achieved  by  exploiting  propagation 
characteristics.  These  include: 

(a)  Use  of  a vertical  fan  beam  or  optimally  elevated  beams  in  a scatter  link,  and 

(b)  Positioning  the  relay  platform  at  a high  altitude  to  take  advantage  of  the 
decrease  in  path  loss  caused  by  the  vertical  decrease  in  aerosol  concentration; 

4.  The  primary  factor  which  determines  the  percentage  of  time  communication  over 
a given  distance  at  a given  bit  rate  can  be  achieved  is  the  meteorological  visibility.  Statisti- 
cal studies  of  the  occurrence  of  visibilities  greater  than  a given  value  are  available.  The  accumu- 
lated statistics  were  analyzed  to  determine  the  general  availability  of  an  optical  scatter-link 
communications  capability  having  bit  rates  of  2400  and  75  bits  per  second  at  a bit-error  rate 
of  less  than  10“^.  When  the  statistics  are  summed  for  all  ocean  areas  the  following  results  are 
obtained  for  daytime: 

Availability  2400  bits/sec  75  bits/sec 

b5%  55  km  85  km 

100  km  165  km 

The  85-percent  availability  isequatable  to  a meteorological  visibility  condition  of  10  km  and 
the  70-percent  availability  to  20-km  visibility.  Higher  availabilities  and  longer  ranges  can  be 
achieved  with  improved  meteorological  visibilities.  For  example,  in  the  Mediterranean,  the 
communication  ranges  listed  for  a visibility  of  20  km  can  be  achieved  with  an  availability  of 
95  percent;  and 

5.  The  communications  ranges  are  considerably  greater  (by  factors  of  2 to  3)  at  night, 
or  if  low  clouds  occur  to  provide  a scattering  layer,  it  is  important  to  note  that  areas  of 
operation  wnth  lower  average  visibility  (eg,  the  North  Atlantic)  also  have  a high  occurrence  of 
low  cloud  cover.  The  interpretation  of  the  visibility  data  on  range,  data  rate,  and  availabili- 
ty requires  further  characterization. 

Performance  predictions  for  state-of-the-art  systems  for  a number  of  conditions  are 
given  in  tables  3 and  4 of  section  3.  .Section  3 discusses  the  systei.i  concepts  and  gives  a more 
detailed  summary  than  that  presented  here. 
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RECOMMENDATIONS 


The  forward-scatter  link  analysis  presented  in  this  report  is  based  upon  a single  scat- 
tering model.  While  the  model  gives  reasonable  agreement  with  the  limited  field  data,  it  is 
strongly  recommended  that  experiments  be  designed  and  carried  out  to  assess  the  validity  of 
the  model.  The  experiments  might  well  include  simultaneous  pulsed  measurements  at  two 
wavelengths  tie,  I .()(>  and  0.53  /Jin)  since  the  analysis  predicts  large  differences  as  a function  of 
wavelength  The  performance  gain  achieved  by  the  use  of  elevated  beams  and  exploitation 
of  clouds  should  be  evaluated.  It  is  recommended  that  fading  rate  and  other  data  be  ob- 
tained to  allow  evaluation  of  time,  spatial,  and  receive-angle  diversity  receiver  techniques. 

It  is  recommended  that  the  development  of  economical  pointing  and  tracking  sub- 
systems be  continued,  since  this  is  one  of  the  major  cost  items  in  a shipboard  system.  On  a 
long-term  basis,  the  development  of  tunable  lasers  in  the  l-um  region  should  be  encouraged 
so  that  spread-spectrum  techniques  could  be  used.  Tunable  filters  are  currently  available. 

It  is  also  recommended  that  careful  consideration  be  given  to  the  potential  economies  and 
resulting  increase  in  practical  system  performance  levels  which  may  be  realized  with  the 
development  of  low-cost  plastic  optical  systems  for  communications  use. 

Because  a primary  factor  determining  link  availability  is  that  of  the  meteorological 
conditions  prevailing  over  the  link,  it  is  strongly  recommended  that  experiments  be  carried 
out  to  verify  the  meteorological  effects  upon  propagation  predicted  by  the  models.  It  is 
also  strongly  recommended  that  an  integrated  and  systematic  capability  for  measuring  and 
predicting  the  meteorological  conditions  predominating  over  a path  be  developed.  This 
might  well  take  advantage  of  a combination  of  optical  measurements  and  existing  micro- 
wave  meteorological  measurement  capabilities. 


APPENDIX  A:  INCORPORATION  OF  THE  PROPAGATION 
MODELS  IN  THE  COMPUTER  PROGRAMS. 


When  work  on  the  task  covered  by  this  report  began  it  was  thought  that  existing  Monte 
Carlo  computer  codes  would  be  employed  to  produce  the  numerical  results  required  for  the 
system  tradeoff  studies.  This  was  not  done  for  several  reasons.  First,  an  analytical  model  for 
the  over-the-horizon  (OTH)  case  was  developed  (see  appendix  A)  which  has  resulted  in  several 
important  insights,  and  considerable  simplification  in  characterizing  practical  communications 
links.  It  was  felt  that  it  was  important  to  incorporate  this  model  into  a computer  program  in 
order  to  exploit  these  advantages.  Second,  the  Monte  Carlo  codes  could  not  model  the  OT1I 
case  without  modifications  which  would  not  be  completed  in  the  time  frame  necessary  to 
meet  task  milestones.  Third,  the  line-of-sight  (LOS)  propagation  mode  (eg,  using  an  RPV) 
could  be  adequately  modeled  without  resorting  to  the  cost  and  complexity  of  using  the  Monte 
Carlo  codes. 

Both  the  OTH  and  LOS  propagation  models  were  incorporated  into  interactive  com- 
puter programs  written  in  the  BASIC  language,  for  use  on  a minicomputer.  The  result  is  a tool 
which  permits  the  scientist  and  system  engineer  to  explore  both  the  propagation  models  and 
potential  system  designs  interactively,  with  high  efficiency.  The  benefits  which  have  resulted 
from  this  effort  have  far  outweighed  the  unanticipated  additional  effort.  This  appendix  de- 
scribes how  the  OTH  and  LOS  models  were  incorporated  into  these  programs,  including  the 
approximations  used.  The  programs  themselves  are  described  in  appendix  C (OTH)  and  ap- 
pendix 1)  ( LOS),  along  with  user  instructions  and  listings.  Because  this  programming  effort 
was  not  originally  planned,  an  effort  has  been  made  to  keep  the  documentation  as  brief  as 
possible,  consistent  with  providing  enough  information  for  other  researchers  to  make  use  of 
the  programs. 

The  OTH  model  is  described  in  detail  in  appendix  E.  However,  in  order  to  use  it,  cer- 
tain parameters  are  required.  The  extinction  coefficient,  /3,  must  be  calculated  as  a function 
of  wavelength,  X,  and  meteorological  range  (also  termed  visibility).  V.  As  discussed  by  Wood- 
man33. the  most  widely  used  expressions  for  calculating  /3(X.  V)  underestimate  /3  by  a factor 
of  from  2 (at  X = 1/a.  V = 23.5  km)  to  greater  than  10  (at  X = 10p.  V = 23.5  km),  with  much 
higher  errors  for  smaller  values  of  V.  Note  that  this  underestimates  the  attenuation  by  fac- 
tors of  from  5 to  >5  X l(P.  Flterman34  lias  done  an  extensive  study  of  available  models 
and  experimental  results,  and  derived  an  expression  for  /3(  V,  X)  which  is  consistent  and  in 
good  agreement  with  experimental  results  for  the  range  .27/a<  X <2. 1 Ip.  for  V>  2 km.  Sur- 
face aerosol  distributions  characteristic  of  hazes  with  2 km  < V < 13  km  are  assumed. 
Flterman34  considers  the  effects  of  such  parameters  as  aerosol  distribution,  mixing  layers, 
scale  height,  and  the  approximations  used  in  deriving  his  results  in  great  detail,  none  of  which 
will  be  repeated  here.  However,  it  is  clear  that,  for  the  system  engineering  purposes  of  the 
work  described  in  this  report,  the  Flterman  model  is  quite  adequate  for  the  above  range  of 
wavelengths  and  for  V<13  km,  which  Flterman  takes  as  the  boundary  between  haze  and 
clear  conditions.  The  numerical  results  presented  in  this  report  include  values  for  V as  high 
as  390  km.  The  Flterman  formula  was  used,  without  modification,  for  these  values  of  V. 


' ’Woodman,  Douglas  P,  “Limitations  in  Using  Atmospheric  Models  for  Laser  Transmission  Estimates," 
Applied  Optics,  v 1 3,  p 2193  to  2195,  1974 

34Elterman,  L,  “Relationships  Between  Vertical  Attenuation  and  Surface  Meteorological  Range."  Applied 
Optics  v 9,  p 1804  to  1810,  1970 
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This  was  felt  to  be  justified  for  two  reasons;  first,  due  to  the  methods  used  by  Elterman,34  it 
should  produce  conservative  values  for  0 for  V >13  km;  second,  the  values  lor  the  extinction 
coefficient  for  13  km  < V < 336  km  (the  Rayleigh  limit)  at  X = 0.55/a  agree  well  with  the 
values  given  by  the  Elterman  formula  for  these  values  of  V.  The  Elterman  formula  is 


0<V.  X)  = 


0(4.  X)  H3.91/V)-0rl  . 
[(3.91/4)  - 0r] 


(1) 


V is  in  kilometers  and  0 and  0r  are  km-'  The  formula  is  normalized  to  V = 4 km,  hence  the 
numerical  factors  of  4.  0r  is  the  Rayleigh  attenuation  coefficient  at  X = 0.55/z.  0r  = 1162  X 
10~“  km-' . The  computer  programs  store  a table  of  values  (given  by  Elterman)  for  0(4,  X). 
Linear  interpolation  on  the  logarithms  of  the  tabular  values  as  a function  of  the  logarithm 
of  the  wavelength  is  used  to  find  0(4,  X)  for  the  particular  value  of  X desired,  and  equation 
( 1 ) then  gives  0(  V.  X).  The  use  of  log-log  interpolation  was  arrived  at  empirically  after  trying 
various  polynomial  interpolation  methods.  This  method  of  interpolation  produces  values  of 
0 for  intermediate  values  of  X which  are  at  least  as  accurate  as  the  model. 

Elterman  gives  values  of  0(4,  X)  only  for  0.27/a<  X < 2 1 7/a.  To  extend  these  values  to 
10.2/a,  some  theoretical  values  of  0(23.5,  X),  computed  by  Bamhardt  and  Street  and  given 
by  Woodman  < B 1 ),  were  used.  These  values  were  adjusted  to  a meteorological  range  of  4 km 
using  equation  ( I ).  These  theoretical  values  agreed  with  Elterman ’s  values  where  the  range 
of  wavelengths  involved  overlapped.  There  is  no  limitation  in  the  theory  used  to  derive  equa- 
tion ( 1 ) which  would  preclude  its  use  at  10. 2/a. 

The  definition  and  choice  of  the  equation  for  the  signal-to-noise  ratio  is  discussed  in 
detail  in  section  4 of  the  basic  report.  As  used  in  the  computer  programs,  the  signal-to- 
noise  ratio  (S/N)  is  given  by 


S/N  = S3  P J / |2qB[SU(Pb  + Pc)  + Id  ]}  . (2) 

Sq  is  the  detector  responsivity,  given  by 

Sj  = T?q/lm,  (3) 

where 


rj  is  the  detector  quantum  efficiency, 
q is  the  electronic  charge  (couloumbs), 
h is  Planck’s  constant  (joule-sec), 

v is  the  carrier  frequency  (Hz), 

Pc  is  the  received  carrier  power  (watts). 

B is  the  electrical  bandwidth  of  the  system  (Hz), 

Pb  is  the  background  power  incident  on  the  detector 
(watts),  and 

Iq  is  the  detector  dark  current  (amperes). 
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The  model  for  the  background  noise  power  is  discussed  in  section  5 of  the  basic 
report.  As  implemented  in  the  computer  programs. 


= PQ  +4.75  X lo/  (x'3/2  [exp  ( 2.44/X )—  1 || 

, daytime. 

(4) 

= P0  +8.3  X 10-3/  jx ' -7/2  [exp  (2.44/X H || 

, nighttime. 

(5) 

PQ=  0,X<  1.8/i 

= 1.2  X 108/|x5  | exp  < 50/X)— 1 )j  , X>1.8p6. 

The  neglect  of  the  PQ  term  avoids  a numerical  problem  injhe  computations  without  altering 
the  results.  In  these  equations,  X is  in  micrometres,  and  P^  is  in  W m~~  Sr-'  p-'.  In  other 
words,  P^  is  the  background  spectral  radiance.  It  thus  is  necessary  to  multiply  Ph  by  the 
receiver  aperture  area,  the  receiver  optical  bandwidth,  and  the  receiver  solid  angle  FOV  in 
order  to  obtain  Pjv  Although  the  aperture  area  and  optical  bandwidth  are  specified,  care  must 
be  exercised  in  the  value  used  for  the  solid  angle  intercepted  by  the  receiver  field  of  view.  It 
can  be  shown  that  for  a circularly  symmetric  FOV  of  full  angle  0.  the  solid  angle  intercepted  is 

2ir  [1- cos  (0/2)1.  (7) 

For  small  0 this  reduces  to  tt0-/4.  In  the  OTH  case,  it  is  reasonable  to  expect  that  the  receiver 
FOV  may  be  broader  vertically  (full  angle  0j ) than  horizontally  (full  angle  0^)  (see  appendix 
Ft.  This  permits  accommodating  the  vertical  fan  beam,  while  minimizing  the  background  power 
received  horizontally.  As  an  approximation,  the  solid  angle  intercepted  by  the  receiver  field 
of  view  for  computing  the  background  power  is  calculated  as  7r0|02/4  for  the  OTH  case.  A cir- 
cularly symmetric  FOV  is  assumed  for  the  line-of-sight  (LOS)  case. 

For  the  OTH  ease,  it  is  necessary  to  know  f (0)  for  0=0°,  where  f (0)  is  the  gain  over 
isotropic  for  aerosol  scattering  at  the  scattering  angle  0 ( ie,  the  ratio  of  the  actual  scattering 
cross  section  at  angle  0 to  that  of  an  isotropic  scatterer  ot  identical  total  scattering  cross 
section).  This  is  available  in  Deirmendjian’s  hook35  as  a function  of  wavelength  and  aerosol 
size  distribution.  Values  for  two  distributions,  a continental  haze  “L  ’ and  a marine  haze  M . 
are  stored  as  a function  of  wavelength  by  the  computer  programs,  and  “log-log”  interpolation 
is  used  to  generate  f (0°)  for  a particular  wavelength. 


35  Dcirmendjian,  D,  Electromagnetic  Scattering  <»)  Spherical  Polydhpersions,  American  I lsevier  Publishing 
Company,  Inc,  I ‘*69 
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APPENDIX  B OTH  PROGRAM  DESCRIPTION 

Hie  Ol  H propagation  model,  as  described  in  appendix  E,  is  incorporated  into  an 
interactive  program.  Ibis  program  is  designed  to  permit  the  user  to  choose  the  range  and 
step  size  ot  the  independent  variable,  to  easily  vary  parameter  values,  and  to  print  out  the 
path  loss  and  signal-to-noise  ratio  as  a function  of  the  independent  variable  in  tabular  form. 

I he  program  will  be  described  by  referring  to  the  line  numbers  of  the  various  program 
modules  and  their  subsets.  A sample  console  listing  and  output  precede  the  program  listing. 

I he  program  is  written  in  BASIC,  and  has  been  run  on  a Data  General  NOVA  800  computer. 
It  should  execute  without  modification  on  any  other  manufacturer’s  computer  under  a 
BASIC  interpreter  which  accepts  standard  Dartmouth  BASIC.  The  only  exception  to  this 
statement  is  the  use  ot  a form  ol  the  INPUT  statement  which  includes  a string  of  text  to  be 
printed  prior  to  the  question  mark  which  requests  a user  response.  These  statements  can  be 
replaced  by  a PRINT  statement  and  the  standard  INPUT  statement  if  required. 

The  computation  portion  of  this  program  is  marginally  compute  bound  for  the  pen- 
cil beam  case  and  strongly  compute  bound  for  the  fan  beam  case  (using  a teletype  as  the  out- 
put device).  For  comparison,  compute  bound  FORTRAN  programs  tend  to  execute  on  the 
NOVA  800  at  from  I to  3 times  the  speed  of  execution  for  the  same  program  (level  G 
compiler)  on  an  IBM  360/67  running  under  OS.  Of  course,  a compute  bound  BASIC  pro- 
gram will  execute  on  the  order  of  a factor  of  ten  times  slower  on  the  NOVA  800  than  a 
FOR  TRAN  program  doing  the  same  computations.  The  response  time  of  the  interactive 
program  is  quite  adequate. 

HM50:  Set  up  parameter  and  variable  values. 

This  section  of  the  program  initializes  the  variable  and  parameter  values.  The  arrays 
dimensioned  in  statements  30  and  35  are  used  as  follows: 

W,  B Wavelength  and  associated  value  (respectively)  of  the  extinction  coefficient 

at  X = 0.55p,  V = 4km.  This  forms  the  table  used  by  the  interpolation  rou- 
tine to  compute  |?(V,A)  using  Elterman’s  model  (see  appendix  A). 

X,  V Temporary  arrays  used  as  the  working  storage  of  the  interpolation  routine. 

P Stores  the  parameter  values. 

N.  M Wavelength  and  associated  value  (respectively)  of  the  forward  scattering  coef- 

ficient, f (0),  for  a marine  haze  aerosol  distribution. 

O,  L 1 he  same  as  N,  M,  but  for  a continental  haze  aerosol  distribution. 

V Stores  the  three  possible  independent  variable  values. 

The  subroutine  at  4130  is  called  to  till  the  tables  used  for  interpolation  and  to  put  the  default 
parameter  values  into  P.  The  aerosol  distribution  to  be  assumed  is  chosen,  and  the  fan  or 
pencil  (0  divergence)  beam  case  is  specified.  For  the  fan  beam  case,  the  default  beam  angle 
is  printed  as  a reminder;  this  is  the  angle  over  which  the  received  power  is  averaged.  For  the 
pencil  beam  case  the  default  receiver  vertical  field  of  view  value  is  printed.  For  both  cases, 
the  program  integrates  the  received  power  over  this  field  of  viev  If  the  user  decides  to  change 
any  of  the  default  parameter  values,  the  subroutine  at  statement  4570  is  called. 
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There  are  three  possible  independent  variables:  wavelength,  path  length,  and 
meteorological  range.  The  user  enters  the  starting  value  for  the  desired  independent  vari- 
able, values  for  the  other  two.  specifies  which  of  the  three  is  the  independent  variable  and 
then  enters  the  increment  and  number  of  independent  variable  values  to  be  looped  over. 
F-inally,  the  choice  is  made  between  day  or  night  background  power  values  to  be  used  in 
computing  the  signal-to-noise  ratio. 


490-630  Begin  the  loop  over  the  independent  variable,  and  compute  path  loss 

This  series  of  statements  controls  calls  to  the  various  subroutines  necessary  to 
compute  the  path  loss.  At  statement  630,  P5  is  the  path  loss,  defined  as  follows: 


P5 


rec 


^trans  \ec  * 


(Bl) 


where  Prec  is  the  received  power,  Ptrans  is  t,K‘  transmitted  power,  Arec  is  the  area  of  the 
receiver  aperture  in  in-,  and  f (0°)  is  the  forward  scattering  coefficient  (see  appendix  A). 

635-971:  C omputation  of  signal-to-noise  ratio 

Depending  on  the  aerosol  distribution  chosen,  the  scattering  coefficient  is  calcu- 
lated in  lines  650-830.  The  background  power  is  computed  in  lines  850-950.  The  test 
on  the  value  of  X at  line  872  and  its  subsequent  results  avoid  an  underflow  problem  at 
line  875.  This  does  not  affect  the  accuracy  of  the  result.  At  line  971,  SI  is  the  signal-to- 
noise  ratio,  defined  as  given  in  equation  2)  of  appendix  A. 

1000-1050:  Output 

As  currently  written,  line  1000  prints  the  value  of  the  independent  variable,  the 
path  loss,  the  path  loss  in  dB,  the  signal-to-noise  ratio,  and  the  signal-to-noise  ratio  in 
dB.  The  program  then' gives  the  user  the  option  of  terminating  execution  or  initiating 
another  run 

2500-2590.  Integration  over  field  of  view 

This  subroutine  integrates  the  path  loss  over  the  receiver  field  of  view,  the  result 
is  Z3.  The  variables  are  as  follows: 

HI,  112  Height  of  transmitter,  receiver  in  km 
Rl  4/3  earth  radius,  in  km 

DI  propagation  range,  in  km 

F2  vertical  field  of  view  of  the  receiver,  rads. 

The  integration  loop  is  from  2550  to  2580.  The  choice  of  0.003  radian  as  the 
integration  step  size  is  arbitrary,  a choice  that  provides  adequate  accuracy  for  most  practi- 
cal fields  of  view  while  maintaining  reasonable  response  time  for  the  program.  For  a much 
longer  run  in  which  execution  time  can  become  quite  significant,  the  step  size  can  be 
chosen  as  F2/5.  This  provides  accuracy  consistent  with  the  accuracy  of  the  model.  The 
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number  of  iterations  through  this  loop,  and  the  number  of  iterations  through  the  loop 
which  starts  at  line  2620,  control  the  primary  execution  speed  versus  accuracy  tradeoff 
possible  for  the  programs  in  this  and  the  following  appendix. 

2600-3220  Average  over  fan  beam  and  compute  path  loss 

This  routine,  and  the  short  auxiliary  routine  in  lines  3500-3520,  compute  the  path 
loss,  Z5,  as  given  in  equation  Bl ).  The  additional  variables  not  already  defined  are  as 
follows: 


PI,  P2 

Azimuthal  pointing  angles  (transmitter,  receiver)  in  radians,  measured 
relative  to  the  great  circle  path. 

B0 

the  extinction  coefficient  in  km-' 

110 

the  aerosol  lapse  rate,  in  km-'  (the  aerosol  density  as 
height,  h,  is  proportional  to  exp  (— h/H0)) 

a function  of 

Al,  A2 

elevation  angles  (transmitter,  receiver)  of  the  beam  in 
to  the  local  horizontal,  measured  positive  downwards. 

radians,  relative 

FI 

vertical  fan  angle  of  the  transmitter,  in  radians 

The  details  ot  the  computations  of  this  routine  can  be  deduced  by  comparison  with  the 
equations  in  appendix  E.  See  also  the  comments  under  “Caution  " 

4000-4450:  Default  model  parameters 

The  parameters  (stored  in  array  P)  their  default  values,  units,  and  the  names  of 
the  variables  to  which  they  are  equated  (if  other  than  the  P array)  are  listed  in  lines  4030- 
4170.  The  units  used  for  these  parameters  were  chosen  to  facilitate  interactive  input. 

4460-45  50:  Equate  Variables  and  Convert  Units 

This  routine  equates  the  P array  values  to  the  variables  used  in  the  other  routines, 
converting  units  as  required.  Note  that  some  parameters  are  used  directly  from  the  P 
array.  Unit  conversion  for  these  values  is  carried  out  where  they  are  used. 

4560-4710:  Changes  parameter  values  interactively 

6000-6170:  Interpolation 

This  routine  performs  “log-log”  interpolation  on  tabular  data.  Arrays  Y and  X 
store  the  dependent  and  independent  tabular  values.  The  values  must  be  monotonieally 
increasing,  but  need  not  be  equally  spaced.  There  are  Ml  table  entries,  and  Yl  = Y(X9) 
is  the  value  returned.  The  dimension  statement  (line  30)  assumes  MK50.  Extrapolation 
is  performed  if  required. 

7000-7170.  Evaluates  fflV,  X) 

I his  routine  evaluates  the  extinction  coefficient,  B^f  (in  km-')  for  a given  value  of 
wavelength  in  microns,  L5,  and  a given  value  of  meteorological  range  in  km,  V5.  The 
method  used  is  described  in  appendix  A. 


Caution 


There  is  a potential  numerical  difficulty  with  this  program  as  it  appears  in  the 
listing.  Care  should  be  taken  when  running  this  program  in  single  precision  arithmetic. 
Care  was  taken  to  assure  that  this  numerical  problem  would  not  and  did  not  affect  any 
of  the  numerical  values  presented  in  this  report. 
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0020 
0030 
0035 
0037 
0040 
0060 
0070 
0080 
0090 
0092 
0094 
009  5 
009  6 
0097 
0100 
0110 
0120 
0130 
01  40 

01  50 
0160 
0170 
0160 
0210 
0220 
0230 

02  40 
02  50 
0260 
0270 
0280 
0290 
0300 
0370 
0380 
0390 
0400 
0410 
0420 
0430 
0440 
0450 
0460 
0470 
0480 
0490 
0495 
0500 
0510 
0520 


REM  OTH  PROPAGATION  MODEL 1 NG  PROGRAM. 

RtM  MEGATEK  CORPORATION.  DECEMBER.  197  5 
DIM  WC21  1.BC21 l.XC  501. YC  501.PU41.Mt  I 31.  LI  I0J.NI  I3J 
DIM  0( 10}. VC  37 

DtP  FNA( X)*10*.43  429  S*L  OG ( X ) 

L tT  R 1 *4*  6370/3 
RtM  FILL  PARAMtTER  ARRAY. 

GOSUB  4130 

INPUT  “AEROSOL  DISTRIBUTION!  MARINE! I > OR  CONTI NtNTAL ( 2) I ". fa  1 
PRINT  “ " 

INPUT  “FAN!  I)  OR  PENC1L(2)  bEAMI".fa4 
IF  04-2  GOTO  0097 

PRINT  "DEFAULT  FAN  ANGLE*"! PC  61! " DEGREES . " 

GOTO  0110 

PRINT  "DEFAULT  VERTICAL  FIELD  OF  VI EW»" 1 PC  7 1 1 " 

REM  ChECK  TO  SEE  IF  ANY  PARAMETER  CHANGES 
INPUT  "PARAMETER  CHANGES  (l»YES.  0*NO)l".fa2 
PRINT  " " 

IF  02*0  GOTO  0210 
GOSUB  4570 

REM  DETERMINE  INDEPENDENT  VARIABLE.  VAhI AbLES 

REM  1 WAVELENGTH  (MICRONS! 

REM  2 PATH  LENGTH  (KM) 

REM  3 METEROLOGICAL  RANGE  (KM! 

PRINT  "INPUT  VARIABLES.  ENTER  START  VALUE  FOR 
PRINT  ”1.  WAVELENGTH  I "i 
INPUT  VI 13 
PRINT  " " 

PRINT  "2.  PATH  LENGTH!"! 

INPUT  VC  2 I 
PRINT  " " 

PRINT  "3.  METEROLOGICAL  RANGE!"! 

INPUT  VC  31 
PRINT  “ " 

PRINT  "INPUT  SUBSCRIPT  OF  INDEPENDENT  VARIABLE! "1 
INPUT  fa 6 
PRINT  " " 

LET  VI  « VC  fab  J 

PRINT  "INCREMENT.  NUMBER  OF  VALUES!"! 

INPUT  V2.V3 
PRINT  " " 

REM  DAY  OR  NIGHT  AFFECTS  bACKGROUND  POWER  VALUE. 

INPUT  "DAK  I!  OR  NI  GHT  ( 2 ) I "«  fa3 

rem*****LOOP  OVER  INDEPENDENT  VARIABLE  STARTS  HERE***** 

FOR  fa  5* I TO  V3 

REM  CALCULATE  P5«PR/( PT*A*F ) . WHERE  PR*POWER  RECEIVED.  PT= 

REM  POWER  TRANSMITTED.  A*KEC VR  AREA.  F*SCATTER1NG  F'N. 

LET  Vl«61*(fa5-1 )* V2*V1 

RtM  CALCULATE  EXTINCTION  COEFFICIENT.  IF  V( 1 ) OR  V(2)  HAS  CHANGED. 
IF  65*1  GOTO  0540 
IF  fa6< >3  GOTO  0540 


DEGREES.* 


ARE! 

L5 

Dl 

V5 

INDEPENDENT  VAR" 


P530  GOTO  AS7P 

0540  LET  V5- VI  3 J 

0550  LET  L5-VIIJ 

0560  GOSUB  7850 

0570  LET  D I ■ VI 2 1 

0580  IF  84“  1 GOTO  0620 

0590  REM  PENCIL  BEAM 

0600  GOSUB  2510 

0605  LET  P5-23 

0610  GOTO  0650 

0415  REM  PAN  BEAM 

0620  GOSUB  2 600 

0630  LET  PS-24 

0635  REM*****AT  THIS  POINT  P5«”PATH  LOSS"  AS  DEf 1 NED  ABOVE***** 

0640  REM  COMPUTE  SCATTERING  FUNCTION.  P5 

0 6 50  IP  6 1 *2  GOTO  0750 

0655  REM  MARINE  MA2E 

0660  POR  11-1  TO  13 

0670  LET  Xllll-NUIJ 

0680  LET 

0690  NEXT  11 

0700  LET  Ml-13 

0710  LET  X9-VI 1 1 

0720  GOSUB  6030 

0740  GOTO  0830 

07  50  HEM  CONTINENTAL  HAc E • 

0760  POR  11-1  10  10 

0770  LET  TCI  I 1-LC1 1 1 

0780  LET  All  1 1-01.11  1 

0790  NEXT  11 

0800  LET  HI-10 

0010  LET  X9-V1 I J 

0000  GOSUB  6030 

0830  LET  PS-M 

0840  REM  P 6-KECE1 VED  ROGER 

0850  LET  P6-P5*Pt  13 J*3. 1 41  59*K  121*PC  12 1*P5 

0860  REM  COMPUTE  S/N  RATIO. 

0870  REM  COMPUTE  P7 -BACKGROUND* 

087  I LET  P8-0 

0078  IF  VI 11-1.8  GOTO  0880 

0875  LET  PS-1 .2E*08/( VI 1 1 » S*< EXP< 50/ VI l 1 J - I > > 

0800  IP  W3-2  GOTO  0950 

0890  REM  OAT  TIME  BACKGROUND 

0900  LET  P7-4750/1 VI  I 1 • 6. 5*< EXP ( 2 . 44/ VI 1 J > - 1 > »*P8 

0930  GOTO  0960 

0940  REM  NIGH TIME  BACKGROUND 

0950  LET  P7-.  0083/1  VI 1 3 » 6. 5*< EXPC 2 . 44/ VI 1 3 > - 1 > ) *P8 

0960  LET  P7-P7-6. 28318*1 1-CQSl F2/8J )*3« 14159*PI I21*PI I2J*PI 141 

0970  LET  51-2*1 . 602I9E-I9*PI I0)*(P19 l*(P7*P6»*Pl 1 1 1) 

0971  LET  S1-P191*P191*P 6*P6/S 1 

0980  R EM* •••• AT  THIS  POINT  Sl-SNR***** 


HA 


jJ 


0990  RIM* ••••OUTPUT  ROUTINE  SHOULD  BE  INSERTED  HERE***** 

1000  PRINT  V106J.P5.  FNA(P5)»S1* FNA< SI > 

1010  NEXT  OS 

1020  PRINT  "END  OF  RUN.  ANOTHER  tl-Y«  0-NO)l"J 

1030  INPUT  08 

1040  1 f 08*0  GOTO  9999 

I0S0  GOTO  0080 

2S00  REM  FOV  1NTEG 

2S10  LET  A3-ATNCSUR(2*H2/R1 ) ) 

2S20  LET  A4-C0S(A3)*C0S(P2)*S1N(A3>*< 01 /<2*R1 > ♦ IH2-H I > /01 > 

2 530  LET  A4-A4/S0R( 1 *(01/(2*21 ) ♦IH2-H 1 ) /D1 ) • ( D1 /( 2*R 1 ) ♦CH2-H 1 ) /D1 ) ) 
2540  LET  A4-ATNC  SOK( 1 -A4*A4) /A4> 

2545  LET  23-0 

2550  FOR  B2-A4  TO  A4*F2  STEP  .003 

2560  GOSUB  3000 

2570  LET  23-23*25*. 003 

2580  NEXT  B2 

2590  RETURN 

2600  REM  FAN  AVG 

2610  LET  A5-ATN(S0h<2*Hl/Kl>) 

2615  LET  24-0 

2620  FOR  A6--A5  TO  -A5*FI  STEP  H/5 

2630  LET  AI--A6 

2640  GOSUb  2510 

2650  LET  24-24*23*. 2 

2655  NEXT  A 6 

2660  RETURN 

3AAA  PEM  X5«PR/<PT*A<SQM>*E< A)  ) 

3A1A  LET  PI -COSf A1 >*C0S<P1 >*SIH<AI >*(D1 /CP*P1 >*<MI -HPl/PI > 

TAPP  LET  PI -PI / SOR( I ♦ ( (MP-mi »/pl-ni/(?*R|»)«(fM?-HI)/m-DI/f?*Rl))) 
3A3A  LET  PI -ATNf  SOPM -PI -PI  ) /PI  > 

3A4A  LET  »1«D1*SIN<RP)/STNfP1*PP> 

3A5A  LET  ?P-DI*S!N(P| >/StN(P| .RP) 

TRAP  LET  A?-SIN(A1 >*C0SfD|/R1 )-COS(PI > -COS( A n *S IN ( PI /R I > 

7A7A  LET  AP- (AP*» I .HP-HI .HI *DI /( P-P1 > >/*P 
3 AAA  LET  AP- ATN ( AP/SOR <|*Ap*A?>> 

3A9A  LET  V9-S0P( P I / ( ?*MA ) > -TAN ( A I > 

3A95  LET  YA-M*C0S(A1  1 / SOR < ?*R 1 *H A > 

3A96  I.ET  XI--Y9 
3A97  GOSUP  35 1 A 
7A9A  LET  GP-GI*.Sr,Nf-Y97 
7A99  LET  YI-YA-Y9 
31AA  GOSUB  35 1 A 

3 1 1 A LET  GP-GP-GI*EYP(P*Y«*Y9-YP*YA)«SGN<YP-Y9) 

3 I PA  LET  GP-GP*EXPf Y9* Y9)*(SGN< YA- Y«)-SGN(-Y9» ) 

3I3A  LET  LI»-BA*S0RfP*PI*MA>*G?«.9P6PP7*EXP(-HI /WA»/COS(A|) 

3145  LET  L 1 -EXP<L 1 ) 

315A  LET  Y9-50R»R1 / < P-HA1 1-TAN< APT 
31AA  LET  Y«-X?*rOSCAP)/S0Ff P*PI «HP > 

3 1 7 A LET  XI--Y9 
31  RA  GOSUB  351  A 
7 I A5  LET  GP-G1«SGNf-Y91 

3 1 9A  LET  xi-y*»Y9 
71*5  GOSUP  75 1 A 

7PAA  LET  GP-G?-G1-E»PfP-YP-Y9-YP»YB»*SGN<Y9-Y9) 


85 


< 


3P05  LET  GP-GP*eXP(Y9*Y9)*<SGN< YR-Y9)-SGN<-Y9> ) 

3P06  LET  L2--P0*SOP(PwPl6M0>6GP*. AP6PP7-EXPY-MI /H0>/COS(A?> 

3P07  LET  LP-EXP(LP) 

3P10  LET  P5-M1**1wP|6C0S<AnwC0SCAl )/(2*PI )-PI*SIN( Al » 

3PPP  LET  93-  • 000001 *L I wLPwB®wCXP< -F5/H0) /<  4*3.t4l39*Dt*SIN<BI>) 
3P30  RETURN 
3300  REM  GirXI) 

3910  LET  Tl-l/(l*.  47047  6ABS( X I > > 

3SP0  LET  ai-((<.T4TB56wTI»-9.3«79Ae-0P)*TI*.34PB24>*TI 
3330  RETURN 


4000 

REM 

SUBROUTINE  TO  SET  MODEL  PARAMETERS. 

4010 

REM 

SUB 

WUANITY 

DEFAULT 

NAME 

4020 

REM 

VAULE 

4030 

REM 

I 

TRANSMITTER  HEIGHT  (M) 

25 

HI 

(KM) 

4040 

REM 

2 

RECEIVER  HEIGHT  (M) 

25 

H2 

(KH) 

4050 

REM 

3 

At.  POINTING  ANGLE  (DEG)  KEL.  TO 

0 

PI 

(RAD) 

40  60 

REM 

4 

GREAT  CIRCLE  PATH  ( 1 -TRANS*  2-KECVR) 

0 

P2 

(RAD) 

4070 

REM 

5 

AEROSOL  LAPSE  RATE  (KM) 

100 

H0 

(KM) 

408  0 

REM 

6 

VERTICAL  FAN  ANGLE(DEG)-FAN  CASE  ONLY 

4 

FI 

(RAD) 

4090 

REM 

7 

VERTICAL  FIELD  OF  VIEW  (DEG) 

1 

F2 

(RAD) 

4100 

REM 

8 

ELEVATION  ANGLE  (DEG)  REL.  TO  LOCAL 

0 

Al 

(RAD) 

41  10 

REM 

H OR 1 2 . OF  TRANS.  (POS«  DOWNWARDS) 

4120 

REM 

9 

DETECTOR  KESP0NSIV1TY 

.065 

KV  ) 

4130 

REM 

10 

SYSTEM  ELECTRICAL  BANDWIDTh  (HE) 

7E7 

PC  10) 

41  40 

REM 

1 1 

DETECTOR  DARK  CURRENT  (AMPS) 

IE- 14 

P(  1 1 ) 

41  50 

REM 

12 

RECEIVER  APERTURE  RADIUS  (M) 

.05 

PC  12) 

4160 

REM 

13 

TRANSMITTER  POWER  (W) 

1E6 

PC  13) 

4170 

REM 

14 

RECEIVER  OPTICAL  BANDWIDTH  (MICRONS) 

1 E-3 

P(  14) 

4180  RESTORE 

4190  FOR  19-1  TO  14 

4200  READ  PC  19} 

4210  NEXT  19 

4220  DATA  25*25,0*0*  100*4*  1*0*  .065»7E'»07.  I E- I 4.  .03.  1 E»06*  .001 
4230  REM  STUFF  W AND  B ARRAYS  f OR  EXTINCTION  COEFFICIENT. 

4240  FOR  19-1  TO  21 
42  50  READ  WC191«BC19} 

4260  NEXT  19 

4270  DATA  .27.2* .28* 1 .89, .3, 1 .78* .32* 1 . 67, .34* | . 56 
4280  DATA  . 3 6,  1 . 45,  • 38  * 1 . 4*  . 4*  1 . 3*  • 4 5*  1 . 1 5*  . 5*  1 . 0 5*  . 55*  .9  66 
4290  DATA  . 6* .8 6* . 65* . 78* • 7 * . 73* .8* . 64* . 9* . 58 * I . 06* . 52 
4300  DATA  I .26. .47, 1 .67* .4*2.17* .36* 10, .312 
4310  REM  STUFF  SCATTERING  COEFFICIENTS  - M*  N- 

4320.  REM  VALUE*  WAVELENGTH  FOR  MARINE  HA2E)  L*  0- VALUE*  WAVELENG  TH 
4330  REM  CONTINENTAL  HA2E. 

4340  FOR  19-1  TO  13 
4350  READ  N!19]*MC19] 

4360  NEXT  19 

4370  DATA  .45* 10.31* .7*5. 176* 1 . 19*2.862* 1 .45*2.423* 1 .61*2.224 
4380  DATA  1 .94* I .939*2.25* 1 .776*3* I • 103*3.9* .921 6* 5.3* . 6851 
4390  DATA  6. 05*. 5279*8. 15*. 4154*10*. 3256 
4400  FOR  19-1  TO  10 
4410  READ  01 1 9 1 *LC l 9 } 

4420  NEXT  19 


X(> 


♦ 


4430  DATA  .45* 3.72b* .7*2.41 5* I • 19* I . 52* I .9 4* .9252* J* >441  I 
4440  DATA  3.9* .3458* 5.3* .2891* 6.05* .243*8. 1 5* . 1971* I 6. 6* . 1393 
4445  GOSUB  4470 
4450  RETURN 

4460  REM  E0UATE  VARIABLES  TO  P ARRAY  AND  CONVERT  UNITS. 

4470  LET  H1«PC  11*. 001 

4460  LET  H2-PC2J4.001 

4490  LET  Pl“P£3)*1.74533E-02 

4500  LET  P2*PC4J4l .74533E-02 

4510  LET  H0-PC  51 

4520  LET  FI*P£6)*1.74533E-02 

4530  LET  F2*PC7 1*1 .74533E-02 

4540  LET  A|*P(8)*I .74533E-02 

4550  RETURN 

4560  REN  SUBROUTINE  TO  CHANGE  PARAMETER  VALUES. 

4570  PRINT  "ENTER  SUBSCRIPT  OK  PARAMETER  TO  BE  CHANGED  lEUR  A LIST" 
4580  PRINT  "OF  SUBSCRIPTS  AND  CURRENT  VALUES*  LIST  LINES  4010  " 

4590  PRINT  "TO  4170)1"! 

4600  INPUT  19 
4610  PRINT  " " 

4620  IF  19*0  GOTO  4700 
4630  PRINT  "NEW  VALUE!"! 

4640  INPUT  PC  19) 

4650  PRINT  M " 

4660  PRINT  "NEAT  SUBSCRIPT  <0  1*  DONE)l"! 

4670  INPUT  19 
4680  PRINT  " " 

4690  GOTO  4620 
4700  GOSUB  4470 
4710  RETURN 

6000  REM  LINEAR  INTERPOLATION  ON  A LOG-LOG  SCALE.  M I ■NUMBER  OF 
6010  REM  POINTS  IN  A AND  7*  WHICH  STORE  TABULAh  VALUES.  Yl« 

6020  REM  Y t X9 ) IS  VALUE  RETURNED.  A IS  MONOTON1 CALLY  INCREASING. 

6030  LET  I 7 -M 1 - 1 

6032  IF  X9*XC17)  GOTO  6120 

6035  FOR  17*H1  TO  1 STEP  -1 

6040  IF  X9>-XC17)  GOTO  6090 

6050  NEXT  17 

60  60  LET  15*1 

6070  LET  1 6*2 

6080  GOTO  6140 

609  0 IF  X9 *Xl 17)  GOTO  6120 

6100  LET  Yl*YC17) 

6110  RETURN 
6120  LET  15*17 
6130  LET  I 6* I 7* 1 

6140  LET  M9*( LOG< YCl6J)-LOG(Ytl5))> /(LOGC XI 1 6))-LOG(XIl 5)) ) 

6150  LET  B9*-M9*LQG( XI 15))*L0G(YC 1 5)) 

6160  LET  YI-EXP(M9*LOG(X9)*B9) 

6170  RETURN 


X7 


7 


7*00  REM  EVALUATES  BETA  ■ B0  (IN  KMt-1)  BY  INTERPOLATION  USING 
7010  REM  ELTERMAN’S  MODEL  (APPLIED  OPTICS*  9*  1804-10*  1970). 
70S0  REM  INPUTS  ARE  VS*  THE  METEORLOGI CAL  RANGE  IN  KM*  AND 
7030  REM  LS*  THE  WAVELENGTH  IN  MICRONS.  B IS  ASSUMED  PILLED 
7040  REM  WITH  THE  21  VALUES  OP  BETA  POR  VS-4  KM.  L STORES 
704S  REM  THE  WAVELENGTHS. 

7050  IP  L 5*  • • 2 GOTO  7090 

7070  PRINT  ’’WAVELENGTH  OUT  OP  RANGE  POh  BETA.  L5-’’*LS 
7000  STOP 

7090  IP  L 5> I I GOTO  7070 
7100  POR  17*1  TO  21 
7110  LET  Xt  1 7 J-Wl 17  J 
7120  LET  Yl 17  J-Bl 17  ) 

7130  NEXT  17 
7140  LET  X9*L5 
7145  LET  MI-21 
7150  G0SU8  6030 

7160  LET  B0-Yl*((3.9l/V5)-.01l6)/((3. 91/4)-. 01t6>*.  0116 
7 170  RETURN 
9999  END 


APPENDIX  C:  LOS  PROGRAM  DESCRIPTION 


The  line  of  sight  (LOS)  propagation  model  is  inherent  in  the  OTH  model  described 
in  appendix  E.  The  basic  quantity  calculated  is  the  fraction  of  power  remaining  in  the 
beam  after  propagating  a distance  z measured  along  the  beam.  In  the  notation  ot  appendix 
E.  this  is  K(z,  a),  given  in  equation  7)  of  that  appendix  (where  a is  the  angle  with  respect 
to  the  local  horizontal).  The  LOS  program  was  derived  from  the  OTH  interactive  program 
described  in  appendix  C.  As  a result,  the  program  comments  given  in  appendix  B for  all 
lines  numbered  4230  and  higher  apply  to  the  LOS  program. 

The  subroutine  from  line  3000  to  line  3158  computes  8(z,  a).  Appendix  F.,  (in 
equation  7 and  the  accompanying  footnote)  gives  the  expression  for  8(z,  a)  and  a numerical 
approximation  for  the  error  function  required  to  compute  it.  However,  certain  values  of 
the  parameters  cause  overflow  to  occur  in  evaluating  the  exponentials  if  this  expression  is 
used  directly.  The  code  given  in  the  listing  takes  advantage  of  the  fact  that  the  exp(-x-) 
term  in  erf(x),  for  the  values  of  x involved,  partially  cancels  the  exponential  term  multi- 
plying the  G function  which  causes  the  problem.  This  eliminates  the  overflow.  Lines 
3500-3530  evaluate  Gj(x),  where 

G(x)  = (\/lrl2)[  l+O-G j(x)exp(-x2)  sgn(x)] 

and  the  G function  is  as  defined  in  appendix  E. 

The  rest  of  the  statements  in  the  LOS  program  carry  out  three  loops  over  the  para- 
meters of  interest,  and  print  a table  of  results.  The  columns  contain  the  path  loss,  and 
the  signal-to-noise  ratio  for  the  four  combinations  of  powers  of  1 MW  and  1 kW,  and 
receiver  fields  of  view  of  1°  and  60°.  Each  row  is  a different  wavelength.  Lines  225-351 
assign  values  to  the  other  system  and  geometric  parameters.  Their  meaning  is  documented 
in  the  program  remarks,  and  the  variables  used  (primarily  the  P array)  are  consistent  with 
the  description  of  appendix  B.  The  three  loops  are,  from  outer  to  inner,  over  the  values 
of  meteorological  range,  propagation  distance,  ana  wavelength.  A separate  table,  giving  the 
columnar  values  as  a function  of  wavelength,  is  produced  for  each  value  of  meteorological 
range.  The  values  of  meteorological  range  (in  km)  used  are  inserted  via  the  data  statement 
at  line  1 10.  The  values  of  distance  (in  km)  are  inserted  at  statement  160,  and  the  wave- 
length values  (in  microns)  are  inserted  at  statements  210-221. 

The  only  real  difference,  other  than  looping  and  parameter  value  initialization, 
between  the  computations  in  the  LOS  program  and  the  analogous  ones  in  the  OTH  pro- 
gram is  found  in  lines  435-480.  Two  LOS  cases  were  considered.  In  the  first,  which 
substitutes  the  lines  of  code  given  below  for  those  in  the  listing,  it  was  assumed  that  oper- 
ational advantages  (primarily  covertness)  might  cause  the  RPV  to  be  flown  at  the  minimum 
altitude  necessary  to  propagate  over  the  desired  range.  The  lines  given  below  compute  this 
minimum  height  and  the  associated  elevation  angle  of  the  transmitted  beam.  The  geometry 
is  as  shown  in  figure  Cl . The  parameters  shown  are  as  follows: 

ht  = altitude  of  transmitter 

hr  = altitude  of  receiver  (RPV) 

a = elevation  angle  with  respect  to  horizontal  (+  downwards) 

a = radius  of  earth  (k=4/3  is  refraction  correction) 

z = propagation  distance  along  path 

D - propagation  distance  over  earth 


Figure  Cl.  Geometry  used  to  compute  minimum  RPV  altitude. 


The  following  geometrical  relationships  hold: 
x = j(ka  + ht)2  - (ka)-j 

y = |(ka  + hr)2  - (ka)~J 

z = x+y  = (2ka  ht)'/2  + (2kahr)’/2  (Cl) 

The  minimum  value  of  hr  for  a given  z and  ht  is  thus 

hr  = fz  - (2kaht)‘/2j  ~/2ka  (C2) 

Ihe  approximation  that  z s D is  made.  This  can  he  seen  to  be  good  for  practical  values 
of  D by  noting  that,  using  the  angles  defined  in  the  Figure, 

D = Oka  = ka0  + ka0^  (C3) 

0 1 = arctan  (2kaht)'/2/kaJ  (C4) 

0i  = arctan  f(2kahr)1/2/kal 


(C5) 


If  0 | and  #2  are  small,  the  small  angle  approximations  for  expressions  (C4)  and  (C5)  may 
be  substituted  into  expression  (C3)  to  yield  expression  (Cl),  ie,  D ^ z.  As  an  example: 
hr  = 3km,  0-  a;  0.035  rad,  for  which  the  approximation  is  excellent.  From  geometry,  it 
also  follows  that  a = 6\.  The  lines  which  should  replace  lines  435480  in  the  listing  are 
given  in  figure  C2. 

The  code  in  the  listing,  lines  435-480,  treats  the  second  LOS  case.  Given  an  over- 
the-ground  communication  range,  D,  and  values  for  the  transmitter  and  receiver  heights, 
the  slant  range,  z.  and  elevation  angle  with  respect  to  the  horizontal  at  the  transmitter, 
a,  are  computed.  The  geometry  is  shown  in  figure  C3.  From  the  law  of  cosines,  with 
R = ka. 


z2  = (R  + h,)2  + (R  + hr)2  - 2(R  + h,)  (R  + hr)cos0  (C6) 

Noting  that  0 = D/R.  expanding  sinO.  and  keeping  only  the  most  significant  terms  in  the 
result  yields 


a = cos"1  |(D/z)|  1-D2/(6R2)  + hr/R]j  (C7) 

Note  that  in  the  LOS  program  listing,  the  arctangent  (ATN)  function  must  be  used  because 
BASIC  has  no  arccosine.  This  requires  some  method  of  determining  the  sign  of  a (which 
is  positive  if  the  LOS  is  below  the  horizontal  and  negative  if  above).  To  do  this  the  length 
of  the  horizontal  path  (a  = 0)  between  the  transmitter  and  a point  above  the  receiver  is 
computed  (Z6,  lines  473475)  and  compared  with  z.  If  z is  larger  than  this  value,  a < 0; 
if  z is  shorter,  a > 0.  For  very  small  values  of  a.  or  for  ht  ~ hr,  cancellation  errors  can 
cause  problems  in  this  comparison.  If  this  is  a problem  for  a case  of  interest,  it  is  sug- 
gested that  SQRlZb)  be  expanded  as  a power  series  in  1/R. 

The  user  is  urged  to  heed  the  "Caution”  in  appendix  B.  It  applies  to  the  LOS  pro- 
gram as  well. 
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0 43  5 REM 
043  6 REM 
0 4 40  REM 
0453  LET 
0460  LET 
0 47  0 LET 
0472  LET 
0 47  3 1 C £1**1 
047  5 LET 
047  6 LET 
0430  REM 


H2  = RFX£!VER  HEIGHT 

CALCULATE  NIN  H2  (KM)  FOR  GIVLN  01.  HI. 

APPKOK  I MAT  I ON  J D1  (OVER  GROUND)  5 (SLANT  RANGE) 

c S-Hl  C2  J 

01=e5 

h2=£S-SQR(2*R1 *H1 ) 

H2=H2*H2/<2*R1 ) 

A 1 = EL  E VAT  I ON  ANGLE*  POSITIVE  DO'  'NV.'ARDS*  RAO. 

A 1 = ATN( SQR( 2*K 1 *H 1 ) /HI ) 
n i =3  5 

LOOP  OVER  WAVELENGTH 


Figure  C2.  Code  to  calculate  minimum  receiver  height  necessary 
to  maintain  tine  of  sight . 


Figure  C3.  Geometry  used  to  compute  elevation  angle 
and  slant  range  front  transmitter  to  KI’V 


Program  Listing: 
LOS  Modeling  Program 


10  REM  LINE  OF  SIGHT  PROPAGATION  MOOELING  PROGRAM. 

20  REM  MEGATEK  CORPORATION,  DECEMBER  1975 

30  DIM  W C 2 1 3 , 8 C 2 1 3»XC503,YC 5B3*GC93,HC 133,1(21 3 » P C 151 

A0  DEF  FNA( X ) = 1 0*  . ^3429 5*  LOG  (X) 

50  RESTOR 

55  REM  R 1=4/3  EARTH  RADIUS  (KM) 

60  LET  Rl=4*6370/3 

70  REM  G=MET.  RANGE  LOOP  VALUES  (KM) 

80  FOR  11=1  TO  9 
90  READ  GC  I 1 3 
100  NEXT  II 

110  DATA  4,  8,  10,  15,  20,  30,  50,  70,  90 

120  REM  H=D1 STANCE  LOOP  VALUES  (KM) 

130  FOR  11=1  TO  13 
1 40  READ  HC  I 1 1 
150  NEXT  II 

160  DATA  20,  40,  60,  60,  100,  125,  150,  175,  200,  225,  250,  275,  300 

170  REM  I =WA VELENG TH  LOOP  VALUES  (MICRONS) 

160  FOR  11=1  TO  21 
190  READ  Kill 
200  NEXT  II 

210  DATA  .4, . 5, . 6, *7, .6, .9, 1, 1 . I,  1 .2,  1 *3 

220  DATA  I • 6,  2 • 1 , 2 • 3,  3 • 5,  3 • 6 » 4 . 1 « 4 • 6,  4 . 9 

221  DATA  8,9, 10.2 

225  REM  Pl=A2.  POINTING  ANGLE  (RAD) 

230  LET  p 1 = 0 

240  REM  H1=TRANS.  HEIGHT  (KM) 

245  LET  PC  1 3=7 

250  LET  H I = .00 1 4P  C 1 3 

266  REM  H0-LAPSE  RATE  (KM) 

270  LET  H0= 1 /. 77 

272  REM  P9*RECEI VER  FIELD  OF  VIEW  (RAD) 

273  LET  PC  1 53  = 1 

274  LET  P9=PC 153*1 .74533E-2 

275  REM  D3=TRANSMI TTER  DIVERGENCE  (RADS) 

27  6 LET  D3=  .001 

280  REM  P(9 ) *DETECTOR  RESPONSi VI TY,  P(I0>=SYSTEM  ELECTRICAL 
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285  REM  BANDWIDTH  (Hi),  PC  1 l > = DETECTOR  DARK  CURRENT  (AMPS). 

287  REM  PC  12) =RECEI VER  APERTURE  RADIUS  ( M ) » PC l 3 ) =TRANSMIT TER  POWER 

300  REM  (W).  PC  1 4)=RECEI VER  OPTICAL  BANDWIDTH  (MICRONS) 

310  LET  PC9I«.065 
320  LET  PC  1 0 I = 7 E*7 
330  LET  PC  1 1 1*IE-14 
340  LET  PC  121*. 05642 

350  LET  PC 13J-1E+6 

351  LET  PC  1 4J=.0001 

352  REM  STUFF  B i W ARRAYS. 

353  GQSUB  4230 

370  REM  LOOP  OVER  MET.  RANGE 

380  LET  Cl*  0 
390  LET  C1*C1-»1 

400  LET  V5=GCCl  1 

401  PRINT 

403  PRINT 

495  PRINT  ~V*"I V5I-KM*  HT«"I PC  1 )J "KM." 

406  PRINT  TAB(25)J”F0V=1")  TAB  ( 3 5)  ) M FO  V*  1 "i 

407  PRINT  TABC 45>)”FOV*60")  TABC  55) J " FO V* 60” 

408  PRINT  TAB(  25)  ) "P  = 1 MW")  TAB<  35)1  *'P=  1 KW") 

409  PRINT  TABC4S))"P*1  MW")  TAB(55))"P*1  KW" 

410  PRINT  "RANGE”)  TAB( 7 )) "LAMBDA”)  TABC 1 5) ) "PATH  LOSS") 

411  PRINT  TABC  25) ) "S/N”)  TAB< 35) I ”S/N”)  T ABC  45) ) "S/N”) 

412  PRINT  TABC 55))”S/N”)  TABC 65) ) "HR" 

420  LET  C2*0 

430  LET  C2-C2+I 

435  REM  H2  = RECEI  VER  HEIGHT 

440  LET  H2*2.5 

445  LET  D1=HCC21 

450  LET  09=D1*DI 

455  REM 

460  LET  08 *H 1 +H2 

463  LET  07»CH1-H2)*CH1-H2) 

465  LET  25*  SQR  C G9*G6*U9/R 1 *07- G9* 09/C  1 2*K 1 *K 1 ) ) 

470  LET  Al*CDl/2  5>*( 1-Q9/C  6*R 1 *K l ) +H2/R 1 ) 

471  LET  Al*  ATN  C SGR  C 1 -A  1 *A 1 ) /A  1 ) 

472  LET  A1--A1 

473  LET  c6=CH2*H2-Hl*Mt»2*Rl*CM2-Hl ) ) 

474  IF  26<  0 GOTO  477 

475  IF  ZS<  SGR  (c 6)  GOTO  477 

476  LET  A I * • A 1 

477  LET  il«25 
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480  REM  LOOP  OVER  WAVELENGTH 

498  LET  C3>0 

580  LET  C3-C3*l 

505  LET  C6*l  5 

510  LET  LS-ICC3J 

520  REM  COMPUTE  B0 

530  GOSUB  7050 

540  GOSUB  3010 

541  IF  Ll»*lE-30  GOTO  0549 

542  IP  C3«»l  GOTO  0546 

543  PRINT  HCC2JI  TAB(7)lL5i  TAB< 1 5>* "«-300"j 

544  PRINT  TAB<  64) J 1NTIH2* 1 000* • S> 

545  GOTO  0784 

546  PRINT  TABC7) JL5J  TAB( 1 5> 4 M<~300" 

547  GOTO  0784 

548  REM  P9-RECE1VER  FIELD  OF  VIEW  (RADS) 

549  FOR  C4>l  TO  2 
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550  LET  PC  t S3-1 

551  IP  C4-1  GOTO  8553 

556  LET  PCI  51*48 

553  LET  P9*PC ISJ*l.74533C-08 

554  REM  LOOP  OVER  POWER 

555  FOR  C5*l  TO  8 

554  LET  Pt 131*1 E*06 

557  IP  C5*l  GOTO  8565 

558  LET  PC  131*1000 

568  REM  P6*RECEIVED  P0WER*Ll*PT*A/<Pl*<25«03/8)  »2> 

563  REM  **«**AT  THIS  POINT  LI  IS  PATH  LOSS*  COMPUTE  SNR>S1**«* 

564  REM  P6*RECEI VEO  POWER»LI*PT*A/tPl*(E5*D3/2> »8> 

565  LET  P6-F5*D3*IAAA/? 

578  LET  P6*L1*PI 131*3*14159*PC121*PCI81/(3*l4l 59*P6*P6> 

57  5 LET  C6*C6*I0 

588  REM  P 6* BACKGROUND* 

598  REM  FORCE  PROGRAM  TO  DAYTIME  BACKGROUND* 

595  LET  03*1 

688  LET  P8*8 

618  IF  L5«l*8  GOTO  8658 

688  LET  P8*l.8E*08/CL5t5*<EXP<50/L5>-l>> 

638  IF  03*2  GOTO  8688 

648  REM  DAYTIME  BACKGROUND 

658  LET  P7  * 47  50/<L5*6*  5*<  EXP <8  • 44/LS)  “ 1 1 > *P8 

668  GOTO  8698 

678  REM  N16HTTIME  BACKGROUND* 

688  LET  P7*«8883/(LS*6»  5*<  LXPt  2*  44/L5) “11 >*P8 

698  LET  P7*P7*<3.l4159/41*P9*P9*3.l4l59*Ptl21*Pll81*PC141 

788  LET  S1*2*1.60219E-19*P£ 101*<Pt91*<P7*P61*PClll> 

718  LET  S1*P191*P191*P 6*P6/S1 

788  REM  *«***AT  THIS  POINT  S1*SNR***** 

738  REM  *****OUTPUT  ROUTINE  SHOULD  BE  INSERTED  HERE***** 

748  IF  SMIE-30  GOTO  8753 

748  IF  C6*55  GOTO  8749 

743  IF  C6«»85  GOTO  8746 

744  IF  C3<»1  GOTO  8751 

745  PRINT  HtC8II  TAB<7ML5l  TAB< 1 5) J FNA'Ll it 

746  PRINT  TAB<C6>* FNAlSllJ 

747  GOTO  8788 

749  GOTO  8791 

751  PRINT  TAB(7) JLSJ  TAB< I 5)1 FNA1L1 >1 

758  GOTO  8746 

753  IF  C6*55  GOTO  8811 

754  IF  C6«»85  GOTO  8757 

755  IF  C3«>1  GOTO  8761 

756  PRINT  HCC81I  TABt7)lL5!  TAB ( 1 5) I FNA< L 1 > I 

757  PRINT  TAB<C6>1"«-300"1 

758  GOTO  8788 

7 59  GOTO  8811 

761  PRINT  TABC71IL5*  TAB! 15>1 FNA(LI) J 

768  GOTO  07  57 
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788  NEXT  CS 
782  NEXT  C4 

704  ir  C3-8I  80T0  MM 

798  GOTO  8588 

791  IF  C3«»!  GOTO  8794 

T9t  PRINT  TAB(C6) J FNA<  SI ) I TAB< 64)1 1NT(H2*1080«. 5) 

793  GOTO  8788 

794  PRINT  TAB(C4>«FNA<SI> 

795  GOTO  8788 

888  IF  C8«13  GOTO  8888 

818  GOTO  8438 

811  IF  C3«»!  GOTO  8814 

818  PRINT  TAB<C6>! •*«“3GG**J  TAB<64)JINTCH2*1000^.5> 

813  GOTO  8788 

814  PRINT  TAB<C6>J"«-300" 

815  GOTO  8788 

888  IF  Cl -9  GOTO  0840 

830  GOTO  8398 

848  PRINT 

858  PRINT 

848  PRINT 

878  GOTO  9999 

880  PRINT  H£C2J1  TAB  <7 it  FNA< LI) 

898  GOTO  784 

3088  REM  COMPUTE  LI -PATH  LOSS  FOR  LINE  OF  SIGHT* 

3810  KEM 

3090  LET  Y9-  SOR  <RI/<2*H0>>*  TAN  <A1> 

3100  LET  Y8*21»  COS  <A1>/  SOR  <2*Rl4H0> 

3110  LET  X I «“Y9 

3120  SOSUB  35*0 

3138  LET  G2*G 1 • SGN  <-Y9> 

3140  LET  X I *Y8“Y9 
3150  GOSUB  3ST0 

3160  LET  G2-G2-GI*  EXP  t 2*Ytt*Y9-Y8*Y8 > * SGn  (Y8-Y9) 

3165  LET  G2-G2*  EXP  <Y9*Y9>4<  SGN  <Y8-Y9>-  SGN  C-Y9>> 

3170  LET  L 1 *-B0*  SOR  < 2*K1 *M0) *G2* *88 6227*  EXP  (-H1/H0)/  COS  <A1> 
3175  LET  LI*  EXP  <L!> 

3180  RETURN 

3500  REM  3 1 (X 1 ) 

3510  LET  Tl*l/tl*. 47047*  ABS  <X1>> 

3520  LET  G 1 *< ( ( *7478 56*T1) “9  * 58  7 98E-2)*T1 ♦•348  0 2 4) *Tl 
3530  RETURN 

4230  REM  STUFF  W A NO  6 ARRAYS  FOR  EXTINCTION  COEFFICIENT* 

4248  FOR  19*1  TO  21 
4250  READ  W£I9J»B£19J 
4260  NEXT  19 


4270  DATA  .27.  2,  ,22,  1.89*  ,3,  1.78*  .32*  1.67*  .34*  1.56 
4288  DATA  .36*  1.45.  .38.  1.4.  .4.  1.3.  .45.  1.15.  .5*  1.85.  .55 
4285  DATA  .966 

4298  DATA  .6.  .86.  .65.  .78.  .7.  .73.  .6.  .64.  .9.  .58.  1.86.  .52 
4388  DATA  1.26*  .47.  1.67.  .4.  2.17.  .36.  10.  .312 
4318  RETURN 

6888  REM  LINEAR  INTERPOLATION  ON  A LOG-LOG  SCALE.  N1 -NUMBER  OF 

6818  REM  POINTS  IN  X AND  Y»  WHICH  STORE  TABULAR  VALUES.  Yl- 

6828  REM  Y<X9>  IS  VALUE  RETURNED.  X IS  MONOTON1CALLY  INCREASING 

• 

6838  LET  I7-M1-1 

6832  IF  X0*XI I 71  GOTO  6128 

6835  FOR  17-MI  TO  1 STEP  -I 

6848  IF  X9—X117)  GOTO  6090 

6858  NEXT  17 

6848  LET  15-1 

6878  LET  16-2 

6880  GOTO  6148 

6898  IF  X9-XC 17 1 GOTO  6120 

6108  LET  Y1-YI171 

6110  RETURN 

6120  LET  15-17 

6138  LET  16-17*1 

6140  LET  H9-(  LOG  (YCI6J)-  LOG  <YCIS])>/(  LOG  (X(16»-  LOG  (X(ISD) 
6150  LET  B9--M9*  LOG  (XI15l>*  LOG  (YCI51) 

6168  LET  Yl-  EXP  <M9*  LOG  (X9)+B9> 

6178  RETURN 

7000  REM  EVALUATES  BETA  - B0  (IN  KM»-1)  BY  INTERPOLATION  USING 

7010  REM  ELTERMAN’S  MODEL  (APPLIED  OPTICS.  9.  1804-10.  1970). 

7020  REM  INPUTS  ARE  VS.  THE  NETEORLOG1CAL  RANGE  IN  KM.  AND 
7030  REM  L5«  THE  WAVELENGTH  IN  MICRONS.  B IS  ASSUMES  FILLED 

7040  REM  WITH  THE  21  VALUES  OF  BETA  FOR  V5-4  KM.  L STORES 

7845  REM  THE  WAVELENGTHS. 

70  50  IF  LS--.2  GOTO  7090 

7878  PRINT  "WAVELENGTH  OUT  OF  RANGE  FOR  BETA.  L5-".L5 
7088  STOP 

7890  IF  L5»U  GOTO  7070 

7180  FOR  17-1  TO  21 

7110  LET  Xt 17 1-W{ 17 ) 

7128  LET  YC 1 7 1-BCI 7 1 

7130  NEXT  17 
7140  LET  X9-L5 
7145  LET  Ml-21 
7150  GOSUB  6838 

7168  LET  88-YI*((3.91/V5>-.8116)/((3. 91/4)-. 0116)*. 8116 
7170  RETURN 
9999  END 
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Appendix  E:  Propagation  Model  for  a Laser-Type 
Beyond-the-Horizon  Communication  Link 

(This  appendix  was  issued  as  a Megatek  Corporation 
Informal  Report.  It  is  included  here  for 
completeness  and  ease  of  reference  in  the  same 
form  as  originally  published) 
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PROPAGATION  MODEL  FOR  A 
LASER-TYPE  BEYOND-THE-HORIZON 
COMMUNICATION  LINK 


GENERAL 

King  and  Kainer1  have  discussed  beyond-the-horizon  propagation  of  laser  beams 
by  means  of  scattering  from  clouds  and  hazes  and  have  specifically  considered  the  use  of 
this  propagation  mode  for  communications.  In  the  present  study  a more  detailed  model 
is  given  for  this  forward  scatter  optical  communications  channel  in  the  cloud-free  case,  ie 
where  scattering  arises  solely  from  atmospheric  “haze”  aerosols.  The  model  takes  explicit 
account  of:  a)  the  vertical  falloff  in  aerosol  concentration;  b)  refraction  and  c)  Earth’s 
curvature,  neglect  of  which  can  lead  to  severe  errors  in  path  loss  estimates  at  ranges  of 
interest  in  beyond-the-horizon  communications.  Indeed,  it  is  found  that  proper  exploita- 
tion of  the  a)  above  can  lead  to  novel  strategies  for  reducing  path  loss,  eg  optimally 
elevated  beams  or  vertical  fans. 


DETAILED  MODEL 

VERTICALLY  INHOMOGENEOUS  ATMOSPHERE 

The  atmosphere  is  taken  to  be  horizontally  homogeneous  but  both  the  local 
refractivity  N,  and  extinction  coefficient  (3  are  assumed  to  decrease  exponentially  with 
height  h: 

1)  N = NQ  exp  [-h/HN] 

2)  0 = 0Q  exp  [-h/Hej 

where  the  zero  subscript  refers  to  the  conditions  at  the  Earth’s  surface.  In  the  CCIR 
standard  atmosphere  NQ  = 289  and  Hj^  = 7.34  km.  We  make  the  simplifying  though  rea- 
sonable assumption  that  extinction  arises  solely  from  aerosol  scattering  with  unity  albedo, 
thus  ignoring  both  molecular  Rayleigh  scattering  and  any  absorption  in  the  aerosol  Mie 
scattering.  Then  for  the  wavelength  range  of  .5  to  10  microns,  0O  will  lie  typically  in  the 
range. 

3)  I0“2  km"1  < 0Q  < 2 X 10"1  km-1 

depending  on  visibility,  although  still  higher  values  of  0Q  are  obviously  possible  in  the  case 
of  ground  tog.  in  equation  2)  is  identified  with  the  lapse  rate  in  aerosol  concentration 
a quantity  which  is  not  well  characterized  statistically  and  whose  magnitude  is  typically  of 
the  order  of  I km. 


103 


■ : W-  - 


REFRACTION 


As  a result  of  the  vertical  decrease  in  refractivity,  light  rays  in  the  atmosphere  are 
curved  (concave  downwards)  rather  than  straight.  Thus  a laser  beam  will  gain  in  altitude 
less  rapidly  with  increasing  range  than  would  be  the  case  in  vacuum,  and  this  effect  may 
be  accommodated  to  an  acceptable  degree  of  accuracy  by  introducing  the  concept  of  a 
fictitious  Earth  radius  R'  related  to  the  true  radius  R by: 

4)  R'  = kR 

where  the  factor  k is  determined  by  the  surface  refractivity,  NQ,  as  shown  in  Table  1. 


TABLE  1.  EARTH’S  RADIUS  CORRECTION  COEFFICIENT  (k) 
FOR  VARIOUS  SURFACE  REFRACTIVITIES  (N^ 


Hi 

k 

250 

1.232 

289 

1.305 

301 

1.333 

350 

1.489 

400 

1.767 

Following  standard  practice,  the  value  k = 4/3  will  be  adopted  in  the  present  study. 
In  this  so-called  “4/3  Earth”  approximation,  it  suffices  to  ignore  refraction  throughout  the 
model  calculations  (ie  assume  straight  line  light  paths)  and  then  to  remedy  this  neglect  by 
the  simple  expedient  of  using  an  Earth  radius  = 8393  km  in  any  numerical  evaluations. 


EARTH’S  CURVATURE  EFFECTS 

Due  to  the  curvature  of  the  Earth  the  altitude  of  a laser  beam  will  not  vary  linearly 
with  range.  Specifically,  if  a laser  beam  is  projected  from  a height  H in  a direction  making 
an  angle  “ below  the  local  horizontal  as  shown  in  figure  1 , then  a simple  geometrical  con- 
struction yields  the  relationship 

5)  (R  + h)2  = (R  + H - Z sin  <x)2  + Z '?  cos2  <* 

where  h is  the  altitude  corresponding  to  a distance  Z along  the  beam.  Solving  for  h,  one 
finds 

6)  h(Z,  oc)  = H-Z  sin  « + Z2  cos2  <*/2R  + 0(l/R2) 


$ 
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BEAM  ATTENUATION 


Consider  a laser  beam  (non-divergent)  projected  is  shown  in  figure  1,  and  let 
B(z,  a)  denote  the  fraction  of  power  remaining  in  the  beam  after  propagating  a distance  z 
measured  along  the  beam.  Then  from  2),  there  follows 


7) 


B(z,  oc)  = exp 


jj  0O  exP  [~h(z,  <x)/He  dzj 

. 

which  after  some  algebra,  may  be  written 

8(z,  oc)  = exp  j-0Q  sec  °c  v/2RHe  exp(-(H  - ~ tan2  a)/He) 

where* 

G(x)  =y^‘JJ  e~U~  du  Crf(x) 


BEYOND-THE-HORIZON  PROPAGATION  MODEL 

As  shown  in  figure  2,  in  the  absence  of  clouds  we  may  distinguish  two  primary 
mechanisms  for  beyond-the-horizon  propagation:  viz.  single  and  multiple  scattering  from 
atmospheric  aerosols.  Since  aerosol  scattering  tends  to  be  peaxed  in  the  forward  direction 
and  to  be  virtually  free  of  absorption  at  the  optical  and  infrared  wavelengths  of  interest, 
it  is  meaningful  to  think  in  terms  of  an  “aura”  of  multiply  (forward)  scattered  photons 
surrounding  the  laser  beam.  As  discussed  in  the  Appendix,  this  aura  will  in  general  atten- 
uate less  rapidly  with  distance  than  the  direct  beam  and,  were  it  not  for  the  effect  of  the 
nearby  Earth’s  surface,  would  be  the  dominant  source /Of  radiance  at  long  ranges.  However, 
the  effect  of  the  Earth’s  surface  is  to  remove  energy  from  the  aura  by  a combination  of 
absorption  and  diffuse  scattering  so  that,  as  shown  in  the  Appendix,  for  near-surface  trans- 
mitter/receiver links  the  aura  will  attenuate  with  range  (D)  approximately  as  e^o^/D.  a 
rate  of  attenuation  more  rapid  (as  will  be  seen)  than  that  of  the  singly  scattered  compo- 
nent. For  this  reason,  the  beyond-the-horizon  propagation  model  is  limited  to  the  singly 
scattered  component.  It  should  be  borne  in  mind,  however,  that  at  ranges  less  than  about 
100  km,  the  aura  can  under  favorable  conditions  (see  Appendix)  - be  the  dominant 
component. 


*A  useful  numerical  approximation  for  the  error  function  is 

2 

erf(x)  = 1 - (a 1 1 + a2t2  + a3t3)  e_x"  + E(x)  x > 0 

where 


I P = .47047 


_ ■ - a,  = .3480242 

I + Px  a2=  -.0958798  a3  = .7478556 

and  | E(x)  | < 2.5  X lO-5  . 
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Consider  the  geometry  shown  in  Figure  3 and  let  the  receiver  aperture  (area  = A)  be 
oriented  normal  to  the  line  of  sight  to  the  scattering  element.  Then  if  PT  is  the  trans- 
mitted power  and  dPR  the  power  received  which  enters  the  receiver  optics  at  elevation 
angles  between  02  and  02  + d02.  there  follows 


10) 


dPr  f(01  + 02) 


Pt*A 


£(zl’al> 


dz 


1 


where  (-)  a | is  the  elevation  angle  of  zj,  « is  defined  in  equation  8),  and  f(0)  is  the  gain 
over  isotropic  for  aerosol  scattering  at  the  scattering  angle  0,  ie  the  ratio  of  the  actual 

scattering  cross  section  at  angle  0 to  that  of  an  isotropic  scatterer  of  identical  total  scat- 
tering cross  section. 

From  the  geometrical  relations 
U)  zj  = D'  sin  02/sin  (0j  + 02) 

dzj 

d02  = ~ sinOi  +/32) 

and  the  result  obtainable  from  8) 

-36 


z2  - D'  sin  /)j/sin03j  + /32) 


A ^ 

Z"  COS*"  a 

z sin  a 1] 

1 lHe  2HeR  " 

He  JJ 

there  follows  from  10) 


13) 


d_Pr  0O  d02f(01  +02)  , 

' 4*D'  sin  0,  £(Z1’°V  «(z2*a2)exP  {* 


H|  + Z j ~ COS'"  a:  | 

ST  2HcR 


Zj  sin  “j 

Hi 


where  Z|  and  z2  are  to  be  evaluated  from  1 1). 

It  now  only  remains  to  relate  the  angles  0,  and  02  to  the  elevation  and  azimuth 
angles  at  the  transmitter  and  receiver  respectively.  To  do  this,  we  set  up  the  geometry 
shown  in  figure  4.  The  transmitter  and  receiver  are  taken  to  lie  in  the  y-z  plane  and  the 
laser  beam  direction  is  specified  by  the  angles  cc,  (measuring  the  negative  elevation  angle) 
and  <Pi  measuring  the  azimuth  relative  to  the  great  circle  path  between  transmitter  and 
receiver  <*2  and  <t>2  are  analogously  defined  for  the  line  of  sight  between  the  receiver  and 
the  scattering  element.  We  find  for  (J,  the  relation: 

14)  cos  0|  = cos  al  cos  *\  s‘h  A (R-*-H2)  - sin  [(R-i-H2)  cos  SI  - R-H,] 

\/  (R+H2)“  sin2  SI  +((R+H2)  cos  SI  - R-H,)2 
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which  can  be  approximated  by 


i 


! 


D 


15)  cos  * cos  cos  0|  + sin  <*j  (—  + (H|-H2)/D) 


./ 1 , (H2-Hl  - D )2 
D 2R 


For  <*2>  we  find 

16)  -sin  = |z(  cos  cos  sinfi  - (R+H2)  + (R+Hj)cosf2  - Zj  sin«j  cosfij  /z2 
which  can  be  approximated  by 


17)  -sin  a2  ~ ^cos0|  cos  “j  -^*-sin«jj 


HrH2 


D2/2R 


In  14-17),  D is  the  ground  range  from  transmitter  to  receiver  and 
18)  SI  = D/r 

Thus  the  basic  propagation  model  is  defined  by  equations  13),  1 1),  14)  and  16),  where 
f(/?l+/)2)  is  to  be  determined  from  a given  aerosol  size  distribution  model.  Since  13)  is  in 
differential  form,  it  must  be  integrated  over  the  field  of  view  of  the  receiver  optics.  To 
do  this,  we  assume  the  field  of  view  is  positioned  as  shown  in  figure  5.*  Then,  if  the 
field  of  view  is  F2  radians,  we  must  integrate  13)  over  the  range 

'9>  0min<02<0min  + F2 
where  (cf.  15)): 

20)  cos  0min  = cos  «2  cos$2  + sin  ocj  + (H2_H1)/D) 


J x , (HI~H2  - D/2R)2 


•Note  that  the  laser  beam  in  the  present  model  appears  to  the  receiver  as  a line  source. 
In  practical  cases,  the  received  power  will  increase  linearly  with  field  of  view  for  fields 
of  view  out  to  about  5 degrees. 
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Figure  5.  Field  of  View  Placement. 


A generalization  is  now  added  to  the  model:  thus  far  we  have  considered  the 
laser  beam  to  be  non-divergent ; we  now  allow  the  beam  to  be  spread  in  the  vertical  plane 
by  F ] radians,  ie  we  consider  the  beam  to  be  a vertical  fan  rather  than  a line.  The  rea- 
son for  this  generalization  is  that  it  turns  out  from  13)  that  the  received  power  will  be  a 
maximum  if  the  laser  beam  direction  is  elevated  above  the  horizon  by  an  angle  which 
depends  on  the  range  D and  the  vertical  lapse  rate  in  aerosol  concentration,  He.  There 
may  thus  be  a practical  advantage  in  spreading  the  beam  in  elevation  and  to  accommodate 
this  possibility  we  introduce  the  vertical  fan  angle  Fj  and  average  13)  over  the  range  of 
beam  elevation  angles: 

2»  “min^l  <ocmin  + Fl 
where 

22)  .an  1 

Finally,  with  the  understanding  that  both  the  field  of  view  (F2)  and  fan  angle 
(Fj)  are  less  than  about  five  degrees,  we  can  approximate  the  scattering  function  f(0]+02) 
in  13)  by  its  value  in  the  forward  direction,  f(0),  thereby  removing  it  from  the  field  of 
view  and  fan  integrations.  Then,  we  may  define  an  overall  transmission  T for  the  beyond- 
the-horizon  link  as: 

23)  PR  = Pt-T 
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NUMERICAL  RESULTS 
EFFECT  OF  VERTICAL  BEAM  WIDTH 

As  a first  application  of  the  model,  we  examine  the  dependence  of  the  transmission 
function  on  the  vertical  beam  width,  ie  the  fan  angle  F|.  To  this  end  we  assume  the  fol- 
lowing typical  parameter  values: 

Receiver  and  Transmitter  height  (Hj,  H2)  = 25  meters 

Azimuth  angles  (0| , = 0 (ie  no  offset  from  great  circle) 

Sea  level  extinction  coefficient  (0O)  = 0.1  km-* 

Receiver  field  of  view  (Fj)  = 1" 

Vertical  lapse  rate  in  aerosol  concentration  (He)  = 1.3  km 

Effective  Earth’s  radius  (R)  = 8393  km  (i.e.  “4/3  Earth”  model) 

and  compute  the  quantity  T/(A*f(0)),  where  A is  the  receiver  aperture  in  square  meters 
and  f(0)  is  the  forward  direction  “gain  over  isotropic”  for  the  aerosol  Mie  scattering,*  as 
a function  of  ground  range  D and  fan  angle  F j . The  results  are  shown  in  figure  6. 

It  is  seen  that  for  ranges  below  1 00  km  the  transmission  is  highest  for  a non- 
divergent  beam,  whereas  at  greater  ranges  transmission  improves  with  increasing  fan  angle. 
(The  effect  saturates  near  4 ; calculations  for  8°  fan  angles  are  very  close  to  those  for  4°). 
Since  from  a practical  standpoint  tracking  and  pointing  requirements  in  the  elevation  plane 
are  far  less  severe  tor  a vertical  fan  than  for  a non-divergent  beam  pointed  at  the  horizon 
it  is  clear  that  a (nominally  4°)  vertical  fan  would  be  the  preferred  mode  in  a communica- 
tion system  - particularly  since  figure  6 indicates  that  the  penalty  paid  at  short  ranges  (ie 
<100  km)  is  less  than  about  3 dB. 


DEPENDENCE  ON  EXTINCTION  COEFFICIENT 

In  figures  7 and  8,  the  quantity  T/A-f(0)  is  plotted  as  a function  of  range  and  sea 
level  extinction  coefficient  for  both  a non-divergent  beam  and  a 4°  vertical  fan.  Parameter 
values  used  are  identical  to  those  listed  above.  Comparing  the  two  figures  it  is  seen  that 
the  vertical  fan  is  particularly  advantageous  at  low  visibilities  (high  0 ).  For  example,  at 
130  km  range  and  0O  = .2  the  gain  is  about  26  dB. 

•Typically  f(0)  can  be  as  large  as  10  while  A will  be  of  the  order  of  .1  meters’  so  that  f(0VA 
is  of  order  unity 
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Figure  6.  Effect  of  vertical  beam  width 
on  Path  Loss  Extinction  Coef- 
ficient = O.lKM'l 


Figure  7.  Laser  Beam  Path  Loss  in  Be- 
yond-the-HoriZon  Propagation: 
Single  Scattering  Model  ( 0 
Vertical  Beam  Width) 


Figure  8.  Laser  Beam  Path  Loss  in  Be- 
yond-the-Horizon  Propagation: 
Single  Scattering  Model  ( 4 
Vertical  Beam  Width) 


Figure  9.  Comparison  of  Path  Loss  Pre- 
dictions with  NRL  Observation. 


The  monotonic  increase  of  transmission  with  decreasing  extinction  coefficient 
(ie  increasing  visibility)  is  noteworthy  since  at  some  point  th*  trend  must  reverse  (T  -*■  0 
as  f}Q  -*  0).  If  -130  dB  is  taken  as  a typical  system  threshold,  then  it  is  seen  that  for  the 
range  of  extinction  coefficients  considered  (.05  to  .25  km-1),  the  maximum  range  for 
communication  varies  from  about  140  to  40  km  respectively. 


COMPARISON  WITH  FIELD  OBSERVATIONS 

Transmission  data  for  beyond-the-horizon  laser  propagation  is  not  apparently 
available.  A limited  test  of  the  model  for  the  4°  vertical  fan  was  possible  once  it  was  real- 
ized that  an  isotropic  (ie  2 w)  source  should  behave  (in  the  elevation  plane)  approximately 
as  a small  angle  vertical  fan  since  aerosol  scattering  is  primarily  in  the  near  forward  direc- 
tion and  consequently  photons  emitted  appreciably  away  from  the  horizontal  will  not  be 
scattered  to  beyond  the  horizon  ranges. 

In  1964,  Curcio  and  Drummeter2  at  the  Naval  Research  Laboratory  reported  the 
results  of  an  experiment  in  which  an  uncollimated  6000  joule  xenon  flash  lamp  was 
mounted  about  3 meters  above  water  level  on  the  cabin  of  a small  boat.  During  night 
time  measurements,  the  boat  traveled  a course  on  Chesapeake  Bay  up  to  a maximum 
range  to  the  receiver  of  41  km.  The  receiving  system  had  a 3°  field  of  view  and  was  sim- 
ilarly located  at  3 meters  height.  The  meteorological  range,  monitored  by  a transmissom- 
eter,  was  13  km. 

In  figure  9,  the  Naval  Research  Laboratory  data  of  visible  light  intensity  (in  the 
4700-7000  A band)  vs.  range  is  compared  with  the  predictions  of  the  4°  fan  model  for 
Hj  = H2  = 3 meters.  The  extinction  coefficient  was  chosen  to  be  0O  = .2575  which  is 
what  one  estimates  at  6000  A for  a 13  km  visibility  using  scaling  laws  discussed  elsewhere3. 
Since  only  relative  transmission  data  were  given  in  the  Naval  Research  Laboratory  report, 
a comparison  of  transmission  vs.  range  was  possible  by  adjusting  the  theoretical  curve  to 
the  experimental  data  at  one  range  chosen  arbitrarily  at  20  km  - and  observing  the  fit  at 
other  ranges.  As  seen  in  figure  9,  the  agreement  is  reasonably  good.  Clearly,  however, 
this  limited  data  cannot  be  construed  as  a test  of  the  model  at  the  ranges  of  practical 
interest,  ie  of  the  order  of  100  km,  and  further  experiments  would  be  needed  for  this 
purpose. 


DISCUSSION 

An  analytical  model  for  beyond-the-horizon  laser  propagation  has  been  developed 
based  upon  the  single  scattering  approximation.  Numerical  results  suitable  for  estimation 
of  received  power  levels  in  specific  systems  have  been  given.  The  model  appears  to  be  in 
reasonable  agreement  with  limited  field  data  on  beyond-the-horizon  optical  propagation 
at  (relatively  close)  ranges  up  to  40  km.  A more  conclusive  field  test  of  the  model,  par- 
ticularly at  longer  ranges,  would  be  necessary  to  assess  its  utility  for  beyond-the-horizon 
optical  communications  system  design. 

A clear  conclusion  ot  the  study  is  the  desirability  of  providing  a vertical  beam 
divergence  of  the  order  of  4e  as  a means  of  increasing  the  received  signal  at  distances 
greater  than  100  km  and  to  relax  beam  pointing  and  tracking  requirements  in  the  elevation 
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plane  This  is  a particular  consequence  of  the  decrease  of  aerosol  concentration  with 
height.  It  is  recognized,  however,  that  a vertical  fan  will  reduce  the  bandwidth  of  the 
link  due  to  the  increased  multipath  delay. 

Implicit  in  the  model  is  the  neglect  of  multiple  scattering,  a semi-quantitative 
justification  for  which  is  given  in  the  Appendix  for  the  limit  of  long  ranges.  The  trans- 
mission estimates  of  the  present  model  should  therefore  be  veiwed  as  conservative,  per- 
haps by  as  much  as  10  dB  at  ranges  below  100  km  where  the  aura  can  under  favorable 
conditions  be  dominant. 
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APPENDIX:  MULTIPLE  SCATTERING 


Consider  a perfectly  collimated  beam  propagating  in  a homogeneous  medium  far 
from  any  boundaries.  Then  from  multiple  scattering  theory4  it  is  known  that  if  the  scat- 
tering law  emphasizes  forward  scattering  so  that  the  r.m.s.  scattering  angle  9 < 1,  then 
an  “aura”  of  multiply  scattered  photons  will  develop  (similar  to  that  shown  in  figure  2) 
whose  radial  extent,  r,  increases  as  the  3/2  power  of  the  distance  z along  the  beam.  More 
specifically,  if  the  mean  free  path  for  scattering  is  8,  then 

A-l ) r = C z3'2  fit'* 


where  C is  a constant  of  order  unity. 

Unfortunately,  this  result  cannot  be  directly  applied  to  the  problem  at  hand  be- 
cause of  the  proximity  of  the  Earth’s  surface  to  the  beam.  A downward  scattering,  even 
at  an  angle  small  compared  to  9,  will  result  ir.  the  photon  striking  the  surface  at  which 
point  it  is  either  diffusely  reflected  into  non-forward  directions  or  absorbed  but  in  any 
event  lost  to  the  aura.  Only  those  forward  scatterings  at  sufficiently  small  angles  to  avoid 
hitting  the  surface  will  contribute  to  the  aura. 

We  can  accommodate  this  situation  by  the  following  approximate  device.  Let  us 
assume  the  aura  to  be  built  up  from  multiple  scattering  events  for  which  the  individual 
scattering  angles  are  all  less  than  or  equal  to  some  critical  angle  0Q.  That  is.  if  the  actual 
differential  scattering  cross  section  is  oa(0),  then  the  effective  cross  section  oe(0)  is 


A-2)  ae(0) 


o.,(0)  0<0C 

0 0>OC 


Assuming  lor  simplicity  a homogeneous  medium  of  aerosol  density  N.  then  the 
extinction  coefficient.  0,  for  the  laser  beam  is  given  by 


A-3)  /Jdirect  = N / 2ir  sm0  a (0)  d0 
Jo 

The  effective  extinction  coefficient  for  the  aura.  is 

tlUld 


A-4) 

^aura  ~ ^ Jq 

where 

f°< 

A-5) 

* ° . 

/ 

I < 

Jo 

7a<0)  (10 
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If  F0(z)  is  the  total  radiant  flux  in  the  direct  (ie  unscattered)  beam  at  a distance  z 
along  the  beam,  and  F|(z)  is  that  in  the  aura  (where  these  fluxes  have  been  integrated  over 
the  beam  and  aura  cross  sections)  then  the  equations  of  transfer  are 

A-6)  dF0 

~dz=0F° 

dF. 

from  which  we  obtain  for  the  aura  flux 

A-7)  F,  = PT  e^< 1 " *,/4)z  [ 1 - e-W1/4zJ 

where  Pj  is  the  initial  flux  in  the  direct  beam. 

Now  since  the  scattering  coefficient  for  the  aura  (ie  those  scatterings  which  do  not  de- 
plete the  aura)  is  given  by 

f0 c 

A-8)  N 2w  sine  oa(0)  d0  = /30V4«?C) 

we  may  rewrite  A-l ) in  the  form 
A-9)  r = fz3?2  ec  01/V/%C) 

where  is  of  order  unity.  Then,  the  power  Pr  received  by  a detector  of  area  A located 
within  the  aura  at  a distance  z along  the  beam  is 

F.-A  PtA*'  e-^1-^)2  [1-e-^zj 

A- 10)  Pr  = — 1 T-=— ' - 

?rr“  vz3  P0c2<l>(0c)f2 

where  j'  is  of  order  unity  and  depends  on  the  position  of  the  detector  within  the  aura,  the 
field  of  view,  detector  orientation,  etc. 

It  remains  to  estimate  <pc  and,  from  the  scattering  law,  0(0C).  To  do  this  we  consider 
a beam  launched  from  the  Earth’s  surface  with  elevation  angle  0e,  and  we  neglect  the  Earth's 
curvature.  C learly  0c~0e  since  a downward  scattering  at  an  angle  >0e  will  intersect  the  sur- 
face. Since,  from  A-*)),  the  aura  radius  at  a given  range  increases  as  0C  increases,  if  the  aura 
is  to  intersect  the  Earth’s  surface  at  a given  range  D,  then  the  beam  must  have  an  elevation 
angle  0e  = 0C  satisfying  the  relationship  D0C  = r or: 

A-  II ) jD3!2  ec  pH  <pVH0c)  = D 0C 


that  is 


Equation  A- 1 2)  determines  0C  once  the  functionaJ  form  of  0(0)  has  been  determined 
from  the  scattering  law.  To  proceed  further,  we  idealize  the  aerosol  Mie  scattering  cross 
section  as 


A-13)  o(0)  constant  x ( 1 - (“—)**)  0 < 0 < 0O 

Ur 


0 >0r 


Then  from  A-5) 

n A ^ ^ C ^ 7 0r* 

A- 14)  0<0)=s(2— ^)  (-i)2  n.  -i-^n] 
n c/q  n » tm  Uq 


If  the  range  D is  large. 


A- 15)  (?)2S-i 


n + 2 


ie  D0»1  than  from  A- 1 2)  and  A*14)  we  find 


Finally,  introducing  A-15)  and  A-12)  into  A-IO)  we  obtain  the  result 


A- 16) 


PR  e-PP B 

PtA  jtD 


f'f2  [eW2-l] 

1 - ^ 

( — 21 — ) 

L n + 2'  J 

n + 2 

n + 2 


To  estimate  f,  we  note  that  in  the  Amush4  calculation 
A- 17)  <r2>  = y pz3  <0> 

In  the  model  A-13),  <0>  ^ 2/3  0C  from  which  we  estimate  f ^ 2/3. 

From  Deirmendjian’s5  Mie  scattering  calculations  for  marine  haze  models,  we 
deduce  n **  1 and  0Q  ranging  from  about  9°  in  the  visible  (4500 A)  to  about  30°  at  1.19 
microns. 

In  order  to  compare  A- 18)  with  the  single  scattering  results  from  the  text,  we  must 
relate  0O  to  the  quantity  f(o)  - the  gain  over  isotropic  in  the  forward  direction.  Since 


A- 18) 


f(o)  = 2f2i 

4ir 


and 

n 


V-19)  [ 

-'n 


2jt  sin  0 d0  = 1 
4jt 


117 


then  from  A- 13)  for  n = 1 and  small  0 


A- 20)  2wi(o) 


)d0  = 1 


or 


A-21) 


nf(o)0o 2 
3 


= 1 


Then  since  the  first  bracketed  term  on  the  right  hand  side  of  A-16)  is  .ibout  I0(for 
f - 1 , j = 2/3,  n - I ) we  arrive  finally  at  the  following  estimate  for  large  D 


A-22) 


PR  _ lO0e-0D 
PjAf(0)  3D 


which  can  now  be  compared  to  the  results  for  T shown  in  the  text.  This  is  carried  out  in 
figure  A-l  where  we  compare  A-22)  with  0 = . 1 to  the  results  from  figure  6 in  the  text. 

Two  features  of  this  comparison  are  noteworthy.  First  it  is  seen  that  at  large 
distances  the  single  scattering  model  will  dominate  by  virtue  of  its  more  gradual  fall-off 
with  range.  When  the  extinction  coefficient  is  . I km-*,  the  cross  over  point  is  between 
130  and  1 50  km  depending  on  wavelength.  Under  poor  visibility  conditions,  this  cross 
over  point  will  shift  to  shorter  ranges.  At  0 = .257  km-*  (visibility  = 13  km  at  6000  A), 
for  example,  the  single  scattering  will  dominate  at  ranges  greater  than  about  75  km. 


Figure  A-l.  Comparison  of  Single  and  Multiple  Scattering. 
Extinction  Coefficient  * 0.1  KM"* 

I 18 


On  the  other  hand,  however,  it  is  significant  that  at  shorter  ranges,  the  aura  can 
exceed  the  single  scattered  component  by  about  10  dB.  Thus  if  conditions  are  right  (speci- 
fically that  the  vertical  beam  divergence  be  such  that  energy  is  radiated  at  the  critical  eleva- 
tion angle,  0C , determined  by  A- 1 2)  then  it  is  expected  that  at  the  shorter  ranges  (nominally: 
under  100  km),  a bright  point-like  source  due  to  the  aura  should  be  superimposed  or.  the 
line  source  due  to  single  scattering.* 

Clearly,  if  communications  system  design  is  to  be  based  upon  detection  of  the  aura, 
a more  detailed  and  rigorous  treatment  of  the  multiple  scattering  problem  - perhaps  via 
Monte  Carlo  simulation  would  be  required.  The  approach  adopted  in  the  present  study 
is  the  more  conservative  one  of  basing  detection  on  the  single  scattered  component  whose 
transmission  properties  can  be  modeled  in  a relatively  straightforward  manner. 


i 


i 


‘Indeed,  this  is  precisely  what  Curcio  and  Drumnieter'  observed  at  a 45  kin  range  when  the  meteorological 
range  was  20-25  km. 


i iy 
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